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Inhomogeneous superconductivity arises when the species participating in the pairing phenomenon 
have different Fermi surfaces with a large enough separation. In these conditions it could be more 
favorable for each of the pairing fermions to stay close to its Fermi surface and, differently from 
the usual BCS state, for the Cooper pair to have a non zero total momentum. For this reason in 
this state the gap varies in space, the ground state is inhomogeneous and a crystalline structure 
might be formed. This situation was considered for the first time by Fulde, Ferrell, Larkin and 
Ovchinnikov, and the corresponding state is called LOFF. The spontaneous breaking of the space 
symmetries in the vacuum state is a characteristic feature of this phase and is associated to the 
presence of long wave-length excitations of zero mass. The situation described here is of interest 
both in solid state and in elementary particle physics, in particular in Quantum Chromo-Dynamics 
at high density and small temperature. In this review we present the theoretical approach to the 
LOFF state and its phenomenological applications using the language of the effective field theories. 
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I. INTRODUCTION 

Superconductivity is one of the most fascinating chapters of modern physics. It has been a continuous source of 
inspiration for different realms of physics and has shown a tremendous capacity of cross-fertilization, to say nothing 
of its numerous technological applications. This review is devoted to a less known chapter of its history, i.e. inho- 
mogeneous superconductivity, which arises when the main property of the superconductor is not uniform in space. 
Before giving a more accurate definition of this phenomenon let us however briefly sketch the historical path leading 
to it. Two were the main steps in the discovery of superconductivity. The former was due to Kamcrlingh Onnes 
l)Kamerlingh Onneslll91l|) who discovered that the electrical resistance of various metals, e. g. mercury, lead, tin and 
many others, disappeared when the temperature was lowered below some critical value T c . The actual values of T c 
varied with the metal, but they were all of the order of a fe w K, or at most of the ord e r of te nths of a K. Subsequently 
perfect diamagnetism in superconductors was discovered ijMeissner and Ochsenfeldl . Il933|) . This property not only 
implies that magnetic fields are excluded from superconductors, but also that any field originally present in the metal 
is expelled from it when lowering the temperature belo w its critical value. These tw o features were captured in the 
equations proposed by the brothers F. and H. London ({London and Londonl Il935|) who first realized the quantum 
character of the phenomenon. The decade starting in 1950 was the stage of two major theoretical breakthroughs. 
First, Ginzburg and Landau (GL) created a theory describing the transition between the superconducting and the 
normal phases ijGinzburg and LandauL 17950'!. It can be noted that, when it appeared, the GL theory looked rather 
phcnomenological and was not really appreciated in the western liter ature. Seven years la ter Bardeen, Cooper and 
Schrieffer (BCS) created the microscopic the ory that bears their name ijBardeen et a^J . ll957j) . Their theory was based 
on the fundamental theorem llCoopeii 1195^1 . which states that, for a system of many electrons at small T, any weak 
attraction, no matter how small it is, can bind two electrons together, forming the so called Cooper pair. Subse- 
quently in l)Gor'kovL Il959|) it was realized that the GL theory was equivalent to the BCS theory around the critical 
point, and this result vindicated the GL theory as a masterpiece in physics. Furthermore Gor'kov proved that the 
fundamental quantities of the two theories, i.e. the BCS parameter gap A and the GL wavefunction if), were related 
by a proportionality constant and if> can be thought of as the Cooper pair wavefunction in the center-of-mass frame. 
In a sense, the GL theory was the prototype of the modern effective theories; in spite of its limitation to the phase 
transition it has a larger field of application, as shown for example by its use in the inhomogeneous cases, when the 
gap is not unifor m in space. Anot her remarkable advance in these years was the Abrikosov's theory of the type II 
superconductors l)Ahrikosovl Il957|) . a class of superconductors allowing a penetration of the magnetic field, within 
certain critical values. 

The inspiring p ower of superconductivity became s oon evident in the field of elementary particle physics. Two 
pioneering papers ijNambu and Jona-Lasiniol Il961albl) introduced the idea of generating elementary particle masses 
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through the mechanism of dynamical symmetry breaking suggested by superconductivity. This idea was so fruitful 
that it eventually was a crucial ingredient of the Standard Model (SM) of the elementary particles, where the masses 
are generated by the formation of the Higgs condensate much in the same way as superconductivity originates from 
the presence of a gap. Furthermore, the Meissner effect, which is characterized by a penetration length, is the origin, 
in the elementary particle physics language, of the masses of the gauge vector bosons. These masses are nothing but 
the inverse of the penetration length. 

With the advent of QCD it was ear ly realized that at high density, due to the asymptotic freedom property 
(G ross and Wilczekl 11973: Politzer, 197 3]) and to the existence of an attractive channel in the color interaction, diquark 
condensates might be formed ijBailin and LoveL 1 19841 iBarroisl Il977t ICollins and Perrvt Il975l iFrautschl Il978j) . Since 
these condensates break the color gauge symmetry, the subject took the name of color superconductivity. H owever , 
only i n the last few years this has become a very active field of research; these developments are reviewed in (|Alfordl 
120011 iHond. l200ll IHsuI l2000t iNardullil l2002al iRaiaeopal and Wilczekl l200l|) . It should also be noted that color 
superconductivity might have implications in astrophysics because for some compact stars, e.g. pulsars, the baryon 
densities necessary for color superconductivity can probably be reached. 

Superconductivity in metals was the stage of another breakthrough in the 1980s with the discovery of high T c 
superconductors. As we anticipated, however, the main subject of this review is a different and separate development 
of superconductivity, which took place in 1964. It originates in high-field superconductors where a strong magnetic 
field, coupled to the spins of the conduction electrons, gives rise to a separation of the Fermi surfaces corresponding to 
electrons with opposite spins. If the separation is too high the pairing is destroyed and there is a transition (first-order 
at small temperature) from the super conducting state to the normal one. In two separate and contemporary papers, 
l|Larkin and Ovchinnikovl[l964|) and l|Fulde and Ferrell Il964|) . it was shown that a new state could be formed, close 
to the transition line. This state that hereafter will be called LOFF 1 has the feature of exhibiting an order parameter, 
or a gap, which is not a constant, but has a space variation whose typical wavelength is of the order of the inverse 
of the difference in the Fermi energies of the pairing electrons. The space modulation of the gap arises because the 
electron pair has non zero total momentum and it is a rather peculiar phenomenon that leads to the possibility of a non 
uniform or anisotropic ground state, breaking translational and rotational symmetries. It has been also conjectured 
that the typical inhomogeneous ground state might have a periodic or, in other words, a crystalline structure. For 
this reason other names of this phenomenon are inhomogeneous or anisotropic or crystalline superconductivity. 

Inhomogeneous superconductivity in metals has been the object of intense experimental investigations especially 
in the last decade; for reasons to be discussed below the experimental research has aimed to rather unconventional 
superconductors, such as heavy fermion superconductors, quasi-two dimensional layered organic superconductors 
or high T c superconductors. While different from the original LOFF proposal, these investigations still aim to a 
superconducting state characterized by non zero total momentum of the Cooper pair and space modulation of its 
wavefunction. At the moment they represent the main possibility to discover the LOFF state in condensed matter 
physics. 

Quite recently it has been also realized that at moderate density the mass difference between the strange and the 
up and down quarks at the weak equilibrium and/or color and electric neutrality lead t o a difference in the Fermi 
momenta, which renders in principle the LOFF state possible in color superconductivity l)Alford et al l l2001tJ> . The 
same authors have pointed out that this phenomenon might have some relevance in explaining the sudden variations 
of the rotation period of the pulsars (glitches). 

The main aim of this review is to present ideas and methods of the two main roads to inhomogeneous supercon- 
ductivity, i.e. the condensed matter and the QCD ways. Our approach will be mainly theoretical and the discussion 
of phenomenological consequences will be limited, first because we lack the necessary skills and second because the 
theory of the LOFF superconductivity is up to now much more advanced than experiment and its main phenomeno- 
logical implications belong to the future. For this reason we will give large room to the theoretical foundations of 
inhomogeneous superconductivity and will present only a summary of experimental researches. Our scope is to show 
the similarities of different physical situations and to present a formalism as unified as possible. This not only to 
prove once again the cross-fertilization power of superconductivity, but also to expose experts in the two fields to 
results that may be easily transferable from one sector to the other. Moreover, by presenting the LOFF phenomenon 
in a unified formalism, this review can contribute, we hope, to establish a common language. To this end we discuss 
the LOFF state both in solid state and in QCD physics starting with Nambu Gor'kov (NG) equations. For the solid 
state part they will be derived by the effective theory of the relevant degrees of freedom at the Fermi surface and in 
the QCD sector by the so called High Density Effective Theory (HDET) that, as we shall see, leads to equations of 
motion which coincide with the NG equations. In this way one is able to get in touch with a dictionary allowing to 



1 In the literature the LOFF state is also known as the FFLO state. 
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switch easily from one field to the other. 

The pla n of this review is as follow s. In Section IT T ] we start descr ibing the general formalism, based on NG 
equations ijrTor'kovilTflR^lNambiiLFmBnl) . As shown by (Polchinskil Fmfll using the Renormalization Group approach, 
the excitations at the Fermi surface can be described by an effective field theory. Its equations of motion are exactly 
the NG equations of ordinary (homogeneous) superconductivity. We will then apply this formalism to fermions 
with different Fermi surfaces. The difference can be due to a magnetic field producing an energy splitting between 
spin up and spin down electrons, or, as in QCD, to a difference in the chemical potential originating from weak 
equilibrium, or color and electric neutrality, or mass difference between the pairing fermions. We will discuss the 
circumstances leading, in these cases, to inhomogeneous superconductivity. The Ginzburg Landau expansion can be 
used, as already mentioned, for the description of the inhomogeneous phase. It will be discussed in Section HTT1 both 
at zero temperature and close to the tricritical point. The T = case is more interesting for QCD applications while 
the finite temperature case might be relevant in condensed matter. In Section lTvl we will switch to QCD. We will first 
give a brief introduction to color superconductivity and then a description of the effective lagrangian for quarks at 
zero temperature near to the Fermi surface. We will also discuss more specifically the LOFF case for QCD with two 
massless flavors. Since in the LOFF phase both translational and rotational symmetries are spontaneously broken, the 
Goldstone theorem requires the presence in the physical spectrum of long wave-length, gapless, excitations (phonons). 
In Section [V] we discuss the phonon effective lagrangians for two crystalline structures, i.e. the single plane wave 
and the cubic structure. We will limit our presentation to the QCD case, though the presence of these excitations is 
obviously general. We will also discuss the gluon propagation inside these two crystalline media. In Section lvTl we will 
discuss the possible phenomenological applications of the LOFF phase. This discussion will go from strongly type II 
superconductor to two-dimensional structures for condensed matter. For hadronic matter we will discuss applications 
both in nuclear physics and in QCD, with particular emphasis on the physics of glitches in pulsars. 

Let us conclude this introduction by apologizing to the many authors whose work is not reviewed here in depth. 
Space limits forced us to sacrifice a more detailed exposition; the extensive bibliography at the end should help to 
excuse, we hope, this defect. 

II. THE GENERAL SETTING 

In this Section we give a pedagogical introduction to inhomogeneous superconductivity. We begin by reviewing 
homogeneous superconductivity by a field theory with effective Nambu-Gor'kov spin 1/2 fields describing quasi- 
particles. The effective field theory considers only the relevant degrees of freedom in the limit of small temperatures 
and high chemical potential; they are the modes in a shell around the Fermi surface. The dominant coupling in this 
limit is the four fermion interaction as first introduced in the BCS model. The dominance of this coupling can be also 
proved in a modern l anguage by using the renormalization group approach ijBenfatto and Gallavottlll990HPolch inski. 
ll99aiShankarl 11994ft . which shows that the BCS coupling is marginal and therefore, in absence of relevant couplings, 
it can dominate over other irrelevant couplings and produce the phenomenon of superconductivity. 

After having derived the Nambu-Gor'kov equations and the gap equation in Subsection 111. Al we discuss the case 
of homogeneous superconductor in Section III.BI and analyze its phase diagram in Section III. B. II We assume from 
the very beginning that the two species participating in the Cooper pairing have different chemical potentials, as this 
is the necessary situation for the LOFF state. In Section Til. CI we discuss the case of anisotropic superconductivity. 
In Section III. C. II we will show that for appropriate values of the difference in chemical potentials an anisotropic 
modulated gap A(r) oc exp(iq • r) leads to a state that is energetically favor ed in comparison to bot h the BCS and 
the normal non superconducting states. This was the state first discussed in l)Fulde and 

FerrellL 11964 - 

A. Nambu-Gor'kov equations 

To start with we consider, at T = 0, a fermion liquid formed by two species, that we call u and d, having 
different Fermi energies. In the electron superconductivity, as in the original LOFF papers ( F ulde and Ferreij Il964t 
lLarkin and Ovchinnikovll 1964ft . the species are the electron spin up and down states, but our formalism is general 
and will be applied later to the case where the fermion forming the Cooper pair are two quarks with different flavors. 
In superconducting materials the difference of chemical potentials can be produced by the presence of paramagnetic 
impurities. All these cases give rise to an effective exchange interaction that can described by adding the following 
term to the hamiltonian 



Hexch = -<5/z^V 3 ^. 



(2.1) 
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In the case of electron superconductivity 5/j, is proportional to the magnetic field and the effect of l|2.1|l is to change 
the chemical potentials of the two species: 



p u = p + Sp, p d = p- Sp 
Adopting a BCS interaction, the action can be written as follows 



.4 



A = A +A B cs 
dp 



dt 



(2tt) 3 



^(p) [id t - E(p) +n + Spa 3 ) i/j(p) 



(2.2) 

(2.3) 
(2.4) 



Abcs = f J dtH ^ tyt(pity(P4)) (V' t (P2MP3)) x (2^( Pl +p 2 -p 3 -p 4 ). (2.5) 

Here and below, unless explicitly stated, ip(p) denotes the 3D Fourier transform of the Pauli spinor ip(r,t), i.e. 
ip(p) = Vv(p, i). For non relativistic particles the functional dependence of the energy would be E(p) = p 2 /2m, but 
we prefer to leave it in the more general form (|2.4(l . 
The BCS interaction (|2.5|) can be written as follows 



Abcs — A con d + A int , 



(2.6) 



with 



Acond = / dt JJ 



dpk 



S(P3, P4)V- T (Pl)CV T (P2) 



S*(pi, P2)^(P3)C^(p 4 ) (2tt) 3 5(pi + p 2 - P3 - P4) , 



^(pi)C^(p 2 ) + s*( Pl ,p 2 ) 



X 



where C = ier 2 and 



■0(p 3 )C^(p 4 ) - 3(p 3 , p 4 ) (2tt) 3 S(pi + p 2 - p 3 - P<t) 



§(p, P ') =< v(p)c^(p') > 



(2.7) 



(2.8) 



In the mean field approximation the interaction term can be neglected while the gap term A conc i is added to Aq. Note 
that the spin condensate S(p, p') is simply related to the condensate wave function 



by the formula 



H(r) =<V(r,i)CV(r,t)> 



^-/(^(If^ p+pV ^-')- 



(2.9) 



(2.10) 



In general the condensate wavefunction can depend on r; only for homogeneous materials it does not depend on the 
space coordinates; therefore in this case S(p, p') is proportional to S(p + p'). 

In order to write down the Nambu-Gor'kov (NG) equations we define the NG spinor 



x(p) 



1 



( <Mp) 



where we have introduced the charge-conjugate field 



We also define 



V/ = cy> f . 



(2.11) 



(2.12) 



(2.13) 
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The free action can be therefore written as follows: 



with 



A = J dt (2nf ( 2 t)3 ^(p)S- l (p t p')x(p'), (2-14) 

5- 1 (p,p / ) = (2 7 r) 3 f (i5t -^+?rt (P_P ' ) r-p^,~A (P,P l mV (2-15) 

V ' V -A*(p,p') (zdj +£ p + £/ZO- 3 )<5(p-p')/ v y 

Here 

£ p = E(p) - p « v F • (p - Pf .) , (2.16) 

where 

fl£(p) 



<9p 



(2.17) 



is the Fermi velocity. We have used the fact that we are considering only degrees of freedom near the Fermi surface, 
i.e. 

p F -S<p<p F + 8, (2.18) 
where 6 is the ultraviolet cutoff, of the order of the Debye frequency. In particular in the non relativistic case 

£ =£^_Zi ; Vf = — (2.19) 
p 2m 2m ' m 

S 1-1 in i|2.15|l is the 3D Fourier transform of the inverse propagator. We can make explicit the energy dependence by 
Fourier transforming the time variable as well. In this way we get for the inverse propagator written as an operator: 

^=(^r _~f yl ), (2.20) 

and 

[G+p 1 = E-£p+fy«73+iesignJ5, 

[Go]- 1 = -E-fr -6(mt 3 - design E, (2.21) 

with e = + and P the momentum operator. The ie prescription is nothing but the usual one for the Feynman 
propagator, that is forward propagation in time for the energy positive solutions and backward propagation for the 
negative energy solutions. As for the NG propagator S, one gets 

S=(% f) ■ (2.22) 

S has both spin, cr, a', and a, b NG indices, i.e. S^, 2 . The NG equations in compact form are 

S~ 1 S = 1, (2.23) 

or, explicitly, 

[G+]- 1 G + AF = 1 , 

-[G ]~ 1 F + A*G = 0. (2.24) 



2 We note that the presence of the factor l/\/2 in 12.111 implies an extra factor of 2 in the propagator: S(x, x') = 2 < T {x( x )x^ ( x> ) } >> 
as it can be seen considering e.g. the matrix element < 7 |v(a.')v^(a:')| >= (idt — — Sfiaz) 5(x — x'), with (x = (t,r)). 
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Note that we will use 

< r | A|r ' >= | 3(r) *(r - r') = A(r) 5{r - r') , (2.25) 

or 

<p|A|p'>=A(p !P ') (2.26) 
depending on our choice of the coordinate or momenta representation. The formal solution of the system (|2.24l) is 

F = G p A*G , 

G = G+ - G+AF , (2.27) 



so that F satisfies the equation 



and is therefore given by 



F = G A*(G+-G+AF) (2.28) 



(2.29) 



A*[G+]-i[A*]-i[G ]-i + A* A ' 

In the configuration space the NG Eqs. (|2.24|) are as follows 

(E-E(-iV) + fj, + Sna 3 )G(r,r',E)+ A(r)f(r,r',J5) = 6(r - r') , 
(-E-E(-iV) + v-6n<73)F(r,r',E)-A*(r)G(r,r\E) = 0. (2.30) 

The gap equation at T = is the following consistency condition 

A*(r)=-ilf^TrF(r,r,E), (2.31) 
where F is given by (|2.29(l . To derive the gap equation we observe that 



. 9 
l 2 



which gives l|2.31|l . 

At finite temperature, introducing the Matsubara frequencies u„ = (2n + 1)ttT, the gap equation reads 



A*(r) = 9 -T J2 TrF(r,r,E) _ . (2.33) 
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B. Homogeneous superconductors 

It is useful to specialize these relations to the case of homogeneous materials. In this case we have 

2A 

H(r) = const. = . (2.34) 

9 

S(pi,p 2 ) - — -^(2^( Pl +p 2 ) . (2.35) 
9 PpO 
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Therefore one gets 

A( Pl ,p 2 )=A<5( Pl -p 2 ) (2.36) 

and from and TH^ 

A(r) = A*(r) = A . (2.37) 
Therefore F(r, r, E) is independent of r and, from Eq. 12.29|) . one gets 

TrF(r, r,E) = -2 A / ^ ( - _ ^ g - A2 (2.38) 

which gives the gap equation at T = 0: 



. dE d 3 p 1 . , 

'^/^(2^ ^-^-a-A^ (2 - 39) 



and at T ^ 0: 



y MR 



d 3 p A 



+ ^) 2 + e(p,A) 2 



with 



We now use the identity 



1 + °° 1 
-[l-n u -n d ] = e(p,A)T £ (2.42) 

2 K, + i(5^) 2 + e 2 (p,A) 



where 



(2.40) 



e(p,A) = ,/A2+a. (2.41) 



n " (p) = e (e+j/T + 1 > ^(P = e ( £ -j/T + 1 • ( 2 - 43 ) 

The gap equation can be therefore written as 

In the Landau theory of the Fermi liquid n u , rid are interpreted as the equilibrium distributions for the quasiparticles 
of type u, d. It can be noted that the last two terms act as blocking factors, reducing the phase space, and producing 
eventually A — > when T reaches a critical value T c (see below). 

Before considering the solutions of the gap equations in the general case let us first consider the case S^i — 0; the 
corresponding gap is denoted Ao- At T = there is no reduction of the phase space and the gap equation becomes 

<««) 



(2.46) 



2 J (2tt)3 e(p,A ) 
whose solution is (we have assumed d 3 p = p 2 F dpdVl) 

5 



sinh — 

gp 



Here 
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is the density of states and we have used £ p w vp(p — pp), see Eqs. (I2.16ll - H2.19ll . In the weak coupling limit (|2.46(l 
gives 

A = 2Se- 2/p9 . (2.48) 

Let us now consider the case S/i ^ 0. By (|2.44|) the gap equation is written as 

g f d 3 p 1 _ g f d 3 p n u + n d 
- 1+ 2j WT-e-2j (2^)3 e ' (2 ' 49) 

Using the gap equation for the BCS superconductor, the l.h.s can be written, in the weak coupling limit, as 

Lh.s^in^E, (2.50) 
2 A y J 

where we got rid of the cutoff 8 by using A , the gap at Sfi — and T = 0. Let us now evaluate the r.h.s. at T = 0. 
We get 



r.h.s. 



9P f S d£ p 



= f / Sft) + 6(-e + dp)} . (2.51) 

t=o 2 J n e 



The gap equation at T = can therefore be written as follows: 



In ^ = 6{5pi - A)arcsinh ^ 2 A A " , (2.52) 



A 
i.e. 

In A ° = 0. (2.53) 

+ y/Sp 2 - A 2 

One can immediately see that there are no solutions for Sfi > Aq. For S pi < Ao one has two solutions. 

a) A - Ao, (2.54) 
6) A 2 = 2 5^0-/^1. (2-55) 

The first arises since for A = Ao the l.h.s. of the Eq. (|2.52fl is zero. But since we may have solutions only for S/i < Ao 
the ^-function in Eq. I|2.52|l makes zero also the r.h.s.. The existence of this solution can also be seen from Eq. H2.39fl . 
In fact in this equation one can shift the integration variable as follows: E — > E + <5/i, getting the result that, in the 
superconductive phase, the gap A is independent of 8 pi, i.e. A = Ao. 

To compute the free energy we make use of the theorem saying that f or small variations of an ext ernal parameter 
of the system all the thermodynamical quantities vary in the same way ijLandau and Lifshitd . Il996j) . We apply this 
to the grand potential to get 

f e -0- <-> 

From the expression of the interaction hamiltonian (see Eq. (|2.5|) ~) we get immediately (cfr. l)Ahrikosov et aZl Il963|) . 
cap. 7): 

n = - \% i ^l A (^)! 2 ■ (2-57) 



For homogeneous media this gives 



Using the result (|2.48|) one can trade the integration over the coupling constant g for an integration over A , the BCS 
gap at 8fi = 0, because dA /A = 2dg/pg 2 . Therefore the difference in free energy between the superconductor and 
the normal state is (we will use indifferently the symbol Q for the grand potential and its density fl/V) 



P f A ° a2 dA 



- = - ^ / A — — . (2.59) 

Z JA f ^0 
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Here Af is the value of Ao corresponding to A = 0. A/ = in the case a) of Eq. I|2.54[l and A/ = 2(5 in the case 
b) of Eq. I|2.55|) : in the latter case one sees immediately that J7a — f^o > because from Eq. 12.55fl it follows that 
Ao < 2(5 /j. The free energies for 5/j, ^ corresponding to the cases a), b) above can be computed substituting (|2.54(l 
and H2.55fl in l|2.59|l . Before doing that let us derive the density of free energy at T = and (5/j. 7^ in the normal non 
superconducting state. Let us start from the very definition of the grand potential for free spin 1/2 particles 



(2.60) 



Q o (0,T) = -2VT J ^m(l + e^P»/ T ) . 
Integrating by parts this expression we get, for T — ► 0, 

n °(°) = I dn pP 3 de9{fi-e). (2.61) 

127T° J 

From this expression we can easily evaluate the grand-potential for two fermions with different chemical potentials 
expanding at the first non-trivial order in <5/i/ /i. The result is 

Su 2 

n (Sn) = n (o)--|-p. (2.62) 

Therefore from 12.54|l . (|2.55jl and Ij2.59jl in the cases a), b) one has 

a) n A {Sti) = n {5^)-^{-2S^ + A 2 ), (2.63) 

b) n A (S^) = Clo(Sn)-^(-4S f j, 2 + 48fiA -Al). (2.64) 

Comparing (|2.63|) and (|2.64|) we see th at the solution a) has lower fl. Therefore, for Sfi < Aq/\/2 the BCS supercon- 
ductive state is stable llClogstonllT96l . At <5/z = A /\/2 it becomes metastable, as the normal state has a lower free 
energy. This transition would be first order since the gap does not depend on 5/i. 

The grand potentials for the two cases a) and b) and for the gapless phase, Eq. (|2.62l) . are depicted in Fig. ^ 
together with the corresponding gaps. 



8|i 

Ao 



A 2 

A2 1 " 4 / 

A 2 / 
1.0 £ 

/ 

/ 

/ 

0.6 / 

/ 

/ 

0.2 / 

^_ / 

6^--0^4 0.6 0.8 1.0 1.2 

-0.2 

-0.4 V \ 

\ 

\ 

a -0.6 \ 

FIG. 1 Gap and grand potential as functions of Sfj, for the two solutions a) and b) discussed in the text, see Eqs. \2. 5J$ , \2. 55)) 
and 12.6$ , \2. &4t - Upper solid (resp. dashed) line: Gap for solution a) (resp. solution b)). In the lower part we plot the grand 
potential for the solution a) (solid line) and solution b) (dashed line); we also plot the grand potential for the normal gapless 
state with S/j, ^ (dashed- dotted line). All the grand potentials are referred to the value fio(0) (normal state with Sfi — 0). 
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A different proof is obtained integrating the gap equation written in the form 

- = (2.65) 

The normalization can be obtained considering the homogeneous case with Sfi — 0, when, in the weak coupling limit, 
from Eqs. Ij2.57jl and Ij2.48jl one gets 

= -|Ao, (2-66) 

see below Eq. (I2.71f) . In this way one obtains again the results (|2.63[1 and 12.64f> . 

This analysis shows that at 8/i = Sfii = An/v^ one would pass abruptly from the superconducting (A ^ 0) to the 
normal (A = 0) phase. However, as we shall discuss below, the real ground state for Sfj, > Sfi\ turns out to be an 
inhomogeneous one, where the assumption (|2.37|l of a uniform gap is not justified. 

1. Phase diagram of homogeneous superconductors 

We will now study the phase diagram of the homogeneous superconductor for small values of the gap parameter, 
which allows to perform a Ginzburg-Landau expansion of gap equation and grand potential. In order to perform a 
complete study we need to expand the grand-potential up to the 6 th order in the gap. As a matter of fact in the 
plane (Sfi,T) there is a first order transition at (<5/ii,0) and a second order one at (0,T C ) (the usual BCS second 
order transition). Therefore we expect that a second order and a first order lines start from these points and meet 
at a tricritical point, which by definition is the meeting point of a second order and a first order transition line. A 
tricritical point is characterized by the simultaneous vanishing of the A 2 and A 4 coefficients in the grand-potential 
expansion, which is why one needs to introduce in the grand potential the 6 th order term. For stability reasons the 
corresponding coefficient should be positive; if not, one should include also the A 8 term. 

We consider the grand potential, as measured from the normal state, near a second order phase transition 

n = laA 2 + i/3A 4 + i 7 A 6 . (2.67) 
2 4 6 

Minimization gives the gap equation: 

«A + /3A 3 +7A 5 = 0. (2.68) 

Expanding Eq. (|2.4l)|) up to the 5 th order in A and comparing with the previous equation one determines the 
coefficients a, f3 and 7 up to a normalization constant. One gets 



^ s 

A = 2gpTRe V / d£_ 



n=0 



A A 3 



(2.69) 



.^i+e) ^i+er w+er 

with 

uj n = uj n + iSp = (2n + l)7rT + i6[i . (2.70) 

The grand potential can be obtained, up to a normalization factor, integrating in A the gap equation. The normal- 
ization can be obtained by the simple BCS case, considering the grand potential as obtained, in the weak coupling 
limit, from Eqs. l(237jl and 

n = -fAj». (2.71) 

The same result can be obtained multiplying the gap equation H2.45J) by An and integrating the result provided that 
we multiply it by the factor 2/g, which fixes the normalization. Therefore 



2 

a = — 
9 



p = 4pTRe^2 / ,- 2 ' 2 , 2 , (2-73) 

7 - - 4prKe ?J W+P?' (2 - 74) 
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In the coefficients (3 and 7 we have extended the integration in £ up to infinity since both the sum and the integral 
are convergent. To evaluate a is less trivial. One can proceed in two different ways. One can sum over the Matsubara 
frequencies and then integrate over £ or one can perform the operations in the inverse order. Let us begin with the 
former method. We get 



2 

a = — 
9 



I .--'dlL j |(tanh 



2T 



tanh 



2T 



(2.75) 



Performing an integration by part we can extract the logarithmic divergence in S. This can be eliminated using the 
result i|2.46[l valid for <5/i = T = in the weak coupling limit 



99, 25 
~2~ 



(2.76) 



We find 



a = p 



log- 



2T 1 



efcc lna; 



COSh 2 ^ COSh 2 



(2.77) 



where 



<5/i 



(2.78) 



Defining 



log 



2T c (y) 



(ix lna; 



cosh 2 cosh 2 5^ 



(2.79) 



we get 



a(v,t) = plog 



tc(«/t) ' 



(2.80) 



where 



t = 



tr. = 



(2.81) 



Therefore the equation 



* = t c (v/t) 



(2.82) 



defines the line of the second ord er phase transition. Perform ing the calculation in the reverse order brings to a 
more manageable result for t c {y) ijBuzdin and Kachkachil Il997[) . In Eq. I|2.72|l we first integrate over £ obtaining a 
divergent series which can be regulated cutting the sum at a maximal value of n determined by 



uj n = S => N 



2ttT 



(2.83) 



We obtain 



a = - I 1 — n g p 
9 \ 



TReY —) 



(2.84) 



The sum can be performed in terms of the Euler's function ip{z) 

N 



x ^ 1 1 

Re ? = 7T- 

^ i7i- O.tt 



n=0 



2ttT 



Re 



1.1/ 
2 +Z 27 



2vrT 



log 



2vrT 



Reip 



1.1/ 
2 +Z 27 



(2.85) 
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Eliminating the cutoff as we did before we get 



By comparing with Eq. (|2.77|) we get the following identity 



^ +i ^) = -^ )+1 ^r^^^ + ^^) ■ (2 - 87) 



The equation (|2.79|) can be re-written as 



47rT c (y) Y \2 2tt 

In particular at Sfi = 0, using (C the Euler-Mascheroni constant 



we find from Eq. 



-l,»j.|4-.i. -,=<'•. f U..-777.... (2.80) 



7rT 

a(0,T/A ) = log — , (2.90) 
7A0 



reproducing the critical temperature for the BCS case 

T c = - A ~ 0.56693 A . (2.91) 

The other terms in the expansion of the gap equation are easily evaluated integrating over £ and summing over the 
Matsubara frequencies. We get 

- „ rJb f;£_-^,W»(! +< i) , (2.92) 

71=0 ™ x 7 



n=0 

where 

^ w = d^w- (2 - 94) 

Let us now briefly review some results on the grand potential in the GL expansion l|2.67|l . We will assume 7 > in 
order to ensure the stability of the potential. The minimization leads to the solutions 

A = 0, (2.95) 



A 2 = Ai 

27 



- _L (-f3 ± V/3 2 - 4a 7 ) . (2.96) 



The discussion of the minima of depends on the signs of the parameters a and (3. The results are the following: 



1. a > 0, (3 > 



In this case there is a single minimum given by (j2.95(l and the phase is symmetric. 
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2. a > 0, (5 < 



Here there are three minima, one is given by (|2.95() and the other two are degenerate minima at 

A = ±A+ . (2.97) 
The line along which the three minima become equal is given by: 




n(0) = fi(±A+) — > p=-AJ-^-. (2.98) 

V ^ 

Along this line there is a first order transition with a discontinuity in the gap given by 

To the right of the first order line we have f2(0) < f2(±A+). It follows that to the right of this line there is the 
symmetric phase, whereas the broken phase is in the left part (see Fig. 0). 



3. 



a < 0, P > 



In this case Eq. 12.95fl gives a maximum, and there are two degenerate minima given by Eq. H2.97fl . 
Since for a > the two minima disappear, it follows that there is a second order phase transition along the line 
a = 0. This can also be seen by noticing that going from the broken phase to the symmetric one we have 

lim A^ = 0. (2.100) 

ct—>0 



4. 



a < 0, P < 



The minima and the maximum are as in the previous case. 

Notice also that the solutions A± do not exist in the region P 2 < icry. The situation is summarized in Fig. [21 Here we 
show the behavior of the grand potential in the different sectors of the plane (a/7, P/j), together with the transition 
lines. Notice that in the quadrant (a > 0, p < 0) there are metastable phases corresponding to non absolute minima. 
In the sector included between the line P = —2^fotpy and the first order transition line the metastable phase is the 
broken one, whereas in the region between the first order and the a = lines the metastable phase is the symmetric 
one. 

Using Eqs. H2.86(l . H2.92[l and (|2.93|) which give the parameters a, P and 7 in terms of the variables v = Sfx/Ao 
and t = T/Ao, we can map the plane a and p into the plane (<5^t/Ao, T/Aq). The result is shown in Fig. [3] From 
this mapping we can draw several conclusions. First of all the region where the previous discussion in terms of the 
parameters a, P and 7 applies is the inner region of the triangular part delimited by the lines 7 = 0. In fact, as already 
stressed, our expansion does not hold outside this region. This statement can be made quantitative by noticing that 
along the first order transition line the gap increases when going away from the tricritical point as 

(2.101) 

Notice that the lines P(v,t) = and j(v, t) — are straight lines, since these zeroes are determined by the functions 
tfj^ and i/^ 4 ) which depend only on the ratio v/t. Calculating the first order line around the tricritical point one gets 
the result plotted as a solid line in Fig. [3] Since we know that S/j = S/ii = is a first order transition point, the 

first order line must end there. In Fig. |21we have simply connected the line with the point with grey dashed line. To 
get this line a numerical evaluation at all orders i n A would be required. This is feasible but we will skip it since the 
results will not be necessary in the following, see l)Sarmalll963l). The location of the tr i critical point is determined b y 
the intersection of the lines a = and P = 0. One finds ijBuzdin and Kachkachiri997tlCombescot and Moral 12002) 




()// 
A~o 



= 0.60822, — 

trie Z\q 



trie 



0.31833. (2.102) 



We also note that the line a = should cross the temperature axis at the BCS point. In this way one reobtains the 
result in Eq. (|2.91l) for the BCS critical temperature, and also the value for the tricritical temperature 

-^^= 0.56149. (2.103) 

Tbcs 
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FIG. 2 The graph shows the first order and the second order transition lines for the potential of Eq. \2. 61\j . We show the 
tricritical point and the regions corresponding to the symmetric and the broken phase. Also shown is the behavior of the grand 
potential in the various regions. The thin solid line is the locus of the points f3 2 — 4cey = 0. In the interior region we have 
/3 2 - 4a7 < 0. 

The results given in this Section are valid as long as other possible condensates are neglected. In fact, we will see 
that close to the first order transition of the homogeneous phase the LOFF phase with inhomogeneous gap can be 
formed. 

C. Gap equation for anisotropic superconductor: One plane wave (FF state) 

Let us now consider again the condensate wave function S(r) of Eq. (|2.9|) : 

3(r) = (vac\i>(r,t)Cip(r,t)\vac) . (2.104) 
Here \vac) is the ground state. We develop it as follows 

oo 

\vac) = c N \N) , (2.105) 

where N is even, the state \N) contains N/2 quark pairs of momenta 

Pi=+P + q, P2 = -p + q, (2.106) 

respectively for up and down species and the sum also implies an integration over the p variables and sum over spin. 
Clearly we have 

S(r) = Y, c *N c M{N\^{v,t)C^{v,t)\M) 

N,M 

= J2 c *NCN + 2(N\ij(r,t)Cij(r,t)\N + 2) = 

N 

= Y, c *NCN+2e 2lclN - r (N\yj(0)Cyj(0)\N + 2) . (2.107) 

N 
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FIG. 3 The curve shows the points solutions of the equation A = in the plane (v,t) — (5/i/Ao, T/Aq). The tricritical point 
at (Sfj,,T) « (0.62, 0.28) Ao is also shown. The upper part of the curve (solid line) separates the homogeneous phase from the 
normal one. Along the dashed line A = but this is not the absolute minimum of the grand potential. 



The homogeneous solution discussed in the previous subsection corresponds to the choice (Cooper pairs) 

q N = (for all N) , (2.108) 

while ^ corresponds to the inhomogeneous state. Let us now assume that the interaction favors the formation 
of pairs with non zero total momentum and suppose that the values qi, q2, ...qp are possible. Clearly this hypothesis 
has to be tested by comparing the values of the free energies for the normal, homogeneous and non homogeneous 
state. In any event, under such hypothesis, since the gap is proportional to S(r), one would get 

p 

A(r) = A m e 2 "i'"- r . (2.109) 

m— 1 

We will call the phase with A(r) given by l|2.109|l inhomogeneous or LOFF superconducting. At the moment we shall 
assume the existence of a single q and therefore 

A(r) = Ae 2lqr . (2.110) 

This is the simple st hypothesis, the one considere d in ( Fuld e and Ferrel1lll964[K see Fig. H it is therefore called FF 
state. The paper (jLarkin and Ovchinnikovl Il964l) examines the more general case (I2.109H ; we will come to it below. 
The assumption (|2.106|l with q ^ produces a shift in energy: 

£ p ± 5fi = v F ■ (p - p F ) ± Sfx -> v F • (p =f q - Pf) ± <5m = C P ± M P , (2-111) 

with 

/ip = £/i — q • Vf = 5[i — qvp cosO , (2.112) 

where the upper (resp. lower) sign refers to the d (resp. u) quasi particle. Using the analogous result for the hole 
with field ip c (—p), one can follow the same steps leading to (|2.44l) from (|2.38|) : therefore the gap equation is still given 
by 12.44f) . but now the quasiparticle occupation numbers are 

n ^ = e («+Pp)/r + 1 ' n ^ p ) = e (<-Pp)/r + 1 ' ( 2 - 113 ) 
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FIG. 4 Kinematics of the LOFF state in the case of one plane wave behavior of the condensate. The Cooper pair has a total 
momentum 2q 7^ 0. 



By <|2.113|) . using the gap equation for the BCS superconductor with gap Ao, the gap equation for the inhomogeneous 
superconductor is written as 

gp A Q g f d 3 p ,„,.,, 

T ln a" = 2 J Mp) Mp)] ■ ( } 

Differently from the case with equal chemical potentials (Sfi = 0), when there is phase space reduction at T ^ 0, now 
also at T = the blocking factors reduce the phase space available for pairing. As a matter of fact the gap equation 
at T = reads 

gp A g f d 3 P 1 _ 

Y ln X = 2 7 (2-7r) 3 e(p, A) W^-M + ^ + Mp)) 

gp /" dfi D , C(6) 

= £l / — E arcs i n h_^ (2.115 

where 



C(fl) = ^(z,-cos») 2 -A 2 (2.116) 

and 

z 9 = ^ = cos^, (2.117) 

where ipo is the angle depicted in Fig. 0] The angular integration is not over the whole Fermi surface, but only over 
region defined by e(p, A) < or 

q 2 v 2 F (z q - cos6>) 2 > A 2 . (2.118) 

Notice that there are no solutions to this inequality for qvp + 5/j, < A (compare with Eq. 12.52|l ). Analyzing this 
inequality in terms of cos# we see that there are three regions, obtained comparing qvf — dp, to ±A, characterize d by 
different domains of angular integration. They are displayed Table ||| As pointed out in l)Fulde and Ferreiill964|) . the 
blocking regions correspond to regions in momentum space where fermions do not pair. In regions E and S fcrmions 
of one type (for instance spin up) do not pair, whereas in region D fermions of both types do not pair. The effect of 
the blocking regions is to reduce the phase space where pairing is possible. The complementary phase space is where 
the pairing is possible and therefore it will be called pairing region. It is formed by two rings that loosely speaking 
are around the two circles of Fig. Since the pairing is possible not only on the Fermi surface, but also for modes 
just below and above it, each ring has a toroidal shape. ijj = 2arcos(z g ) is the aperture of the cone, with vertex at 
the origin of the spheres, intersecting the Fermi surfaces along the rings. 

Once fixed the integration domain, the remaining integral in cos 9 is trivial and the result can be expressed, for the 
three cases, in the following uniform way 



ln^= A 



A 2qv F 



q , qvF + Sfi \ +G fqv F - Sp 



(2.119) 
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Region 


Definition 


Domain of integration in cos 9 


E 


qVF — Sfi < — A 


(-!,+!) 


S 


-A < qv F -Sfi< +A 


(-1,COS0_) 


D 


qVF — Sfi > +A 


(-l,cos0_)|J(cos $+,+!) 



TABLE I In the table the three blocking regions are shown. Here we have defined cosO± = z q (1 ± A/ Sfi) 



where the function G(x) is defined as follows: 



G(x) — x arccosh(a;) — \/ x 2 — 1, |ac| > 1, 
G(x) = 0, \x\ < 1, 

G(x) = -G(-x), x<0. (2.120) 



1. Second Order phase transition point 

The reduction of the available phase space implies a reduction of the gap, therefore one expects in general smaller 
gaps in comparison with the homogeneous case. In particular we see from Eq. (|2.115|l that increasing Sfi the effect 
of the blocking terms increases; eventually a phase transition to the normal phase occurs when Sfi approaches a 
maximum value Sfi2- Therefore the anisotropic superconducting phase can exist in a window 

Sfix <Sfi< Sfi 2 . (2.121) 

One expects that 5fi\ is near the Chandrasekhar-Clogston ljChandrasekharlll962UClogstonLll962h limit Ao / V2 because 
Eq. (|2.63|) shows that near this point the difference in energy between the isotropic superconducting and the normal 
phases is small and one might expect that the LOFF state corresponds to the real ground state. We shall discuss 
the gap equation and prove this guess below. For the moment we determine Sfi2- For Sfi — > Sfi2 the gap A — * 0, and 
in the blocking regions E and D the domain of integration in cos# is (— 1, 1) (the region S disappears in the limit). 
Expanding the function G{x) for i-*cowe get from Eq. I|2.119|l 

m ^ - = ~ i + 77 l n ? m ^ , » — — (2.122) 

A 2 qv F qv F ~ Sfi 2 4(g 2 u| - Sfi 2 ) v ' 



which can be re-written as 



a {qVF , S,) = -1 + i^iln^±^ - \ in - 2 2 A ° = . (2.123) 
2 qvF qvF — Sfi 2 4{q z Vp — Sfi z ) 



In terms of the dimensionless variables 



V=^, (2.124) 



A ' A 



the condition a = is equivalent to the equation 



K + '-IlF?) ■ ( 2 - 125 ) 



2 \z-y 

The critical line is plotted in Fig. [S] 

Notice that the equation l|2.122|l can be written also in the form 

ln A = 1^ / gv F \ = 1 | 1 Sfi ^ qvp + Sfi 1^ Sfi 2 (2 126) 

ISfi 2 \ Sfi J 2 qvF qvF ~Sfi 2 {q 2 v 2 F — Sfi 2 ) ' 

where 

f (x)= du\n(l + xu). (2.127) 
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FIG. 5 The critical line for the LOFF phase at T = in the plane (qvF / Ao, Sfi/ Ao) ■ Also the line determining qvF as a 
function of 5^2 is given. 



We can fix q by minimizing the function a with respect to it. This is equivalent to minimize the grand potential close 
to the second order phase transition. This is obtained for x solution of the equation 

x = coihx, (2.128) 

i.e. at 

x = ^ = 1.1997 = x 2 . (2.129) 

0^2 

This can be also obtained from Fig. intersecting the curve a = at its maximum value 8^2/ Aq with a straight line 
passing from the origin. 

The value of at which the transition occurs is obtained by substituting this value in (|2. 126(1 and solving for 5/i2- 
One gets in this way 

6^2 = 0.754 A . (2.130) 

Since b\i2 > $/J>i ~ 0.71Ao, there exists a window of values of Sfi where LOFF pairing is possible. We will prove 
below, using the Landau- Ginzburg approach, that the phase transition for the one-plane wave condensate at T = 
and 5fi = 5fi2 is second-order. 



III. GINZBURG-LANDAU APPROXIMATION 

The condensate wave function acts as an order parameter characterized by its non vanishing value in the super- 
conducting phase. At the second order phase transit ion it vanishes and one can apply the general Ginzburg-Landau 
(GL) approach there l|Ginzburg and Landaul 119501) . We will begin by performing the GL expansion at T = 
for a general inhomogeneous gap function l)Bowers and Raiaeopaj. l2002t lLarkin and OvchinnikovL [1964]) . From this 
we will derive the grand potential measured with respect to the normal state and we will evaluate it explicitly 
for several cases. Next, in Section flll.DI we will pe rform an analogous expan s ion at T ^ around the tricritical 
point that we have shown to exist in Secti on III.BI ( Alexand er and McR.tag ueL Il978t iBuzdin and KachkachiL Il997t 
ICombescqt and Moral l2002l iHouzet and Buzdinl . l2000bHHouzet et aZJ . 120021 11999^ . We will follow in this discussion 
the Ref. nfc^mb^s^cTan^Moran2002|) . These authors have made a rather general analysis with the conclusion that 
in the generic case the favored state corresponds to a pair of antipodal wave vectors. 



A. Gap equation in the Ginzburg-Landau approach 



We will start this Section by considering the Ginzburg-Landau expansion of the Nambu-Gor'gov equations. Let us 
perform an expansion in A of the propagator F in Eq. (|2.28|) . It is depicted in Fig. 
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FIG. 6 Ginzburg- Landau expansion of the propagator; the lines represent alternatively G and Gq, see Eq. fS.lp . Full (resp. 
empty) circles represent A* (resp. A.). 



Formally it is written as follows 



Go A*Gq 
G A*G+AG A*G+ 

Gq A*Gq AGq A*Gq AGq A*Gq 



(3.1) 



The gap equation has an analogous expansion, schematically depicted in Fig. 




+ ^ + ^ — ^ + 



FIG. 7 Ginzburg-Landau expansion of the gap equation; the lines represent alternatively G and Gq , see Eq. iU.tJtl . Full (resp. 
empty) circles represent A* (resp. Aj. 



It has the form 



A* = 



3 

/ n^ j G -(r,r 1 )A*(r 1 )G+(r 1 ,r 2 )A(r 2 )G -(r 2 ,r3)A*(r 3 )G+(r 3 ,r) 

J 3 = 1 
5 

J n dr j G o _ ( r ' r i) A *( r i) G J( r i' r 2)A(r2)G (r 2 ,r3)A*(r 3 ) G+(r 3 ,r 4 ) 



xA(r 4 )G -(r4,r5)A*(r 5 )G+(r 5 ,r) 



Substituting (|2.109|l we get 



A n = (^n(q k ,q n )A^(q k -q n ) 

k 

+ J (qk,q^q m ,qn)AfeA^A^(5(q k - q e + q m - q n ) 



+ J2 ^(qk,q£,qmqj,qi,q„)At.A^A; i A J A* 

x ^(qk - qi! + q m - qj + qi - q n ) 



Here <5(q k — q n ) means the Kronecker delta: 6 n ,k and 



(3.2) 



(3.3) 



(3.4) 



J(qi,q 2 ,q 3) q4)= / ^ d£ | + J g f[ /,(£, f/*, {q}) 



(3.5) 
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#(qi,q 2 ,q 3 ,q4,q5,q6) = + '4r f [ g * / ^ Jl^^'^'^- 



We have put w = »pw and 

fi{E,5n,{q}) = 

Moreover the condition 



E + iesignE -6^+ (-!)<[£ -SELiH)^ ■ q k ] 



^(-l) fc qk = 



(3.6) 



(3.7) 



(3.8) 



fe=i 



holds, with M = 2, 4, 6 respectively for II, J and K. 
For 11(g) = II(q, q) one gets 



igp f dw 

47T 



+<5 



-6 



dE 



2tt (E + iesignE ~ jT) 2 - £ 2 



(3.9) 



where p, = dfi — vpq ■ w is defined in Eq. I|2.112|l and is identical to the function C(0) of Eq. (|2.116|) with A = 0. In 
performing the energy integration in (|3.9|) we use the fact that there are contributions only for |£| > \p,\. Using the 
gap equation for the homogeneous pairing to get rid of the cutoff 5 we obtain the result 



n( ff ) = i + f 



1 + x log 



A 2 



5fi 



log 



qvF + S/j, 



2 4| (qvF) 2 — Sfi 2 \ 2qvF qvF — Sfi 
11(g) can be rewritten in terms of the function a introduced in (|2.123|) as follows: 

1 - U(q) 



a(q) = 2- 



fJP 



(3.10) 



(3.H) 



Clearly the gap equation in the GL limit, 1 = H(q), coincides with Eq. I|2.123|) . which was obtained in the one plane 
wave hypothesis. The reason is that, since II depends only on |q|, it assumes the same value for all the crystalline 
configurations; therefore II does not depend on the crystalline structure of the condensate and the transition point 
we have determined in Sec. III. C. H is universal. 

For the evaluation of J and K we have to specialize to the different LOFF condensate choices. This will be discussed 
below. 



B. Grand potential 

The grand potential £1 is given in the GL approximation by 



= - 



,6, 



qk-q n 



i( WqicQn) - 1]A£.A„ 

& k.n=l 
I P 

i q^ i qm j qn )A* k A i A ni A n 6 l 

fc,^,m,n— 1 

i p 

Qm j Qj ) Qi i Qn ) A£ A t A* m Aj A* A n <5 qk _ q£ +qm _ ^ +qi _ q „ 



J qk-q«+q m -qn 



k,£,m,j,i,n— 1 

where P is the number of independent plane waves in the condensate. Let us assume that 



(for any k) 



(3.12) 



(3.13) 



so that we can rewrite (|3.12f) as follows: 
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where a is related to 11(g) through (|3.11|l and 

2 P 

P = J ( f lk,q^qm,qn)<5q k -q (! +q m -q n > (3-15) 

^ k,t,m,n=l 

2 P 

7 = J! if (qk,qf,qm,qj,qi,qn)6q k -q 4 +q m -q j +q i -q„ • (3.16) 

It follows from the discussion in subsection III . C . Tl that . at S[i = &[i2-, a vanishes; moreover a < for 5/i < S/i2, see 
below, eq. i|3.23|l . Exactly as in Section lil.B.ll we distinguish different cases: 

1. /3 > 0, 7 > 0. In this case A = is a maximum for f2, and the minima occur at the points given in Eq. 1)2. 96[) . 
which now reads: 



Near the transition point one has 



A 2 = ZI±VEIF^ . (3 . 17 ) 



A 2 «-^. (3.18) 



A phase transition occurs when a — 0, i.e. at Sfi = Sfj,2- The transition is second order since the gap goes 
continuously to zero at the transition point. 

2. f3 < 0, 7 > 0. Both for a < and for a > A 2 in (j3.17|l is a minimum for fl. In the former case it is the only 
minimum, as A = is a maximum; in the latter case it competes with the solution A = 0. Therefore the LOFF 
phase can persist beyond 8^2 , the limit for the single plane wave LOFF condensate up to a maximal value 5/j>* . 
At 5[i = <5/i* the free energy vanishes and there are degenerate minima at 

A = 0, A 2 = ^. (3.19) 
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The critical point S/i* is obtained by Eq. (|2.98|l that in the present case can be written as 

33 2 

a(qv F = 1.1997V, S/i*) = . (3.20) 

The phase transition from the crystalline to the normal phase at Sfi* is first order. 

3. j3 < 0, 7 < 0: In this case the GL expansion Ij3.14|l is inadequate since is not bounded from below and another 
term C(A 8 ) is needed. 

In the case f3 < 0, 7 > we can select the most favored structure by computing the free energy at a fixed value of 
Sfi. We choose Sfx = 8^2 where the FF state has a second order phase transition and a — 0. One has there 



7 p 12j z 



C. Crystalline structures 

For any crystalline structure the function a in the first term of the GL expansion is given by 

qvF + 5jji 



1 A 2 , 

a = 2 l0g 41(^)2- VI ' 2qv F 



8u , 
+ ^-log 



A 1 ( qv F 
= - l0g 26il + 2 f °{-Sj 



qv F — Sfx 



(3.22) 



where we have used Eqs. (I2.12(j|l . 13.10|1 and (|3.11l) : a vanishes for Sfi — 5/i 2 , which characterizes the second order 
transition point at T = 0, see (I2.123f) or l|3.18|l : therefore we can write 



JL 

SfJ-2 



(3.23) 
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where 

T) = <5/i 2 - 5(i (3.24) 

and we have expanded a around S(i2 and used the property of minimum of fo(x) at Sfi = 5(i2- We observe that, for 
Sfi < Sfi2, ol is negative; therefore the transition at T = is always second order if ft > 0. 

As to the other terms, we can use the results of Appendix El to get the first terms of the GL expansion for any 
crystal structure. The exception is the one-plane-wave case, where the free energy can be computed at any desired 
order. 

1. One plane wave 

Using the results of the Appendix El on e gets for the Fulde-Ferrel one plane wave condensate: 
t t 9P 1 „ _ u _ 9P (qvf) 2 + 3<5/i 2 

J - Ja ~ T {qv F y - 6j? ' K - Ka ~ 64 {{qv F y - «5 M 2 ] 3 ' { ' 

From (|3~25|) we get (x 2 = qv F /S^ 2 = 1.1997): 

1 _ 0.569 

3 + x| _ L637 
7 8^1(x| - 1)3 " + 5ni ■ (6 - 2b > 

Since @ > the 7-term is ineffective near the transition point and Eq. I|3.18|l gives 

A 2 = 4?7 (x\ - l) 5fi 2 ~ l.757r)8(j,2 . (3.27) 

We can get fi from Eq. (|3~T4) with P = 1 and from Eqs. $T%ty and The result is 

2 

H = -777 = -0-439 p (5(i ~ 5fi 2 ) 2 , (3.28) 
The same result could also be obtained using Eqs. (|2.59() and (|2.62l) . 

2. Generic crystals 

In the general case P ^ 1 and the evaluation of J and T^T is more complicated. First, one introduces Feynman 
parameterizations, then the integrals over energy, longitudinal momenta and angles are performed, along the lines 
sketched in Appendix lAl mainly based on ijBowers and Raiaeopall 12002^ . Next, one has to perform the integration 
over the Feynman parameters. To do this it is useful to draw two pictures: a rhombus, with lines formed by the four 
vectors appearing in J(qk, q^, q m , qn), implementing the condition qt — q^ + q m qn = 0, and an hexagon, with 
lines formed by the six vectors appearing in K(q^, q^, q m , qj, qi, q n ) that satisfy qk — q^ + q m — qj + qi — q n = 0, 
see Fig. |SJ Note that the rhombus and the hexagon need not be in a plane. The simplest example is provided by two 
plane waves. 

3. Two plane waves 

In this case P = 2; let the two vectors be q a , qt,, forming an angle ip; a simpler case is provided by an antipodal 
pair, q a = — qb = q and 1/3 = tt, with 

A(r) = 2Acosq- r . (3.29) 
To get j3 from (|A5|) one may notice that the integral J assumes two different values 

Jo = ^(qa,q a ,qa,q a ) , = J(q a ,q a ,qb,qb) (3.30) 

corresponding to Fig. [5] 
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FIG. 8 Rhombic and hexagonal configurations for the vectors qi. The vectors are assumed of the same length q and such that 
qi — q2 + q3 — q4 = for the rhombus and qi — q2 + q3 — q4 + qs — qe = for the hexagon. The vectors need not be all in 
the same plane. 



a a 




FIG. 9 The two rhombic structures corresponding to the integrals Jo and of Eqs. 13.3U6 . The indices a and b refer to the 
vectors q a and q;, respectively. 



Jq has been already computed, see Eq. (|3.25[) : on the other hand 



2?2\/cOS lb— 1 

arctan . v 
26 fi* 2x 2v /(cosV>- l)(2-a;|(l + cos V)) 



which for %j) = n gives 



Using rotation and parity symmetry of the integrals one gets 

0(ip) =-— (2J +4J^) . (3.33) 
9P 

The result for /3(t/>) as a function of V> is reported in Fig. In the case of the antipodal pair (q, — q), when %j) — 180°, 
one gets 



o 1/1 \ 1 38 



For if we have three possibilities (see Fig. Illfl : 
K = K(q a , q a , q a , q a , q a , q a ) , Ki(ip) = K(q a , q a , q a , q a , q b , q b ) , K 2 (ip) = K(q a , q a , q b , qb, qb, qb) ■ (3.35) 
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FIG. 10 /3 = /3 5/j,2 as a function of the opening angle ip between the two plane wave vectors q a and qt,; ip — 00 = 67.07° is the 
angle defining the LOFF ring; (3(ipo) = —1.138. 




FIG. 11 The three hexagonal structures corresponding to the integrals Kq, K\ and K2 of Eqs. IS. 356 . The indices a and b refer 
to the vectors q a and q& respectively. 



Therefore we have 

7(V>) = -— (2K + 12^(0) + 6K 2 (ip)) . (3.36) 
39 

Kq has been already computed in (|3.25|) . whereas Ki and K 2 can be evaluated using the results given in Appendix 
A. j(ip) is plotted in Fig. In the case of the antipodal pair, when ip = 180°, the result for 7 is in Table 2. 



Figs. ITUI and H2*l show a divergence at 

if} = tp = 67.07 = 2 arccos — . (3.37) 

qv F 

■00 is the opening angle depicted in Fig. In this case, differently from the one plane wave situation, we have two 
different rings for each Fermi surface. For ip > ipo the two rings do not intersect, at tp = ipo they are contiguous, while 
for ip < ip they overlap. The structure with ip < ip is energetically disfavored because, bei ng large and positive, 
the f ree energy would be smaller according to Eq. (|3.18(l . According to the discussion in l|Bowers and Rai agonal 
2002), this behavior seems universal, i.e. structures with overlapping rings are energetically disfavored in comparison 
with structures without overlaps. We will use this result in Section IV. Dl 

For ipo < ip < 132° (3 is negative. Therefore according to the discussion above we are in presence of a second order 
phase transition (7 is always positive as it can be seen from Fig. 1121) . As it is clear from Eq. H3.21fl the most favorable 
case from the energetic point of view occurs when 7 assumes its smallest value and \(3\ its largest, i.e. at ip = ipo, 
when the rings are tangent. The values for this case are reported in Table ITT1 

For comparison, at ip = 90° we have 5^2/3(90°) = —0.491, <5/i2 7(90°) = 1.032; the first order transition takes 
place at 5^* — 0.771Ao, only marginally larger than 6^,2, and the dimensionless free energy f2 = J7/(pA§) assumes at 
6/j, = S/i2 the value ft = —0.005, which is larger than the value obtained for ip = ipo, see Table ITT1 
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FIG. 12 7 = 7<5/i2 as a function of the opening angle ip between the two plane wave vectors q a and qt,; ip = ipo — 67.07° is the 
angle defining the LOFF ring; j(ipo) = 0.249. 



Structure 


P 


P 


7 




5/x*/A 


FF state 


1 


0.569 


1.637 





0.754 


antipodal plane waves 


2 


0.138 


1.952 





0.754 


Two plane waves (ip — ipo) 


2 


-1.138 


0.249 


-1.126 


1.229 


Face centered cube 


8 


-110.757 


-459.242 







TABLE II Candidate crystal structures with P plane waves. /3 = j3 t 7 = 5^1 > ^ — ^V(pAo), with p = p%/(tt 2 vf), is the 
(dimensionless) minimum free energy computed at 8/j, — S/j, 2 , obtained from W.21\) . The phase transition (first order for (3 < 
and 7 > 0, second order for f3 > and j > 0) occurs at Sfi* , given, for first order transitions, by Eq. \3. 2(A) . 




4. Other structures 



One could continue in the same way by considering other structures. An extensive analysis can be found in 
(Bow ers and Raiagopalll2002|) where 23 different crystalline structures were considered. We refer the interested reader 
to Table I in this paper, as well as to its Appendix where the technical aspects of the integration over the Feynman 
parameters of the K integrals for the more complicated structures are worked out. From our previous discussion we 
know that the most energetically favored crystals are those which present a first order phase transition between the 
LOFF and the normal phase. Among the regular structures, with 7 > 0, examined in ijBowers and Raiaeopall |.2002) 
the favored one seems to be the octahedron (P = 6), with 6 /a* = 3.625A . Special attention, however, should be given 
to th e face centered cube; we have reported the values of its parameters, as computed in l|Bowers and Raiaeopall 
2002), in our Table [n] We note that 7 < for this structure. The condensate in this case is given by 

8 8 

A(r) = A * ( r ) = I] A exp(2z g n k • r) , (3.38) 

k=l k=l 

where rik are the eight unit vectors defining the vertices of the cube: 

&i = -j=(+i,+i,+i), fi 2 = -j=(+i,-i,+i), 

fia = -^(-l,-l,+l), n 4 = -j=(-l,+l,+l), 
fi 5 = -^(+l+,l,-l), &e = (+1,-1,-1), 

n 7 = -L(-l,-l,-l), fi„ = -L(-l, +1,-1). (3.39) 

Strictly speaking, since both (3 and 7 are negative, nothing could be said about the cube and one should compute 
the eighth order in the GL expansion, given by i5A 8 /8; the transition would be first order if 6 > 0. However 
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iBowers and Rai agonal I2002D argue that, given the large value of 7, this structure would necessarily dominate. 
Reasonable numerical examples discussed by the authors confirm this guess. 



D. LOFF around the tricritical point 



The LOFF phase can be studied analytically around the tricritical point l)Buzdin and KachkacM Il997t 
ICombescot and Moral l2002fl that we have considered in Section III.BI Here we will follow the treatment of Ref. 
!|CombesccTandMoralT2002^ . The tricritical point is the place where one expects the LOFF transition line to start. 
Close to it one expects that also the total pair momentum vanishes, therefore one can perform a simultaneous expan- 
sion in the gap parameter and in the total momentum. Starting from the expressions given in Section IlIII ( see Eqs. 
(|3.4|) . I|3.5|) and l|3.6l) ) and proceeding as in Section III. Bl we find 



n = 



q 



fi(q) |A q | 2 + i^/3(q,)A qi A; 2 A q3 A; 4 + \ £ 7 ( qi ) A qi A q2 A q3 A; 4 A q5 A^ . 



(3.40) 



Here we have used the momentum conservation in the fourth order and in the sixth order terms 



qi + Q3 = q2 + q4, qi + Q3 + qs = q2 + Q4 + qe , 



(3.41) 



with 



a{q) 

/%*) 



7 



where a, (3 and 7 were defined in Eqs. 



7(0? 
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Qi 



and 



Q = qv F 



Ql + Qi • Q 3 + Q2 • Q 4 ) 



and 



(3.42) 



(3.43) 



In Appendix [B] we show, as an example, how a can be obtained from the expansion of H(q) around Q = 0. In order 
to get a coherent expansion one has to consider the modulus of the pair total momentum of the same order of the 
gap. In fact, as we shall sec, the optimal choice for Q turns out to be of order A. Correspondingly one has to expand 
the coefficient of the quadratic term in the gap up to the fourth order in the momentum and the fourth order term in 
the gap up to the second order in the momentum. In the form given in Eq. (|3.4U|) one can easily apply the general 
analysis around the tricritical point used in Section III.BI In particular for vanishing total momenta of the pairs we 
are back to the case of the homogeneous superconductor studied in Section lll.BI 

It is interesting to write the expression for the grand potential in configuration space, because it shows that around 
the critical point the minimization problem boils down to solve a differential equation, whereas in a generic point the 
Ginzburg-Landau equations are integral ones. By Fourier transformation we get from Eq. (|3 .401) 



d 3 r 



a|A(r)| 2 + -/3|VA(r)| 2 + - 7 |V 2 A(r)| 2 



/3|A(r)| 4 + - 7 (2(V|A(r)| 2 ) 2 + 3(VA 2 (r))(VA* 2 (r)) 



i 7 / d 3 r|A(r)| 6 . (3.44) 



Let us now recall from Section III. Bl see Eq. (|2.98|l . that the first order phase transition is given by: 



3first 




(3.45) 



with a discontinuity in the gap given by 



-4^ = 





3/3 
I7 



(3.46) 



see Eq. (|2.99|l . Let us now consider the possibility of a second order transition in the general LOFF case. Only the 
quadratic term in the gap is necessary for the discussion, and we have to look at its zero, given by a = 0. Since we 
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are considering only the quadratic term we can choose an optimal value for Q 2 by minimizing this term with respect 
to |Q|. We find 



Q 2 = -\-, (3-47) 

o 7 



requiring /? < 0. The corresponding value for a turns out to be 



^ R 2 

"-24T' (3 - 48) 



or 



/24 

/Second = -y ~^ a l ■ (3.49) 

The LOFF second order transition line is higher than the first order transition line of the homogeneous case, since 
/^second > /ferst (see Fig. showing the relevant lines in the plane (5/j,/A , T/A j). Therefore the second order 
transition to the LOFF case overcomes the first order transition to the homogeneous symmetric phase as it can be 
checked by evaluating the grand potential for the LOFF state along the first order transition line. 

The situation considered before corresponds to the physics of the problem only when the second order transition 
is a true minimum of the grand potential. This is not necessarily the case and we will explore in the following this 
possibility. 



tricritical point 




0.61 0.615 0.62 0.625 0.63 A 



FIG. 13 The graph shows, in the plane (Sfi/Ao,T/Ao), the first order transition line (dashed line) from the homogeneous broken 
phase to the symmetric phase. The solid line corresponds to the second order transition from the LOFF phase to the symmetric 
one. The lines start from the tricritical point and ends when the Landau- Ginzburg expansion is not valid any more. 



1. The LO subspace 



The second order term in the grand potential requires that the vectors Q have the same length along the second 
order transition line. It is natural to consider the subspace (LO) spanned by plane waves corresponding to momenta 
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with the same length Qo : 



A(r) = 22 Aqe 2iq ' r . (3.50) 
IQI=Oo 

The authors of Ref. ijLarkin and OvchinnikovL Il964f) have restricted their consideration s to periodic solutions, but 
this is not strictly necessary, although the solution found in l|Combescot and Moral 12002'! is indeed periodic. We will 
see that within this subspace the usual LOFF transition (the one corresponding to a single plane wave) is not a stable 
one. We will show that there is a first order transition that overcomes the second order line. In order that the LOFF 
line, characterized by Eqs. (|3.47[1 and l|3.48|l is a true second order transition the coefficient of the fourth order term 
should be po sitive. However, in the actual case, the mixed terms in the scalar products of the vectors Q could change 
this sign. In ijCombescot and Moral l2002|) the mixed terms are studied by defining the following quantity 

2 b Q 2 ]T A qi A q2 A q3 A q4 = J2 A qi A q2 A q3 A q4 (Q x Q 3 + Q 2 Q 4 ) . (3.51) 
q; q; 

Clearly 

— 1 < & < 1 , (3.52) 

where b = 1 is reached in the case of a single plane wave. With this definition and for the optimal choice of Qo (see 
Eq. H3.47JI ). we get for the coefficients appearing in the expression of the grand potential (see Eqs. H3.40J) and Ij3.42|l ^: 

a = a-|A = -1/3(56 + 1), 7 = 7 . (3.53) 



Therefore, for any order parameter such that 



it follows 



&<-~, (3-54) 
5 



0<O (3.55) 

and the LOFF line becomes unstable (we recall that close to the second order line (3 < 0). In fact, since 5 = along 
this line and (3 < 0, a small order parameter is sufficient to make ft negative. In other words one gains by increasing 
the order parameter as long as the sixth order term does not grow too much. But then we can make Q = (equal 
to its value in the symmetric phase) by increasing a. Therefore we have a new t ransition line in the plane ( a, (3) 
(or, which is the same in the plane (S/i/Ao, T/Aq)) to the right of the LOFF line. l)Combescot and Moral f2002) also 
shows that necessarily 

&>--' (3 ' 56) 

The equality is reached for any real order parameter A(r). In order to get a better feeling about the parameter b it 
is convenient to consider the following quantity 

2cQ 2 ^ A qi A* 2 A q3 A q4 =£(Qi - Q 2 ) 2 A qi A q2 A q3 A q4 . (3.57) 
qi m 

Expanding the right hand side of this equation and using Qi • Q 2 = Q3 • Q4 we find 

c=l-b. (3.58) 

To minimize b is equivalent to maximize c. To this aim it is convenient to have opposite Qi and Q2 because then 
(Qi — Q2) 2 reaches its maximum value equal to 4Q 2 . In this case we have 

A q = A_ q . (3.59) 



This is equivalent to require that the order parameter is real in configuration space. Of course, it is not necessary 
that the amplitudes of the different pairs of plane waves are equal. 
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To proceed further one can introduce a measure of the size of the order parameter 

J d 3 r|A(r)| 2 = ]T A qi A; 2 ee N 2 A 2 (3.60) 



and 



-L f d 3 r|A(r)| 4 = Y, A qi A; 2 A q3 A* 4 ee 7V 4 A 4 , (3.61) 

i / d 3 r|A(r)| 6 =J2 A qi A; 2 A q3 A; 4 A q5 A; 6 ee iV 6 A 6 . (3.62) 
^ J q. 

In the case of a single plane wave we get 

N 2 = N 4 = N 6 = 1 . (3.63) 

For a real gap the set of the vector s Qi is made of N/2 pairs. If all the plane wave have the same amplitude one can 
show ijCombescot and Moral l2002[) that the quantities N 2 , N4 and N@ assume the following values 

N 2 = 2N, Ni = 3N(N- 1), iV 6 = 5N(3N 2 - 9N + 8) . (3.64) 

With these notations the grand potential becomes 

n = N 2 (a + ^(3Q 2 + ^iQo) A 2 + ~iV 4 U + 7 Q^j A 4 + ^N 6l A 6 , (3.65) 



where 



a=^, l< a <i- ( 3 - 66 ) 



Minimizing this expression with respect to Qo we find 



<2o = - o- ~ 7 a 7T A . 3.67) 

87 4 7V2 



Therefore a non zero solution for Qq is obtained if 



1 0N 2 _ - A2 



A 2 <--^-^ = A 2 iax . (3.68) 
2a 7 JV4 

The corresponding expression for Q becomes 

n - ° N > (• - IS) A! + (' - r) A ' + r w « ^ ™) A *- <" 9 ' 

In order to have a transition from the symmetric phase we must allow fi/A 2 to become negative. The zero of ft is 
reached for 

5a^ 2 



5 /3 2 , 3 /3 2 , ( .,_ m 



The value of A corresponding to the zero of Q is given by 

5a 



2 _ 3(3N 2 V 1 3, 
A ~ 27^4 ^^V6_5 fl2 ^ ' ^- 7ij 
A/f 2° 
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showing that this is a first order transition. We have shown that it is convenient to have b as small as possible and 
that the minimum is reached for b = —1/3 or a = 5/9. In correspondence with this value we have A 2 > and 



We get 



a ~ YA~ + 972^ N 2 N 6 125- (373) 
Nj 162 



We see that it is convenient to take N 2 Nq/N% as small as possible. However, notice that 

N 2 N 6 



> 1 (3.74) 



as it follows from the Schwartz inequality 



[/AA 2 (r)] [KrA6(r)] ^ 
[Jd3rA4(r)] 2 



(3.75) 



Let us look for example at the case of pairs of plane waves with opposite q without any further constraint. In this 
case we get 

N 2 N 6 _ 5 3iV 2 - 97V + 8 
Nl "9 (iV-1) 2 ■ ( ' 

This expression has a minimum at N — 2 (here N is even), where it holds 10/9 and then it increases monotonically 
up to the value 15/9. Notice that for crystalline structures the situation could be different. For instance, in the 
case of the cube the values o f N4 is not the one given by Eq. (|3.64() . i. e. 168, but rather 216. However, in Ref. 
l|Combescot and Moral 12002^1 it is shown that this expression gets indeed its minimum value for N = 2. This result 
is obtained assuming that the plane waves form a generic set of antipodal vectors, which means that the only way to 
satisfy the momentum conservation is through the cancellation of each momentum with the opposite one in the same 
pair. This exclud es special configura tions where other arrangements of vectors could give a zero result. The authors 
of Ref. l|Combescot and Moral [2002) argue that this result should hold in general, but a complete proof is lacking. 
The value obtained for a for the two plane wave case is 

»^(^K(I+— ■)• 

In this case we have three lines, the first order line just found, the second order LOFF transition for a = 5(3 2 /(24 r y) 
and the first order transition to the homogeneous broken phase for a = 3/3 2 / (I67). The distance of these two last 
lines is given by 

P 2 f 3 5 \ 1 P 2 01., ,n~2P 2 



— = — = -2.1 x 10 -2 — . (3.78) 

7 V 16 24/ 48 7 7 

We know that these two lines stay close one to the other up t o zero temperature. It turns out that the same is 
true for the new first order line as it has been shown in Ref. ijMatsuo et all 11998^1 . These results are illustrated 
in Fig. El We notice that whereas in the expansion around the tricritical point the favored state seems to be the 
one corresponding to a pair of plane waves, with a first order transition between the LOFF and the normal state, 
at zero temperature one has a second order phase tr ansition. Therefore the first order transition line must change 
into a second order line at low temperatures. In ref. l|Matsuo et all Il998l) it has been shown that this happens at a 
temperature T/T BCS = 0.075. 

It is also interesting to se e how things change varying the spatial dimensions. In fact it has been found in 
l|Burkhardt and Rainerl Il994|) that the first order transition found previously is second order in two spatial dimen- 
sions . This result is confirmed bv llBuzdin and Kachkachl Il997t iBuzdin and Polonskil Il987l iBuzdin and Tueushevl 
Il983t iMachida and Nakanishil Il989fl which show that the transition is second order in one spatial dimension and 
furthermore it can be given an exact solution in terms of the Jacobi elliptic functions. This solution has the property 
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that along the second order transition line it reduces to the two plane wave case considered here. It is very simple 
to obtain the dependence on the number of dimensions. In fact the only place where the dimensions enter is in the 
angular integration as, for instance, in Eq. (|B2|) . In general this is an average over the D-dimensional sphere and we 
need the following equations for the terms of order Q 2 and Q 4 respectively: 

rfw „ „ Sij f dw „ „ „ A 1 

-g—WiWj = -^j-, J —WiWjWkWi = D - 2 ^ {dijdki + dikdji + dudjk) , (3.79) 
where Sx> is the surface of the sphere with unitary radius in D-dimensions 

s -rW (3 ' 80) 

Therefore the terms proportional to Q 2 are multiplied by 3/D and the ones proportional to Q 4 by 15/D(D + 2). We 
get 

«(q) ^ + >Q 2 + ^^7Q 4 , 

0(qO - /5 + ^7(Qi + Q2 + Q3 + 04 + Qi -Q3 + Q 2 -Q 4 ), 

7(<k) = 7- (3-81) 
Proceeding as before the grand potential evaluated at the optimal value of Q 2 is 

« - ■* (' - w£) * + >• (> - ^«) * + ^ (> - i6 ■ 

We see from here that the critical dimension for the transition to change from first to second order is at the zero of 
the second term, that is 

2a 

D = 3.83) 

1 — a 

and for a — 5/9, D — 2.5. This means that we have a first order transition for D > 2.5 and a second order one for 
D < 2.5 (remember that (3 < 0). The location of the transition is at (a = 5/9) 

a= £±lt + ^l (2^-5)' 

8D 7 20 7 D{7D - 10) y ' 

The value of the gap along the transition line is given by 

A' = f|^. (3-85) 
5 7 ID — 10 

We see that A 2 > for D > 2 5 . 

ijCombescot and Moral l2002h have considered also the possibility of solutions around the tricritical point not be- 
longing to the LOFF subspace. In fact, the antipodal solution does not satisfy the Euler-Lagrange equation obtained 
from Eq. (|3 .441) . due to the non-linear terms. If these are small it is reasonable to look for solutions close to A cos(q-r). 
This assumption simplifies the problem because the antipodal soluti on is essentially a one-dimensional solution char- 
acterized by the direction of q. Then l|Combescot and Moral l2002|) have found that the corrections at A cos(q ■ r) 
expressed in terms of higher harmonics are indeed very small. Of course this is only a consistency argument, but it 
is an indication that the choice of the LOFF subspace is a good approximation to the full problem. 

To conclude this Section let us say that in our opinion the status of the LOFF phase is not yet settled. Up to now 
we have considered the Ginzburg-Landau expansion both at T = and at the tricritical point. The results in the 
three-dimensional case can be summarized as follows: 

• Zero temperature point: In ijLarkin and Ovchinnikovl Il964|) it is found that the favored phase has a gap 
with a phase modulation cos(q- r) corresponding to a structure with two antipodal vectors. This phas e and the 
normal one are separated by a second order-transition line. However in ijBowers and Raiaeopall 12002). where a 
rather complete study of the possible crystalline structures has been done, it is argued that the most favorable 
structure would be the face-centered cube. The transition between the corresponding phase and the normal one 
should be first-order. 
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• Tricritical point: In l)Buzdin and Kachka ch1lll997ft the non-uniform phase has been studied in different dimen- 
sions with the result that the space modulation related to a single wave vector (i.e. exp(2zq • r)) is always unfa- 
vorable. These authors find also hat the solution with two antipodal wave vectors is the preferred one. In 1 and 
2 space dimensions the transition to the nor mal state is second-o rder, wher eas it is first-order in 3 spac e dimen- 
sions. Analogous results have been found in l)Houzet et aZJ . 11999^1 (see also l|Agterberg and YangM2o6"H) '). where 
the study has been extended to space modulations such as cos(qx) + cos(qy) or cos(qx) + cos(qy) + cos(qz). These 
authors argue that there could be various tr ansition lines at temperatures lower than the tricritical point. Finally, 
as th oroughly discussed in this Section, in ijCombescot and Moral 12002ft the results of ijBuzdin and Kachkachil 
11997ft are confirmed. 

Th ere have been also numerical investigations about the full phase space. In particular (Bur khardt and Raineit 
1994) proved that in the two-dimensional case (layered superconductors) the phase transition from the normal phase 
to the one characterized by two antipodal vectors is second order. The second-order trans i tion line from the phase 
with a single plane wave to the normal phase has been studied in llSaint- James et a/.Ul969ilSarmaLll963ft . The two- 
dimensional case for type II superconductors has been studied in l)ShimaharanT998a|) . The author has considered states 
corresponding to single wave vectors and antipodal pairs together with configurations corresponding to triangular, 
square and hexagona l states. It is found that, according the temperature, all these states may play a role. Finally 
l|Matsuo et a^J . fl998^ . as already mentioned, make a numerical analysis based on the use of quasi-classical Green's 
functions. They find that in 3 dimensions the transition line between the normal phase and the antipodal vectors phase 
starts being first-order at the tricritical point and becomes second-order at T = 0.0075 Tecs- The two-dimensional 
case will be discussed again in Section IYl.DI 

In conclusion the question of the preferred non-uniform state cannot be considered settled down yet. Let us discuss 
by way of example the 3 dimensional case. We have seen that there are strong indications that the favored state 
around the tricritical point is the one corresponding to two antipodal vectors. This being the case the natural question 
is: how the transition line extends down to zero temperature? If at T = the preferred state were the antipodal 
pair a further tricritical point in the plane {5fx,T) would arise. In fact, rec all that the transition is first o rder at the 
tricritical point and second-order at T = 0. However, if the conjecture in ((Bowers and Ral agooal [2002) is correct, 
the cubic phase would emerge in t he path going to T = 0. A possibility is that one goes from one structure to another 
in analogy to what suggested by ijShimaharal Il9 98a1 for the two-dimensional case. The other logical possibility is 
that the Combescot and Mora result at the tricritical point might be evaded by an exceptional arrangement of the 
wave vectors as, for instance, in the case of the face-centered cube. Therefore we think that more theoretical work is 
necessary in order to fill-in these gaps in our understanding of the non-uniform superconducting phase. 



IV. SUPERCONDUCTIVITY IN QUANTUM CHROMODYNAMICS 



Co l or superconducti vity (CSC) is an old subject ((Bailin and Lovel ^984; B arroist Il977t ICollins and Perrvl 
119751 IFrautschl | 1978ft that has recently become one of the most fascinat i ng research fields in Quantum 
Chromo Dynamics (QCD); these developments can be found in ( Agasian et al, 1999; Alford et all 119994 2000 > 
19981 Il999bl ICarter and Diakonovi Il999t iHong et all l2000t iPisarski and Rischket Tl999al l2000bl iRapp et all ll99S : 
Schafer and Wilczek ^ Th999alb]| cflcl|: Shuster andSonl ~ l2000ft : for reviews see ijR^agop^lMio^WiTczekl 12001ft . ijllong , 



20011. iAlfordl200ll .lHsu i l2000ft . It offers a clue to the behavior of strong interactions at high baryonic densities, 



an issue of paramount relevance both for the understanding of heavy ion collisions and the physics of compact stars. 
Color superconductivity arises because for sufficiently high baryon chemical potential \i and small temperature, the 
color interaction favors the formation of a quark-quark condensate in the color antisymmetric channel 3. In the 
asymptotic regime it is also possible to understand the structure of the condensates. In fact, consider the matrix 
element 



(4.1) 

1, • • • , N are flavor indices. Its color, spin 



W£<|o> 

where a, @ = 1,2,3 are color indices, a,b = 1,2 are spin indices and i, j 
and flavor structure is completely fixed by the following considerations. 

• Antisymmetry in color indices (a, (3) in order to have attraction. 

• Antisymmetry in spin indices (a, b) in order to get a spin zero condensate. The isotropic structure of the 
condensate is favored since a larger portion of the phase space around the Fermi surface is available. 



• Given the structure in color and spin, Pauli principles requires antisymmetry in flavor indices. 
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Since the quark and spin momenta in the pair are opposite, it follows that the left(right)-handed quarks can pair only 
with lcft(right)-handed quarks. In the case of 3 flavors the favored condensate is 

3 

WS^IO) = - WS^IO) = A £ £ Q/3 %c . (4.2) 

c=i 

This gives rise to the so-called color-fiavor-locked (CFL) phase ijAlford et all . Il999r3 ISchafer and Wilczekl Il999a|l . 
However at moderate densities other less attractive channels could play a role lAh^rdeT'aZn^fhSl) . The reason for 
the name is that simultaneous transformations in color and in flavor leave the condensate invariant. In fact, the 
symmetry breaking pattern turns out to be 

SU(3) C (g) SU(3) L ® SU(3) R ® U(l) B SU(3) C+L+R O Z 2 , 

where SU(3) c +l+r, is the diagonal subgroup of the three SU (3) groups. Both the chiral group and the color symmetry 
are broken but a diagonal SU (3) subgroup remains unbroken. The Z 2 group arises from the invariance of the 
condensate when the quark fields are multiplied by -1. We have 17 broken generators; since there is a broken gauge 
group, 8 of these generators correspond to 8 longitudinal degrees of the gluons, because the gauge bosons acquire 
a mass; there are 9 Nambu Goldstone bosons (NGB) organized in an octet associated to the breaking of the flavor 
group and in a singlet associated to the breaking of the baryonic num ber. The effective theory describing the NGB 
for the CFL model has been studied in l)Casalbuoni and Gattol Il999|) . 

This is the typical situation when the chemical potential is much larger than the quark masses m u , m<i and m s (in 
these considerations one should discuss about density depending masses). However one can ask what happens when 
decreasing the chemical potential. At intermediate densities we have no more the support of asymptotic freedom, but 
all the model calculations show that one still has a sizeable color condensation. In particular if the chemical potential 
|U is much less than the strange quark mass one expects that the strange quark decouples, and the corresponding 
condensate should be 

WS^JJQ) = Ae^%- , (4.3) 

since due to the antisymmetry in color the condensate must necessarily choose a direction in color space. Notice that 
now the symmetry breaking pattern is completely different from the three-flavor case, in fact we have 

SU(3) C <g> SU{2) L <8> SU(2) R ® U(1) B SU(2) C ® SU(2) L ® SU(2) R $ 17(1) <8> Z 2 . 

The chiral group remains unbroken, while the original color symmetry group is broken to SU(2) C , with generators 
T A corresponding to the generators T 1 , T 2 , T 3 of SU(3) C . As a consequence, three gluons remain massless whereas 
the remaining five acquire a mass. Even though the original U(1)b is broken there is an unbroken global symmetry 
that plays the role of U(1)b- As for U(1)a, this axial symmetry is broken by anomalies, so that in principle there is 
no Goldstone boson associated to its breaking by the condensate; however at high densities explicit axial symmetry 
breaking is weak and therefore there is a light would be Goldstone boson associated to the breaking of the axial 
U(1)a- One can construct an effective theory to describe the emergence of the unbroken subgroup SU(2) C and the 
low energy excitations, much in the same way as one builds up chi ral effective lagrangian with effective fields at zero 
density. For the two flavor case this development can be found in IjCasalbuoni et all 2000; R ischke et all 120011. 

It is natural to ask what happens in the intermediate region of \x. It turns out that the interesting case is for 
fi rj Mf/A. To understand this point let us consider the case of two fermions, one massive, mi = M s and the other 
one massless, at the same chemical potential \x. The Fermi momenta are of course different 

p Fl = vV -M^, Pf 2 =V- (4-4) 
The grand potential for the two unpaired fermions is 

n unpai , = 2 f 1 ^ ( V^TMf - M ) + 2 If, (Ip I " m) ■ (4-5) 

For the two fermions to pair they have to reach some common momentum p^, mm , and the corresponding grand 
potential is 

<W. - 2 ^ ^ ( V^TMl - m) + 2 $L (bl - M ) - ^ , (4-6) 
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where the last term is the energy necessary for the condensation of a fermion pair, see Eq. i|2.71|) . The common 
momentum pf omm can be determined by minimizing f2 pa i r . with respect to Pcomnr The result is (expanding in M s ) 

M 2 

Pfomm = M - • (4-7) 

It is now easy to evaluate the difference ^unpair. — impair, at the order , with the result 

ftpair. " ^unpair. « ^ (M? - 4A 2 /?) . (4.8) 

We see that in order to have condensation the condition 

M 2 

must be realized. The problem of one massless and one massive flavor has been studied by ijKundu and Raiaeopal 
2002). However, one can simulate this situation by letting the two quarks being both massless but with two different 
chemical potentials, which is equivalent to have two different Fermi spheres. The big advantage here is that one can 
use the LOFF analysis discussed in Section Hi. CI 

Color superconductivity due to the non vanishing of the condensates lj4.Hl or Ij4.3|l results from a mechanism analo- 
gous to the formation of an electron Cooper pair in a BCS superconductor and, similarly to the BCS superconductivity 
the only relevant fermion degrees of freedom are those near the Fermi surface. Therefore a two-dimensional effec- 
tive field theory has been developed. We shall briefly review it below, but our main interest is to delineate another 
development of color superconductivity, i.e. the presence of a LOFF superconducting phase. Also in this case the 
condensation is generated by the attractive color interaction in the anti t riplet channel . This phase of QCD has been 
main l y studied at small temperatures, see e.g . l)Alford et all l2001albl : iBowers et all l200lt iBowers and Raiaeopall 
120021: iLeibovich eTal\ . l200ll iRaiagopall l200l|) . Similarly to the CFL and 2SC phases the QCD LOFF phase can 
be studied by the effective theory, as shown in l|Casalbuoni et all l2002bl l2001al [2002clel iNardullit l2002bl) . This de- 
scription is useful to derive the effective lagrangian for the Goldstone bosons associated to the breaking of the space 
symmetries, i.e. the phonons. It is based on an analogy with the Heavy Quark Effective Theory and is called High 
Density Effective Theory (HDET). To describe these developments we organize this Section as follows. In Subsection 
IIV.AI we give an outline of the HDET. We specialize the formalism to the CFL phase in Subsection IIV.BI and to the 
2SC phase in Subsection IIV.CI The final Subsections are devoted to the LOFF phase in QCD. In Subsection IIV.DI 
after a general introduction to the subject, we consider a Nambu-Jona Lasinio coupling for a QCD liquid formed by 
quarks with two flavors. Given the similarities with the BCS four-fermion interaction arising from the electron phonon 
interactions in metals, we can apply the same formalism discussed in previous Sections. In the present case, however 
the two species we consider are quarks of different flavors, up and down, with different chemical potentials [i u , \id- We 
limit or analysis to the FF one plane wave state. However the results of Subsection IIII. Cl are valid also for the QCD 
LOFF state; in particular the guess on the favored structure at T = discussed in l)Bowers and Raiagopalll2002|) and 
reviewed in Subsection IIII. C .41 should point to the cubic structure as the most favorable LOFF crystal. In Subsection 
HV.EI wc discuss the differences induced by considering the one-gluon interaction instead of the effective four fermion 
interaction. LOFF superconductivity in QCD can be induced not only by a difference in the quark chemical potential 
but also by mass differences among the quarks. This situation is discussed in Subsection llV.FI that shows the role the 
strange quark mass can play in favoring the LOFF phase. 



A. High Density EfFective Theory 



At very high baryonic chemical potential [i and very small temperature (T — > 0) it is us e ful to adopt an effective 
description of QCD k nown as High Den sity Effective Theory (HDET), se e l)Beane et all l2000t ICasalbuoni et all 
l2001btlHongll2000albD . and, for reviews, llCasalbuonl 1200 It iNardulTH l2002al) . Let us consider the fermion field 



M*) = I 7^4 e ~ lP ' X ^P) ■ ( 4 - 10 ) 



d A p 
(2^ 

Since the relevant degrees of freedom are those near the Fermi surface, we decompose the fermion momentum as 

p" = fiv" + , (4.11) 

where w M = (0,v), v the Fermi velocity (for massless fermions |v| = 1) and V 1 is a residual momentum. We also use 
V = (1, v) ,W = (1, -v). 
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We now introduce the velocity-dependent positive-energy tp v and negative-energy ^> v left-handed fields via the 
decomposition 

*K*) = J ^ e ~^ V ' X &v( x ) + *v(i)] • (4.12) 



Here 



and 



'vi-D =,"-- P + i •(./•)= / -gt-e-^P+m (4.13) 



* v (x) = e^-'P-^x) = / -_ e -**P_V>(*) . (4.14) 
P± are projectors that for massless quarks are defined by 

P± = P±(v) = i(l±a-v) . (4.15) 

The extension to massive quarks is discussed in l)Casalbuoni et qfl l2002al 12003). The cut-off 5 satisfies 5 -C /U while 
being much larger than the energy gap. 
Using the identities 

^7^ = ^7^* v , tfv7 M Vv = *v7^v , (4.16) 

and substituting into the Dirac part of the QCD lagrangian we obtain 

— [i>tiV ■ Dip v + Vl(2fi + iV • D)* v + (<0 v i#i* v + h.c.)] ; (4.17) 

p± = £> M 7^ and is the covariant derivative: = + igA 1 ^ 1 . We note that here quark fields are evaluated at the 
same Fermi velocity; off-diagonal terms are subleading due to the Riemann-Lebesgue lemma, as they are cancelled 
by the rapid oscillations of the exponential factor in the /i — > oo limit (Fermi velocity superselectio n rule). A similar 
behavior occurs in QCD in the mq — > oo limit, when one uses the Heavy Quark Effective Theory dEichten and Hill 
19901 lGeorgill990t llseur and Wisel Il989l ll99o|) , and for reviews, llCasalbuoni et all Il997t iManohar and Wisel 1200(1 



Ne 1 ihertlll994f 



We can get rid of the negative energy solutions by integrating out the \l/ v fields in the generating functional; at tree 
level this corresponds to solve the equations of motion, which gives 

iV-Dip v = (4.18) 

and 

*v = "- ^—70^^, (4.19) 

2/i + iv ■ D 

which shows the decoupling of ty v in the fi — > oo limit. In the resulting effective theory for ip v only the energy and 
the momentum parallel to the Fermi velocity are relevant and the effective theory is two-dimensional. 
It is useful to introduce two separate fields 

V>± = ^±v ; (4.20) 
therefore the average over the Fermi velocities is defined as follows: 

£ - / 1 ■ <«" 



The extra factor 1 /2 occurs here because, after the introduction of the field with opposite velocity ip_ , one doubles 
the degrees of freedom, which implies that the integration is only over half solid angle. 
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In conclusion, if Co is the free quark lagrangian and C\ represents the coupling of quarks to one gluon, the high 
density effective lagrangian can be written as 



where 



C D = C + C 1 +C 2 + (L — > R) . 

C Q = [iP+iV ■ dip+ + ipUv ■ dip- 

V 

Ci = ig^2 [i'+W ■ + tp-iV ■ M 



and 



£ 2 = -^P^ Ul 



1 



2[i + iV D 



D^D^ + ip 



1 



2(j, + iV ■ D 



(4.22) 

(4.23) 
(4.24) 

(4.25) 



C2 is a non local lagrangian arising when one integrates over the ^ v degrees of freedom in the functional integration. 
It contains couplings of two quarks to any number of gluons and gives contribution to the gluon Meissner mass. We 
have put 



l 



(4.26) 



This construction is valid for any theory describing massless fermions at high density provided one excludes degrees 
of freedom far from the Fermi surface. 



B. CFL phase 

Even though we shall consider the LOFF phase only for two flavors, for completeness we present HDET for the 
3-flavor Color Flavor Locking phase as well. In the CFL phase the symmetry breaking is induced by the condensates 

(^JC^p) = -« T Cty&> = ~ ea/nem , (4.27) 

where ip L,R are Weyl fermions and C — iai. Eq. (|4.27|) corresponds to the invariant coupling (tp = ipL)' 

^C^i ~{L^R) +h.c, (4.28) 

1=1,3 

and (ez)„(, = £/a&- Neglecting the negative energy components, for the Dirac fermions ip± we introduce the compact 
notation 

J_ / i> 

in a way analogous to Equation (|2.11() . We also use a different basis for quark fields: 



* = ^{ok> (4 ' 29) 



(Aa) 

\ 

The CFL fermionic lagrangian has therefore the form: 

9 



fav*>=£;^fe*±- («o) 



v A,B=1 

where 



A AB = ATridT^Te] (4.32) 
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and 

TA = ^ ' (4 ' 33) 
Here A9 = Ao = x 1. We use the identity (<? any 3x3 matrix): 

ej<7 T e/ = 5 - Trg ; (4.34) 

we obtain 

A^ = A A S AB , (4.35) 

where 

A x = • • • = A 8 = A (4.36) 

and 

A 9 = -2A . (4.37) 
The CFL free fermionic lagrangian assumes therefore the form: 

£0 + £A = E E X At f Z T A f T^) X A + (L^R). (4.38) 



v A=l 



Clearly the equations of motion following from this lagrangian are of the same type as the NG eq uations, see Eq 
(|2.15() . For applications of HDET to the CFL phase we refer the reader to ijCasalbuoni et 



3 JNtj_equ 
l l2001bh . 



C. 2SC phase 

For the two flavor case, which encompasses both the 2SC model and the existing calculation in the LOFF phase, 
we follow a similar approach. The symmetry breaking is induced by the condensates 

(tfJC^) = - {^JCiifp) = y e af33 e ij3 , (4.39) 

and the invariant coupling is (ip = 

£ A = ~ip T Ceipe - (L -> R) + h.c. , (4.40) 

where 

e = i<r 2 . (4-41) 
We use a different basis for the fermion fields by writing the positive energy effective fields ip± v .i a as follows: 

= E %=r ^ • (4.42) 

A=Q V * 

The X A matrices are defined in terms of the usual A matrices as follows: 

Ao = ^A 8 + y|A , ~X A = X A (A = 1,2,3), A 4 = X 5 = . (4.43) 

We also define e = 1X2- After the introduction, analogously to 14.29fl . of the fields x A , the 2SC fermionic lagrangian 
assumes the form: 

Cd = £q + £1 + Ca 
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E E a 

v A,B=0 



Ai(iTr[f A V-Df B ] -A AB \ B ^ 



Here 



A 
2~ 

A AB =0 (A,B = 4,5) . (4.45) 



A AB = -Tr[i\ T A i\ B ] (A,B = 0,...3), 



and 



Ta = ^ (A = 0,...,5). (4.46) 

Analogously to (|4.34|) we use the identity: 

lg T e = g-Trg; (4.47) 

we obtain 

A AB = A A 5 AB , (4.48) 

where A^ is defined as follows: 

A A = (+ A, -A, -A, -A, 0, 0) . (4.49) 
Therefore the effective lagrangian for free quarks in the 2SC model can be written as follows 

£o + Ca = E E X* ( *V 7v A B ) X A + (L^R). (4.50) 



v A =0 



D. LOFF phase in QCD 

We shall assume here that in the most interesting phenomenological applications, i.e. in compact stars (see Sec- 
tion [^^GJ, there is a significant difference between the Fer mi momenta of different flavors. Since this produ ces a 
difference in the densities, the BCS phase may be disrupted l|Alford et all Il999al ISchafer and Wilczekl . Il999cj) and 
a phase analogous to the LOFF phase might ar ise. The case of a LOFF phase in QCD was also discussed in Ref. 
l)Son and Stephanovl 120011 ISplittorff et a/.U200l|) in the context of quark matter at large isospin density. Differences 
in the Fermi momenta in these examples arise both from the difference in the chemical potential, due to the weak 
equilibrium, and from the mass difference between the strange an d the up and down quarks. A c omplete study requires 
to take into account both effects. This has been made in Ref. l)Kundu and RaiaeopalL l2002|) . We will discuss this 
paper below. Here w e will consider a simpler case where all quark are massless but have different chemical potentials 
( Alfo rd et all l2001bj) . To simplify further the problem we will restrict ourselves to the case of two massless quarks 
with chemical potentials y, u and /id given by 

jtt« = A + 5(1, (j,d = (J, - 5(1. (4-51) 

These equations are the same as (|2.2f) but now up and down refer to flavor. 

Everything goes according to the discussion made in Sections III.BI and III. 01 except that now the density of gapped 
states at the Fermi surface is multiplied by a factor 4, coming from the two colors and the two flavors. In fact, the 
condensate has the form 

(V?^) « e a ^e l3 , (4.52) 

where a, (3 — 1,2,3 and i,j = 1,2 are respectively color and flavor indices. Other differences are in the value of the 
Fermi velocity, which is vp = 1, since we deal with massless fermions, and in the Fermi momentum which is given by 
Pf = fj,. As a consequence the density of states is now 



(4.53) 
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It follows that the first order transition from the homogeneous phase to the normal one, in the weak coupling limit, 
is given, using l|2.63|) . by 

2 

^a(^) - Wfi) = ^(25fi 2 - A 2 ,) = - A 2 ,) = . (4.54) 

Applying the results obtained in the Section ril.Cl to color superconductivity requires some care. For instance, although 
only two colors are ga pped, in order t o describe the mixed phase it is necessary the use of a proper treatment of the 
two ungapped quarks ijBedaauei l2002|) . Another situation that can be present in QCD but not in condensed matter is 
the case of equal chemical potentials with differen t Fermi momenta due to unequal masses. This is discussed in Ref. 
(Alfo rd et all Il999at iSchafer and WilczekL Il999cft . However in the realistic case different chemical potentials must be 
considered. 

We will describe also the LOFF phase using the formalism of fields close to the Fermi surface, although in the 
present case the corrections to the leading order are expected to be larger since we are not considering asymptotic 
values of the chemical potential. This formalism is very close to the NG formalism developed in Section [H] We 
consider a four-fermion interaction modelled on one-gluon exchange, that is 

Ci = --Gty^A'V VV'A'V i ( 4 - 55 ) 
8 

where A a are Gell-Mann matrices. We then introduce the fields ip"^ through the procedure outlined in Section llV.CI 
We perform the same transformation exp(— ifiv ■ x) for both flavors. For simplicity, in the rest of this Section we will 
denote the fields ij)f + by i/>f . Separating the left-handed and the right-handed modes the previous interaction can be 
written as 

Ci = ~ - S$P S ) e hi e bd i/>% i> S jd = V$ ^ ^ tftf , (4.56) 

where in the last expression the sum over the spin indices, d, c, b, d is understood and 



In obtaining this result we have used the identities 

8 2 

5I(A a ) a/3 (A a ) 57 = -(MarySfu ~ 5 aP 5 lS ) (4.58) 

a=l 6 

and 

{<Jn)db{^)dc = Zeactbd ■ (4.59) 

Here 

^ = (l,<r), cr" = (l,-<r), (4.60) 
with cr the Pauli matrices. As in Section [H] we divide Cj in two pieces 

C cond = VfftyWjitfl'l) + tf^MW) + (L^R) (4.61) 

and 

C mt = V$(1>tM S - (^W>)(WVj 4 - (V-f + (L -> R) . (4.62) 
The first piece can be written as 



where we have defined 



^a/sse^fVfAe 2 **' + c.c.) + (L -> R) , (4.63) 



rse 24q - r = -i £ Q/3 %(C^) (4.64) 
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and 

A = GTg . (4.65) 
The quadratic part of the lagrangian C CO nd in terms of the NG fields can be written as 

4L = 5Exi ta ( s-1 )54> ( 4 - 66 ) 



2 

where, in momentum space 

(q -upi = (5 aP (6 tJ V ■ £ + Sp(a 3 ) l3 S(£ - £')) -e aP3 e tJ A6(£ - £' + 2q) \ 

1 )aj \ -e am e ij AS(£-£'-2q) S^faV ■ £ + 6^)^6(1 - £' )) J ' { ' ) 

and g M = (0,q). Using this expression and performing the same derivation as in Section^ we find the gap equation 

A = 2 l9 A [ ^fE^Tr i . (4.68) 

J An 2tt "2^ 2 (V r -(^ + q)) + ^3)(^-(^-'Z))+^3)-A 2 V ; 

Comparing this result with Eq. (|2.39|l . one notices a factor 2 coming from the trace on the color; the remaining trace 
is on the flavor indices where the matrix 173 acts. Performing the trace and making explicit the ie prescriptions for 
the energy integration we find 

gp A f dv f +S „dE 1 

where p = 4/i 2 /it 2 is the relevant density at the Fermi surface and we have defined E = £0 and £ = £» to emphasize 
the similarity with the original LOFF equation (see Section Pl.Cty . Moreover 

jtZ = fy-v-q. (4.70) 

Performing the integration over the energy we get 
Since 

6(e - \fi\) = 1 - 8(-e - fi) - 6(-e + p) , (4.72) 

we get exactly the LOFF gap equation (compare with Eq. I|2.115|l ). except for the different definition of the density 
of states. 

We have already shown that in the present case there is a first order transition in 5p, between the homogenous 
state (that from now on will be referred to as the BCS state) and the normal state. Furthermore, from Section Til. CI 
we know that there is a second order transition between the LOFF state and the normal one. These results are valid 
also in the present case with the only change in the density of gapped states at the Fermi surface, which, as already 
stressed, is now a factor of four larger than the one for electrons. We recall from that analysis that around the second 
order critical point Sp2 = 0.754A (with A the BCS gap) we have (see Eq. (|3.27|l ) 



A i( ,/ / = x 'l.7r)7.VM<V<2 -'V> : L1~.A„t/ ^ 2 (4.73) 



As for the grand potential, we have from Eq. 

^BCS ~ ^normal = ^p{25p 2 - A 2 ) (4.74) 

and from Eq. (|3~25|) 

^LOFF ~ ^normal = ~0.439 p(Sp - Sp 2 ) 2 ■ (4-75) 



These results are summarized in Fig. 1141 where we plot the grand potentials for the different phases. 
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FIG. 14 The figure shows the differences between the grand potential for the BCS and the normal state and that between the 
LOFF and the normal state plotted vs. Sfi. The grand potentials are normalized to pA§. The inset shows the intersection of 
the two curves close to Sfii. The solid lines correspond to the LOFF case, whereas the dotted ones to the BCS case. 
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FIG. 15 The figure shows the condensates of the BCS and LOFF phases vs. Sfi. 



Since the interval (<5/ii, 8^,2) is rather narrow there is practically no difference between the values of S/i corresponding 
to the BCS-normal transition (S/ii = Ao/\/2) and the value corresponding to the BCS-LOFF transition. This can 
be visualized easily from Figs. 1141 and 1151 The figures have been obtained by using the previous equations in the 
Ginzburg-Landau expansion around 5fi2, but th ey are a very good approximation of the curves obtained numerically 
(|Alford et all l2000t iTakada and izuvamaLll969j) 



All t he discussion here has been done in the weak coupling limit. For a more correct treatment see (|Alford et all 
2000) where the results from the numerical integration of the gap equation are given. In particular we want to stress 
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the results on the size of the window obtained by these authors. If Sfii and 5fi2 are evaluated for a general coupling 
and at the same time one takes into account corrections from the chemical potential in the measure of integration, 
the windows get smaller and smaller for increasing BCS gap Ao. The corrections in the chemical potential arise from 
the momentum integration, which is made on a shell of height 25 but with an integration measure given by p 2 dpdfl, 
rather than ppdpdfl as usually done in the treatment of the BCS gap in the weak coupling limit. The results are 
illustrated in Fig. ^3 where the behavior of the critical points vs. the BCS gap are shown. The curves are plotted for 
a range of values of the cutoff A ranging from 0.8 to 1.6 GeV. The cutoff dependence is not very strong, in particular 
for 5[i2, moreover the window closes for A between 80 to 100 MeV, according to the chosen value of A. 
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A 

FIG. 16 The figure shows the critical points S[ii (solid lines) and 8fi2 (dashed lines) in Ao u nits vs. Ap, for cuto ff values 
A = /i + S ranging from 0.8 GeV up to 1.6 GeV and for /i = 0.4 GeV. This figure is taken from \Alford et all \2001li) . 

In llAlford et all . [2000) the presence of a vector condensate in the QCD LOFF phase is also discussed. The reason 
why this condensate can be formed in QCD but not in condensed matter is the following. Both in the BCS and in 
the LOFF phases the coupling is between fermions of the same helicity. In the BCS phase the fermions have also 
opposite momentum giving rise to a J — pair. On the other hand in the LOFF phase momenta are not exactly 
aligned, therefore a small component of J = 1 condensate may arise. A spin 1 state is symmetric in the spin indices 
and therefore Fermi statistics forbids it for electron pairing. On the other hand in QCD with two flavors one can form 
a state antisymmetric in color and symmetri c in flavor, and the Pauli principle is satisfied. Therefore the structure of 
the vector condensate is llAlford et aZl EoOO) 

|q| 

The ratio Ty/Ts is practically constant within the LOFF window varying between 0. 121 at and 0.1 33 at 5/i2- 
However this condensate does not contribute to the grand potential in the present case ijAfford et all Eoool) : therefore 
it does not change the original LOFF results. The situation is different if, instead of using the NJL interaction 14.55fl . 
use is made of the following interaction 

Cj = -f [G B (^7°A»(^7°A» - Gm(^7 1 A»(^7 1 A»] . (4.77) 

o 

This expression is not Lorentz invariant, but since we are trying to model QCD at finite density, there is no reason to 
use a Lorentz invariant effective action. For instance, at high density the electric glu ons are exp ected to be screened, 
whereas the magnetic ones are Landau damped. In particular, it has been shown by ijSonl Il99^) that at high density 
the magnetic gluon exchange is dominating in the pairing mechanism, which can be simulated assuming Ge <C Gm- 
For the following discussion it is convenient to introduce the quantities 




Ga — ~7\Ge + 3Ga/ ), Gb — ~^(Ge ~ Gm) . 



(4.78) 
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For G E = G 



M 



G, as discussed above one has 



Ga = G, Gb = . 

At zero density we expect Gb = 0, whereas at high density we expect Ge = or Gb/Ga 
relevant physical region for Gb/Ga should be given by 

-!<^<o. 

3 - G A - 



The gap parameters are now defined by 



Ga^s, Ay — Gb^v 



(4.79) 

-1/3. Therefore the 
(4.80) 

(4.81) 



Since the grand potential and the quasi particle energy are determined by the gap, for the Lorentz invariant case, 
Gb = 0, there is no contribution from the vector condensate. For Gb ^ 0, one has to solve two coupled gap equations 
(Alfo rdi p.t q.ZJ . I20T)CT) . The most interesting result found by these authors is about the LOFF window which is modified 
by the presence of the J = 1 gap. The result is shown in Fig. 1171 
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FIG. 17 The figure shows the var iation of the critical points 5/ii/Ao (dotted line) and 8^,2/ Ao (solid line) with Gb/Ga- This 
figure is taken from lAlford et alLlBOOla) . 



The LOFF window closes at Gb/Ga — —3 and increases with increasing Gb/Ga- For Gb/Ga inside the physical 
region (see Eq. (|4.80|l ^ the maximal opening is for Gb = 0. However, inside the physical region the variation of the 
window is rather small. Since 8[ii is essentially defined by the BCS-normal state transition, it is given by Aq/>/2 
independently of the vector condensate. However, the second order critical point 6[i2 is rather sensitive to Gb/Ga, 
the reason being that, for Ga = 0, the J = 1 channel is attractive for Gb > and repulsive for Gb < 0. Therefore 
the stability of the LOFF state is reinforced by the vector condensate in th e region Gb > - 

Let us close this Section considering a different pairing discussed in llDervagin et all . Il992t iPark et all . l200fi 
iRapp et a^.l200lUShuster and Sonll2000(l . It is a quark-hole pairing with non-zero momentum at large baryon density. 
This produces a (qq) condensate varying in space with a wave number 2/x, to b e contrasted with 2|q| ~ 2|<5/x|. This 
state is energetically favorab le only for very large values of the number of colors ijDervagin et all Il992[) . ijPark et all 
l2000UShuster and SorJ.l2000() found that N c should be larger than about 1000. 



E. One-gluon exchange approximation 

The previous results have been obtained in the case of a NJL interaction. In ijLeibovich et a/.U200lj) the case of the 
one-gluon exchange inte r action has been studied. Of course this would be a realistic case only at very high densities 
(Rai agopal and Shusteit l200d|) where, presumably, the CFL phase dominates over the LOFF phase. However the 
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study of different interactions allows to understand the model dependence of the LOFF window. In l)Leibovich et all 
1200 lj) the standard QCD vertex is used in conjunction with the following propagator for the gluon 



pT pL 



where 



Here 



describes the Landau damping and 



D » v ~ p 2 - G(p) + p 2 - F(p) ' (4 - 82) 



pT-X.._P±El pT_pT_r ) 
r vj ~ °V |p|2 ' 00 — r Oi — u ! 

1% = ~9r + ^ - Plu • (4-83) 



G{ P ) = Jm 2 ^ (4.84) 



F(p) = m 2 , (4.85) 
where m 2 is the Meissner mass evaluated for 2 flavors 

u 2 

m 2 =g 2 ^. (4.86) 

The expressions for F(p) and G(p) are obtained in the hard-loop approximation l(Le BelraoUl996l) and evaluated here 
for po <C |p| ~ fi (we recall that fi is the average chemical potential). Solving the gap equation it is found that the 
LOFF window is enlarged of about a factor 10 at average chemical potential fi = 400 MeV. In fact, whereas 6[ii, as 
already noticed, is essentially fixed by the BCS-nor mal state transition at t he value 5fi\ = Ao/y/2, 6fi2 increases in a 
dramatic way. At n = 400 MeV the authors of Ref. l|Leibovich et all l200l|) find 

6fi 2 = 1.24A => 6fi 2 - Sm = 0.55A , (4.87) 

to be compared with the NJL case, where 

<5/i 2 = 0.754A ^ 5^2 - 6m = 0.05A . (4.88) 

At /.t = 10 3 MeV the window is about 60 times larger than the window for the point-like case. In general, by increasing 
fi, <5/i2 increases as well. The interpretation of these results, according to l|Leibovich et oiJ . l200lh . goes as follows. For 
weak coupling the q — q scattering via one-gluon exchange is mostly in the forward direction. This implies that, after 
the scattering has taken place, quarks remain close at the angular position possessed before the scattering, meaning 
that the theory is essentially 1+1 dimensional. In fact, in this case the only possible value for 2|q| is fid — = 26 fi. 
This is not exactly the case in 3+1 dimensions. 
26u- Furthermore, it is known from the 1- 
Il983j) that in the weak coupling limit 6ii 2 /A n — > c 

A similar analysis has been done in Ref. ijGiannakis et al\ . 12002). The results found here are somewhat different 
from the ones discussed before. In particular it is found that at weak coupling: 



ai. in tact, m tms case tne only possioie value tor z\q\ is fid — ^u — zOfi. 
ms. In fact, as it ca n be seen from F ig. 01 2 |q| i s generally bigger than 
-1 dimensional case (jBuzdin and Polonskii Il987l iBuzdin and Tugusheyl 
o — > oo. Both these fea tures have been found in l)Leibovich et all. 2001). 



6fi 2 = 0.968A 5fi 2 - 6fn = 0.26A , (4.89) 

with an enhancement of the window of a modest factor 5 with respect to the point-like interaction. However, the 
evaluation made in this paper consists in an expansion around the tricritical point (called by these authors 6 on ) 
implying, in particul ar, an expansion i n |q|. In order to compare the results of these two groups one should extrapolate 
the results of Ref. ijGiannakis et all l2002|) to zero temperature. It is not evident, at least to us, that thi s can be 
safely made. Of course, also the physical interpretation is different. According to l|Oiannakis et all l2002j) in 3+1 
dimensions, increasing Sfi implies a reduction of the phase space and therefore a smaller gap and a smaller 6fi2- This 
reduction effect, according to these authors, overcomes the enhancement due to the 1+1 dimensional effect discussed 
before. 

In our opinion the case of the one-gluon exchange in the LOFF phase deserves further studies. In fact a sizeable 
increase of the LOFF window would make the LOFF state very interesting as far as the applications to compact 
stellar objects are concerned. 



4G 



F. Mass effects 



In Ref. ijKundu and Raiaeopall l2002|) the combined effect of having two quarks with different chemical potentials 
and one of the two quarks being massive has been studied. In the free case the Fermi momenta are given by (assuming 
that the pair contains an up and a strange quark) 



p s F = V(/x + 5^) 2 - M 2 . 
Assuming both Sfi/n and m s / fi to be much smaller than one, one finds 



IP? 



Pf\ 



M 2 



(4.90) 



(4.91) 



The effect of M s ^ amounts to something more than the simple shift £/i — > 5/j, — M 2 /4/i. In fact, let us recall 
that at M s = t he BCS condensate is not cha nged by 5u as long as 5u < S a-j . However Aq decreases with M s 



see the results in (|Kundu and Raiagopall l2002j) and l|Casalbuoni et all . l2f)02a|i . Furthermore, for M 2 //! 2 <C 1 the 
decreasing is practically linear. This produces corrections in the grand potential of order Aq(0)M 2 . For small values 
of 5[i — M 2 /4/j, we have BCS pairing, whereas for large values there is no pairing and the system is in the normal 
phase. Therefore the BCS-normal transition is M s dependent and occurs for Sfi approximately given by 



M 2 

5 n s - 

* 4 M 



A (M a ) 
V2 



(4.92) 



It can be noted that differently from the case M s = this condition is not symmetric for 5/j, — * — Sfj, and the LOFF 
phase can exist in two different windows in 5(i, above (5^2 > and below (5fj,2 < 5fi\) the BCS region; in any case, 
for M s = one gets back the Clogston-Chandrasekar limit. In order to discuss the size of the window the correct 
variable is 



Sfi 2 (M a ) - StujMs 
A (M S ) 



(4.93) 



At weak coupling (small A o (0)) it is found that the window is essentially the same as for the case M s = 0. Otherwise 
the window generally increases with M s as shown in Fig. Elf° r various values of M s . We have plotted both the cases 
5(i > (left panel) and 5fi < 0. This shows that the LOFF phase is rather robust for M s ^ 0. 
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FIG. 18 In the left (right) panel the width of the LOFF window above (below) the BCS region is reported. The four curves 
correspond to the followi ng values of Ao(0)(MeV) : 10 (solid line), 40 (dotted line), 80 (dashed line), 100 (dash-dotted line). 
This figure is taken from fKundu and Raiaaoval. \ 
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V. PHONON AND GLUON EFFECTIVE LAGRANGIANS 

Translational and rotational invariance are spontaneously broken within a LOFF phase. The energy gap is not 
uniform and actually is expected to vary according to some crystalline structure, as result of the analysis developed 
in the previous Sections. The crystal defined by the space modulation of the gap can fluctuate and its local deforma- 
tions define phonon fields (p 1 - 1 ' that are the Nambu-Goldstone bosons associated to the breaking of the translational 
symmetry. The number of the phonon fields is equal to the number of the broken generators of the translation group . 
The existence of long wavelength oscillations with phonon dispersion law was already noted in (Fulde and Ferrell, 
l^J). More r ecently a n effective lagrangian for phonons in a QCD medium has been developed in ijCasalbuonie^n . 
2002bl l2001all2002cleft and we wish to review it in this Section, dominantly dedicated to the QCD LOFF phase. For 
color superconductivity only the T — > case is physically interesting and we shall consider only this limit. However 
the theory developed in this Section could be extended to T ^ as well as to other physical cases (solid state, nuclear 
physics) . 

Being long wavelength oscillations of the crystalline LOFF structure, the phonons exist only if the quarks of the 
Cooper pair are in the pairing region. This is a portion of the phase space around t he Fermi surface and is forme d by 
a few annular rings, that are likely to be contiguous, according to the discussion in ijBowers and R,aiagopa]ll2002|) . see 
also the discussion in Section Till. CI The effective theory for the phonon fields <j>™ has to display this behavior and 
therefore the phonon-quark coupling must vanish outside the pai ring region. The mathem atical formalization of this 
behavior is rather involved and some approximation is needed. IjCasalbuoni et all l2002blcr) write the phonon-quark 
interaction using the HDET discussed in Section HV.AI They introduce effective velocity dependent fermion fields and 
the lagrangian as a sum of terms, each characterized by its own Fermi velocity v. Also the quark-phonon coupling 
constant becomes velocity-dependent and is proportional to 

ZL^( v .n fc )] . (5.1) 

k v 

Here iifc are the vectors defining the LOFF crystal, R is a parameter and Sji[h(x)] is a function that vanishes outside the 
pairing region. More precisely, it reaches its maximum when the pairing quarks are on the Fermi surface and decreases 
when they leave it. By this approximation an evaluation of the phonon effective lagrangian is possible. In Subsection 
IV. Al we consider the HDET for the inhomogeneous LOFF state and write the quark-phonon lagrangian. Below, we 
discuss two crystalline structures. First we consider the Fulde-Ferrell one plane wave, which is the benchmark case for 
the whole LOFF theory; next we shall examine the cubic st ructure, already studied in Sub section lIII.C.4l because this 
seems the most favored crystalline structure according to ijBowers and Raiagopall 12002). On the basis of symmetry 
arguments one can write down the effective phonon lagrangians for the two cases. This is done in Subsections lV.Bl and 
IV.DI whereas in Subsections IV. CI and IV. El we show how the parameters of the effective lagrangian can be computed 
by the HDET. Let us mention here that the parameter R appearing in (|5.1|) should be fixed by a comparison of the 
gap equation computed in the HDET and the approach discussed in Subsection III. CI for the FF state and in Section 
lllll for generic structures. This comparison has not yet been done and therefore, in the discussion below, we leave R 
as a parameter, even though, in the case of a cubic structure, the requirement that the annular rings are contiguous 
can be used to fix its value. We conclude this Section with a discussion in IV. Fl on the modifications induced by the 
LOFF pairing of quarks on the gluon lagrangian. 

A. Effective lagrangian for the LOFF phase 

Let us begin by writing the gap term in the lagrangian in presence an inhomogeneous condensate. As in Section 
111. CI we write the following formula for the LOFF condensate 

p 

A(r) = ]T A m e 2i *»- r . (5.2) 

m— 1 

We will consider only two cases below: 

a) One plane wave: P = 1 

b) Cubic structure: P = 8. 

In the former case we shall take into account the possibility of having both a J = and a J = 1 condensate as 
discussed above. In the case of the cubic structure we will consider only spin zero condensate; we will take A m = A, 
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real, q m = n m q with n m the eight unit vectors defined in Eq. I|3.39|) . To describe the quark condensate in the case 
of the single plane wave, we consider the lagrangian term: 

C A = £ { l ] + C { 1 ] = e aP ^J a {x)C (a(% + a • nAW^) ^(x) - (L -> R) + h.c. , (5.3) 

which includes both the scalar and the vector condensate. 

We introduce velocity-dependent fields as in Eq. (|4.13(l . with factors exp(ijiijt?j • x), and we take into account only 
the positive energy part that we write as ip Vl -i a for a quark with flavor i and color a; we keep track of the velocities 
of the two quarks that are not opposite in the LOFF phase. We have: 

£a = ~\ ex P{* r • f ( v *> V J' <Jk)Krfv i;i Jz)C (aW^ + a ■ nA^) ^- Vj;J p(x) - (L - R) + h.c. , (5.4) 
where 

f(v ; , vj, q) = 2q-^ 4 v; - Mi v J • ( 5 - 5 ) 

We also define 

/i= , 5(i= . (5.6) 

Since q= 0(A2 SC ) <C /i, the condition 

Pi + P2 = 2q (5.7) 

gives in the /x — > oo limit 

vi +v 2 = O (-] . (5.8) 



Taking into account that P+(— v)Ca fc P+(v) = v k P + (— v)CP+(v) we can rewrite (|5.4|) as follows 

VI, Vj 

These equations can be easily generalized to the case of the face centered cube. We shall discuss this generalization 
below. 

B. One plane wave structure 

Let us rewrite l)5.9|l as follows: 

C A = -~ e 2 ^ e-^+^>r e a ^ } . ia (x)C (A ~ , ;/ - vj • nA^^) ^-v j;j -/j(») — (L R) + h.c. , (5.10) 

Vl,Vj 

There are two sources of space-time symmetry breaking in (|5.1U|I . one arising from the exponential term exp(2ir • q), 
which breaks both translation and rotation invariance, and another one in the vector condensate breaking rotation 
invariance. On the other hand the factor exp(— i/ZiVj — i/^jVj) breaks no space symmetry, since it arises from a 
field redefinition in a lagrangian which was originally invariant. For definiteness' sake let us take the z-axis pointing 
along the direction of q. As a consequence of the breaking of translational invariance along the z-axis, Goldstone's 
theorem predicts the existence of one scalar massless particle, the Nambu-Goldstone Boson (NGB) associated to the 
spontaneous symmetry breaking (SSB). The symmetry breaking associated to the vector condensate is not independent 
of the SSB arising from the exponential term exp(2ir • q), because the direction of q coincides with the direction n of 
the vector condensate. For this reason, while there are in general three phonons associated to the breaking of space 
symmetries here one NGB is sufficient. The argument is sketched in Fig. ^|and follows from the fact that rotations 
and translations are not independent transformations, because the result of a translation plus a rotation is locally 
equivalent to a pure translation. 
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FIG. 19 In the point P the effect of the rotation n — » n' and the effect of the translation r — » r + a tend to compensate each 
other. 



The lagrangian l|5.10[l induces a lattice structure given by parallel planes perpendicular to n: 

nr= — (k =0, ±1, ±2,...) . (5.11) 

We can give the following physical picture of the lattice structure of the LOFF phase: Due to the interaction with the 
medium, the Majorana masses of the red and green quarks oscillate in the direction n, reaching on subsequent planes 
maxima and minima. The NGB is a long wavelength small amplitude variation of the condensate A(r); formally it is 
described by the substitution 



with 



and {4>) =0. We assume 



A(r) = e 2iq& - r A -► e**/ f A , (5.12) 



$ 6 

j = 2q(h + 5n) ■ (r + 5r) = j + 2gn • r (5.13) 



|n + 5n| = 1 , (5.14) 

(Sn) = 0. (5.15) 

Let us introduce the auxiliary functions R and T, 

R = n + <5n, T = 2qR-5r. (5.16) 

The lattice fluctuation tf> describes, in second quantization, the phonon field. Since it must be small, T and R are not 
independent fields and T must depend functionally on R, i.e. T = F\R\, which means 

j = 2 9 R • r + F[R] = G[R, r] . (5.17) 

The solution of this functional relation has the form 



R = h[*] , 



(5.18) 
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where h is a vector built out of the scalar function By this function one can only 3 form the vector V<£ ; therefore 
we get 



R 



(5.19) 



which satisfies (|5.14|) . In terms of the phonon field (f> the vector field R is given (up to the second order terms in 
by the expression 



R = n 



2/<Z 



n(n • V0)] 



f2 q 2 



[3(n- Vc, 



4/V 



(n- V<j 



(5.20) 



We stress that the only dynamical field is </>, $ is an auxiliary field with a non vanishing vacuum expectation value 
($)o = 2q • r; as to Sn, R and <5r, they can all be expressed in terms of <f>. In conclusion, the interaction term with 
the NGB field is contained in 

C mt = -\ e**" Y, e~^+^)- r e am tf } . ia (x)C (M^e.^ - Vj • RA^V*.) V- v s ; j0 (x) (L - R) + h.c. , (5.21) 



where the fields <& and R have been introduced in such a way to reproduce Eq. I|5.10[) in the ground state. At the 
first order in the fields one gets the following three-linear coupling: 



O f 



ir f(vi, vj,q) 



2/ 



A(% - Vj • nA^.l Vj:j , 



1 £ e <r.f(v»,v 4l q) ( _ Vj) . [ V 0- n(n - V0)] AW^e^^^C^^^ - (L - i?) + h.c. (5.22) 



4/g 



We also write down the quadrilinear coupling: 



C^l I - — V e «'f(vi,Vj,q 



2/ 



8/V 



A«e« - Vj • fiA(*>4j Ca/SS^.fa CV-- v j; ,73 
E e ir - f ( v - v ^'(-vj)- [V0-n(n- V^)]AW^ Eo/) 3^ |ta C^ Vjltf 
^ (3(n • V0) 2 - |V0| 2 ) + (vj ■ V0)(n ■ 



e ir-f(vi,Vj,q) 



xA^ 3 £ ;1Q C^_ Vj; ^-(L-,]?) + 



(5.23) 



Through a bosonization procedure one can derive an effective lagrangian for the NGB field. This will be done below. 
For the moment we derive the general properties of the phonon effective lagrangian. It must contain only derivative 
terms. Polynomial terms are indeed forbidden by translation invariance, since 4> is not an invariant field. In order to 
write the kinetic terms it is better to use the auxiliary field $ which behaves as a scalar under both rotations and 
translation. To avoid the presence of polynomial terms in the phonon lagrangian one has to exclude polynomial terms 
in the auxiliary field $ as well; therefore the lagrangian should be constru cted only by derivative te rms. The most 
general invariant lagrangian will contain a tower of space-derivative terms l|Casalbuoni et oil l2001a|) . In fact, since 
(V$) = 2q is not a small quantity, we cannot limit the expansion in the spatial derivatives of $ to any finite order. 
Therefore we write 



P 

C^dpQ) = — 



4> 2 -E c »(i v$ i 2 )' 



(5.24) 



In this lagrangian $ must be thought as a function of the phonon field 



3 In principle there is a second vector, r, on which R could depend linearly, but this possibility is excluded because R is a vector field 
transforming under translations as R(r) — > R'(r') = R(r) 
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Since 



| V$| 2 = V + yn • Vq 



P 



(5.25) 



with similar expression for higher powers. At the lowest order in the derivatives of the phonon field <p we get, neglecting 
a constant term: 



- «f|V||^| 3 -« 2 (4«/V|^ + |V ? 



(5.26) 



where V||^> = n- V0, and vjj, v are constants 



C. Parameters of the phonon effective lagrangian: one plane wave 

In order to derive the parameters of the phonon lagrangian (|5.2(i|) it is useful to make an approximation. We 
assume that <5/i ~ A2 SC <C 5 <C /i. Clearly we cannot take simply the /i — * oo limit in the exponential term 
exp{ir • f(vi, Vj, qk)} in Eq. (|5.9[) : therefore we consider a smeared amplitude as follows: 



lim exp{«r-f(vi, vj, q k )} = lim / dr' exp{ir ' • f (v ; , vj, q k )}.g(r, r ') 

j.L — >OC ' fl — >00 / 

We assume the following smearing function: 

~nq{r k -r' k )~ 



3 sin 



ff (r,r')=ff(r-r')=n 



R 



and we evaluate l|5.27|l in the ^ — > oo limit by taking q along the z— axis, and using the following identity: 

d 3 r'exp{ir' • f}g(r - r') = exp{ir • f} (^) 3 5% (J- 



where 



R . , , 7T 



elsewhere. 



For the components x and y of f we get 



i.e. approximately (for Sfi <C jJt) 



nq 



R 



\V\ + V 2 )x,y\ < — 



From this, in the high density limit, it follows 



vi = -v 2 + 0{6ft//i) . 



(5.27) 



(5.28) 



(5.29) 



(5.30) 



(5.31) 



(5.32) 



(5.33) 



We used already this result in Eq. 1)4. 17|) , in connection with the Riemann-Lebesgue lemma and in Eq. I|5.8|) . A more 
accurate result is as follows. If 6\ and 9 2 are the angles of vi and v 2 with respect to the z-axis one gets 



h = V 2 + 7T 



25n 



tan 6o 



(5.34) 



For the z component we get 



f z = 2qh(cos6 2 ) ■ 



(5.35) 
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where 



1- 



and, neglecting corrections of order Sfi/n, 

h(x) = 1 



cos 8„ 



(5.36) 



(5.37) 



Notice that q — ipo/2 where ipo is the angle depicted in Fig. 0] see also Eq. (|2.117|l . The two factors ir/R arising from 
the x and y components are absorbed into a wave function renormalization of the quark fields, both in the kinetic 
and in the gap terms. As for the z component one remains with the factor 



(5.38) 



in the gap term, whereas for the kinetic term we get a factor of 1. We have assumed exp[i2qhz] = 1 in Eq. (|5.38() 
owing to the presence of the 5r function, that, in the R/t: — > oo limit, enhances the domain of integration where 
h = 0. We will discuss this approximation below. 

Eq. (|5.3l)|l defines a region where Sr ^ 0, i.e. a domain where pairing betwe en the two quarks can occur; it 
correspond to the pairing region in the analysis of ijFulde and Ferreli Il964|) and ((Bowers et a/1 l200lj) . in contrast 
with the blocking region, where 6r = 0. The pairing region intersects the Fermi surface with a 'ring' whose size 
depends on the value of R. As we noticed above, R = oo implies the vanishing of the pairing region and therefore 
one expects R — > oo at the second order phase transition l)Casaibuoni et all 1200 2b). The precise value of R should 
be fixed by the gap equation; since this calculation is still missing, for the purpose of this paper we leave R as a 
parameter. 

In conclusion we can approximate Eq. I|5.9|) as follows 



£a = ~\ E ^IMv • n)] e a ^. ia (x)C (a<%- - v • nA^ 4) ^_ ViW (x) - (L -> R) + h, 

V 

Using the same notations as in Section TlV. CI we can write the effective lagrangian as follows: 

' iTr[f\v ■ DT B ] 



£i+£a=Y1 E X 



At 



v A,B=0 



-A 



AB 



-A AB 
iTr[f\ V ■ D* f B ] 



(L^R). 



Here 



and 



V2 \ot 



Xa 
V2 



(A = 0,...,5) 



v • n A; f ' } a A B 



The matrix Aab vanishes for Aor B — 4or5, while, for A, B — 0, 3, is given by 
with 



TAB 



(5.39) 



(5.40) 



(5.41) 



(5.42) 



(5.43) 
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> \ 
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-l\ 





-1 


















—i 













, GAB = 


















-1 










+i 








\o 












\ 


+1 








/ 



(5.44) 



and 



A^tt 
R 



6 R [h(v ■ n)] 



A 



(v) 
eff 



A^TT 

R 



6 R {h(v ■ n)] 



(5.45) 
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In the present approximation the quark propagator is given by 

/ V-15ac _A AC \ 
Dn r (£) D c 

D AB (£,£") = (2ir) 4 5 4 (£ - £") x ]T 



c 



D CB (£) D cb {£) 
A AC V-£S A c 



(5.46) 



V Dcb{£) D cb {£) J 



where 



d cb (£) = (y-£v-t- AAt ) CB > d cb {£) = (y-iv-i- At a) 



CB 



On the other hand the propagator for the fields x does not contain gap mass terms and is given by 



D(£,£') = (2ir) 4 6 4 {£-£') 



{v-iy 



(V ■ £)~ 1 



(5.47) 



(5.48) 



For the other fields x A , A = 0, • • • , 3, it is useful to go to a representation where AA^ and A T A are diagonal. It is 
accomplished by performing a unitary transformation which transforms the basis \ A into the new basis \ A defined 
by 



X A = Rabx" , 



with 



Rab =t* 



( 1 1 \ 

1 -i 
+i -1 

V i o o -l / 



In the new basis we have 



where 



(AA*)^ = a A 5 AB 
(A^A) AB = & A 5 AB 



ao = a 2 = a 2 =a 3 = (A ( e y f - vnA^) 2 , 
For further reference we also define 



oti = a 3 = a Q = ai = (A ( e y f + V'M^' f ) 



"effJ 



■ W a( s ) -u v .nA (t,) 



» c = (A- - v< A« + vnA<J, A% - vnA- , Kff - ^ W 
In the basis x the 3-point and 4-point couplings 1)5. 22[) and l|5.23[) are written as follows: 

3 / 



o 



-93 - A \ ~B 



v A=0 



-93 ~94 



Here 



!):; 



.94 



id>A (s } f 

" TAB + 0[4>] GAB 



f 

b 2 A (s) 
2/2 



tm; + [ j6[cj>] + Q[<f>] j ,r w 



with 



1 



(5.49) 



(5.50) 



(5.51) 
(5.52) 
(5.53) 

(5.54) 

(5.55) 
(5.56) 
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FIG. 20 Self-energy (a) and tadpole (b) diagrams. 



A 



(v) 
eff 



4fV L 2 



vn 



(3(i 



(v • V0)(n ■ Vc, 



(5.57) 



Terms in 173 and 34 that are proportional to tab arise from the expansion of exp i<f>/ f alone, whereas terms proportional 
to gab get also contribution from the expansion of R in the vector condensate. The effective action for the NGB is 
obtained at the lowest order by the diagrams in Fig. I2U1 



The result of the calculation of the two diagrams at the second order in the momentum expansion is as follows: 

,2 3 



n(p). 



EE 

v C=0 

2 V-pV-p 



16tt 3 / 2 
4a 



A 



4«c 



4a c V ■ I V ■ £ 



D 2 C {£) 



n(p) 



■ 2 3 

;i EE 



(v) - 
eff 



16tt 3 / 2 



d 2 e 

D c (£) 



4ac 



Dm 

lac , 1 
DUl) + D C (£) 



v C=0" 



where 



/Lt c defined in (|5.53() and 



w(p) = p v (p n)(v • n) 



(5.58) 



(5.59) 



(5.60) 



To perform the calculation we will take the limit R — ► 00, when the Sr function becomes the Dirac delta. We 
handle the Sr functions according to the Fermi trick in the Golden Rule; in the numerator, in presence of a product 
of two 5r, we substitute one Sr function with the Dirac delta and for the other one we take 



n8 R [h{x)} tt6r(0) 



1. 



R R 

A similar substitution is performed in the denominator. Moreover we use 

d 2 £ in 



[D c {l)f 



2a 2 c 



Therefore one has 



^6 R [h(vh)} ^ fl [h(vn)] 



1 - 



q v n 



= kRd 



vn 



with 



kR = 



qR 



(5.61) 



(5.62) 



(5.63) 



(5.64) 
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At the second order in the momentum expansion one gets 



Here 



n(p) 



QW(p) 



2tt 2 / 2 



E 5 



A(") 



2 — - > arcsmhr — r 



AW ,*, ^ • v. <* 

— r *(p) 2^ x arcsmh— 

z<? c=o IMc 



AW 



26 



1_log AW ) ( w2 (P)+ v ' fi$ (P)) 



where we have used the result llAlford et a*ll2001bD AW < A^ s ) and 

$(p) = (3p 2 - p 2 ) v • n - 2p • yp z 

From 



£eff(p) 



after averaging over the Fermi velocities we obtain 

1 



£eff - r 



One obtains canonical normalization for the kinetic term provided 



f = 



H 2 k R 
4tt 2 



On the other hand 



v\ = -sin 2 ^ + (1 -3cos 2 ^) (l-log — 



25 \ /AW 



9 



Uij = cos 



In conclusion we get the anisotropic phonon dispersion law 

E(p) = ^v\{p 2 x + p 2 y ) + v 2 p 2 



(5.65) 



(5.66) 



(5.67) 



(5.68) 



(5.69) 



(5.70) 



(5.71) 



(5.72) 



Besides the anisotropy related to v± ^ t>|| , there is another source of anisotropy, due to the fact that p z , the component 
of the momentum perpendicular to the planes l|5.11|l , differently from p x and p v is a quasi momentum and not a real 
momentum. The difference can be better appreciated in coordinate space, where the effective lagrangian reads 



The effective action for the field (f>, S[4>], is obtained by the lagrangian as follows 



S 



f 7T + °° 

dtdxdy — y C(4>(t,x,y,kiT/q) 

J 1 



(5.73) 



(5.74) 



In the action bilinear terms of the type <pk<Pk' with k ^= k' may arise. In the continuum limit these terms correspond 
to derivatives with respect to the z direction. However, in the long distance limit I ^> ir/q, the set of fields (f>k(x,y) 
becomes a function <f)(x,y, z) and the last term in i|5.73[l can be approximated by v 2 {d z (f)) 2 . 
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D. Cubic structure 



The space dependence of the condensate corresponding to this lattice is as follows 

8 

A(r) = A ^2 exp{2ign fc • r} , (5.75) 

k=l 

where the eight unit vectors n k are given in (|3.39|) and 

q = n/a . (5.76) 

In Fig. |^some of the symmetry axes of this cube are shown: they are denoted as C4 (the three 4-fold axes), C3 (the 
four 3- fold axes) and C2 (the six 2-fold axes). 




FIG. 21 Symmetry axes C2, C3 and C4 of the cube. 



To describe the quark condensate we add a term £a completely analogous to 1)5. 10|) . By the same procedure used for 
the plane wave condensate one has 

A 8 7T 

£a = -^J2Y, R 5 RHv-n k )]6 rj e a ^l. ia {x)Ci)- v , jl3 {x) -(L-fii) + h.c. (5.77) 

fc=l v 

£0 + £1 + £a is still given by Eq. (|5.40|) but now 

fc=i 

the quark propagator is given by 1)5. 46[) with A e / / given by (|5.78|) . 

An interesting point should be noted. This equation shows that the pairing region for the cubic LOFF condensate is 
formed by eight distinct rings; each ring is associ ated to one vertex of the cube and has as its symmetry axis one of the 
threefold axes C3. According to the analysis of l|Bowers and Raiaeopall l2002[) . the LOFF vacuum state corresponds 
to a situation where these domains have at most one common point. Given the symmetry of the cubic structure we 
can limit the analysis to one pair of rings, for example those associated to the vertices ni , 115 . The common point 
between these two rings lies on the axis C2 and has v = 1/V2(1, 1,0). Since it must also belong to the boundary of 
the two pairing regions we have the condition: 



Mvnt)] = — , 



(5.79) 
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which implies 

R = j^- r . (5.80) 

2h f\/273J 

Using Eq. (|5.36|) one gets 

(5.81) 

The condensate (|5.75ll breaks both translations and rotations. It is however invariant under the discrete group Oh, 
the symmetry group of the cube. This can be seen by noticing that the condensate is invariant under the following 
coordinate transformations 



R 1 : 


Xi - 


-> Xi, 


X2 — 


X3, X 3 -> - 


■a;2, 


i? 2 : 


Xi - 




X2 


-> x 2 , x 3 - 


► Xi, 


R 3 : 


Xi - 


-> x 2 , 


x 2 -» 


-Xt, x 3 - 




/ : 


Xi - 


-> -Xi, 


X2 


-> -x 2 , £3 





a*, (5.82) 

that is rotations of 7r/2 around the coordinate axes, and the inversion with respect to the origin. Since the group Oh 
is generated by the previous 4 elements the invariance follows at once. 

The crystal defined by the condensate H5.75JI can fluctuate and its local deformations define three phonon fields 
0w that are the Nambu-Goldstone bosons associated to the breaking of the translational symmetry. They can be 
formally introduced following the same procedure discussed for the single plane wave case. One effects the substitution 
in f5~75) 

2^^ = ^ + ^M, (5.83) 
fa f 

where the three auxiliary scalar fields <&W satisfy 

2tt 



, ,, Xi, (5.84) 

/ /o a 

whereas for the phonon fields one has 

(0«(aO)o = O. (5.85) 
We have therefore three fluctuating fields 4> k { k2 k3 for any elementary cube defined by discrete coordinates 

k\1T k 2 TT k 3 TT 

Xk ± = , Vk 2 = , z k3 = , (5.86) 

I q q 

i.e. 

4>f lk2k3 ^<j> (€> {t,x kl ,y k2 ,z k3 ) . (5.87) 

The interaction term with the NGB fields is therefore given by an equation similar to (|5.74() : 

- / \ 3 +°° 8 
S m t = -jdt(-\ Yl EE A ex P{* VfcTfa k 3 //} W^l, iot C 1>- v; jfi — (L —> R) + h.c. , (5.88) 

* ' ki,k2,k 3 = — 00 v m=l 

where 

*L fc 3=E4 m) ^ fc2fe3 (5.89) 

i=l 

and the eight vectors e^ m ' are given by 

(e 4 (m) ) = 73n m . (5.90) 
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The complete effective action for the NGB fields cjr- 1 ' will be of the form 

dt(-J Y, C{^{t,k^/q,k 2 T:/q, k 3 n/q)) . (5.91) 

In the low energy limit, i.e. for wavelengths much longer than the lattice spacing ~ 1/q, the fields <f>j® ks vary 
almost continuously and can be imagined as continuous functions of three space variables x, y and z. 
The coupling of the quark fields to the NGB fields generated by the condensate will be written as 

AV> T CV> ex P {*(ei$ (1) + £2* (2) + e 3 $ (3) )} , (5.92) 

£i=± 

making the theory invariant under translations and rotations. These symmetries are broken spontaneously in the 
vacuum defined by Eq. I|5.84|l . In order to write down the effective lagrangian for the phonon fields <f>^ it is useful to 
start with the effective lagrangian for the auxiliary fields $W that has to enjoy the following symmetries: rotational 
and translational invariance; Oh symmetry on the fields The latter requirement follows from the invariance of 

the coupling l|5.92|l under the group Oh acting upon <S>( 1 \ The phonon fields <f>( l \x) and the coordinates x l must 
transform under the diagonal discrete group obtained from the direct product of the rotation group acting over the 
coordinates and the Oh group acting over $W(x). This is indeed the symmetry left after the breaking of translational 
and rotational invariance. The most general low-energy effective lagrangian displaying these symmetries is 

£ = Y E (* (0 ) 2 + >Cs(/2(V$W),/4(V$W),/6(V$W)) , (5.93) 

j=l,2,3 

where 

h(Xi) = X\ + Xl + Xl h{Xi) = X\XI + XIX\ + Xjxl h(Xi) = XlX 2 2 Xl (5.94) 

are the three basic symmetric fu nctions of three variables . At the lowest order in the fields and at the second 
order in the derivatives one gets l)Casalbuoni et all 12002^ 

C=\ E (W"? E l^ W | 2 -^ E (d^f-c £ d^d^K (5.95) 

i=l,2,3 i=l,2,3 i=l,2, 3 i<j=l,2,3 

which depends on three arbitrary parameters. 

E. Parameters of the phonon effective lagrangian: cubic crystal 

The parameters a,b,c appearing in (|5.95|l are computed by a method similar to the one used in Section IV. CI One 
puts 

3 

^ m \t,r) = Y J ^ n U [l) M , (5.96) 

i=i 



which allows to write the 3-point and the 4-point couplings as follows: 



.93 -94 



£a+A = £Ex At ( _ °_ "^-^W, (5-97) 



v A=0 

Here 



7rA „ iif^ 



93 = E -fttniKv ■ 7 T~AB ) 

m=l ■* 
8 



.94 = - E V W v ■ ^m)]^-r AB , (5.98) 

m=l J 
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to be compared with Eqs. (|5.55|) and (|5.56(l valid for the one-plane wave form of the condensate (we have here 
neglected the vector condensate). To perform the calculation one employs the propagator given in Eq. I|5.4(j|l with 
A e f f given in ((5.78(1 and the interaction vertices in ((5.97(1 . The result of the calculation of the two diagrams in Fig. 
1201 at the second order in the momentum expansion is 



4x4 



16tt 3 / 



' v m,fe=l ^ ' 



dH 



D{£)D{£ + p) 



2A 2 eff + v ■ ev ■ (e+p) + v ■ ev ■ (i+p) 



(5.99) 



^ff( P ) tad = *^E E / m^Wv-U^ 1 ) 8 , (5-100) 

J v m=l J ^ ' 

where 

D(£)=£ 2 -q-A 2 eff + ie, (5.101) 

and, analogously to 1(5 .45(1 . 

A 8 

fc=l 

From 1)5.991 and (|5.1U1|I one can easily control that the Goldstone theorem is satisfied and the phonons are massless. 
As a matter of fact one has 

. 4 x 4fi 2 it A 

Lrnass = £ e ff[0)s.e. + £eff{0)tad = I ^3 p x 

x E/^j[- E ^^[/ 1 (v.n m )]^(™)^[Mv-n fc )]^)+A e// ^5 fl [Mv-n ro )](^W) 2 ] .(5.103) 

v ^ ^ ' m,fc=l m=l 

In the double sum in the r.h.s. of Eq. ((5.103(1 only the terms with m = k survive and one immediately verifies the 
validity of the Goldstone's theorem, i.e. the vanishing of 1(5.103(1 . Notice that in this approximation the masses of the 
Goldstone bosons vanish because the pairin g regions are not overlapping , signaling that when they do overlap one is 
not at the minimum of the free-energy, see l|Bowers and Raiaeopall l2002|) . 

At the second order in the momentum expansion one has 

%W=^i: £ \ (^) 2 ^[fe(v-n ro )](^W)^[Mv-n fe )](^( fe )) j d 2 l 2Alf ^ 3 V ' P ■ (5.104) 
Using the result 

d 2 £ in 



[D(£)}3 2Ai ff ' 

and the absence of off-diagonal terms in the double sum, we get the effective lagrangian in the form 



(5.105) 



*«M - 4p ? (£) 3 1 (v . tWH v. p ) ^ . „ 106) 

To perform the calculation one can exploit the large value found per 7?, and use the same approximations of Section 
IV. Gl The sum over k in ((5.106(1 gives 

D 8 9 



^28 R [h(yn k )}^6 R [h(v-n k )}^ -> J ^^(v-n t )] (ip 

fe=l k=l 
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= *y s 

k=l 



1 - 



S/i 



gv • n fe 



^ 2 = %"*E* 

fe=i 



v • n fe 

9 



with 



Therefore one gets 



qR 



^ i,i=lfc=l v ^ y J 

The integration over the Fermi velocities requires special attention. We use the result 

8 



E 

k=l 



M (k) _ or . 



(5.107) 



(5.108) 



(5.109) 



(5.110) 



this fixes the constant multiplying the time derivative term in the effective lagrangian at the value (taking into account 



8fi 2 k R 
2 x 4tt 2 / 2 ' 

Therefore one obtains canonical normalization for the kinetic term provided 



4tt 2 



J 2 = 



(5.111) 



(5.112) 



The parameters a, b, c of the effective lagrangian l|5.95|l can be now evaluated and one finds ijCasalbuoni et q/ll2002b|) : 



= l (rP 2 4» 2 ^ 6 (i)2 3 COS 6g 1 
2V ^ 12 ^ 6 



i.e., comparing with Eq. 



a =l2' 5 = ° 



i<j=l,2,3 



3 COS 2 q - 1 
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(5.113) 



(5.114) 



F. Gluon dynamics in the LOFF phase 

1. One plane wave structure 

In this Section and in the subsequent one we wish to derive the effective lagrangian for the gluons of the unbroken 
SU(2) C subgroup of the two-flavor LOFF phase. To begin with we assume the crystal structure given by a plane wave 
and we neglect the vector condensate, so that we write 



A eff = -^S R [h(v-h)] 



(5.115) 



The effective action allows the evaluation of the one loop diagrams with two external gluon lines and internal quark 
lines similar to those in Fig. [3U| If one writes 



then the Meissner mass vanishes 



Ktip) = rC(o) + flC(p) ; 
ir^(o) = o 



(5.116) 
(5.117) 



61 



and the Debye screening mass is non-vanishing 



5M L . cos Og - cos b ,. 11fi > 
m D = — yH , (5.118) 

where cos# a and cos6*b (— 1 < cos6* a < cos0(, < 1) are the solutions of the equation 

|M«*0)| = ^. (5-119) 

Next consider 511^ (p). The only non vanishing contribution to 5H^ b (p) comes from the pairin g region, i.e. where 
A e /r 7^ 0. In the approximation of small momenta (\p | <C A) one finds IjCasalbuoni et adl2002q) 

r v^(g ■ p) 2 - gn^(7 ■ g g ■ p) + y ~ g 

v; pairing 

That is 

ryrOO / \ r ffV ^ c ffV^ 2 f d cos 0d(f> Vi v 3 ,,.„.> 

-OT ab (p) = <U^- a" = S *b » A2 , / —77- ^aP«Pi, (5-121) 

37F v; P ^„ S A e// 3A% firing 8tT (<5fl[/l(cOS 6)]) 

where v = (sin#cos</>, sin0sin0, cos#). The integration domain is defined by cos# a < cos 6* < cos 9b- Therefore we 
get 

-5B°M = 6 ab k(f(R)pl+g(R)pj) , (5.122) 
where k is given by (Ao the homogeneous condensate): 

and 

f(R) = J f dcos6»(l-cos 2 6») , (5.124) 

J pairing 

9( R ) = \ ( dcos6cos 2 9 (5.125) 



2 



pairing 



are functions of the parameter R and are reported in Fig. 1221 

It is interesting to note the anisotropy of the polarization tensor exhibited by these results. One has always g > /; 
for large R, and neglecting 5p/p corrections, one finds approximately 



Let us finally write down the remaining components of the polarization tensor. From (|5.12U|I we get 

" ^lip) = E = k Pl (nmWji+faW+giRWja) (5.127) 

v; pairing 

and 

~n° a l( P ) = kp pJ (fiRXSnSji + SaSrf+gWSaSja) ■ (5-128) 
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FIG. 22 Plots of the functions f(R) and g(R). 

These results complete the analysis of the LOFF model in the one plane wave approximation. From we get the 
dispersion law for the gluons at small momenta. The lagrangian at one loop is 4 

C = --F» v F"--n t "?A"A b v . (5.129) 

(sum over the repeated color indices a,b — 1,2,3). Introducing the fields E? = and B" = iEijkF^, and using 
(|5.122|l . (|5.127|) and l|5.128|l we can rewrite the lagrangian (|5.129|l as follows 

C = \ (ey E?Ef ~ B?B?) + ±m 2 D A a A° a , (5.130) 

where 

/l + kf(R) \ 

eij = 1 + kf(R) . (5.131) 

V l + kg(R)J 

This means that the medium has a non-isotropic dielectric tensor e and a magnetic permeability A = 1. These results 
have been obtained taking the total momentum of the Cooper pairs along the z direction. Therefore we distinguish 
the dielectric constant along the z axis, which is 

e|| = 1 + kg{R) , (5.132) 

and the dielectric constant in the plane perpendicular to the z axis 

ex = 1 + kf(R) . (5.133) 

This means that the gluon speed in the medium depends on the direction of propagation of the gluon; along the z 
axis the gluon velocity is 

' (5.134) 



" VkgjR) 

while for gluons which propagate in the x — y plane we have 

1 



VkJW) (5 ' 135) 



4 We do not include here the 3 and 4-gluon vertices that however can be handled as in i'Casalbuo ni et all. I2002dh . with the result that 
the local gauge invariance of the one-loop lagrangian is satisfied. 
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and in the limit of large R, and neglecting 5fi/ fi corrections, 

U|l — > — y= tan 9 q v± . (5.136) 

with cos q defined in Eq. I|5.37|) . 
2. Cubic structure 

The condensate in this case is given by Eq. I|5.75[l . so that we will use the results of Section IV. Al with A e ff given 
by (|5.78|) . The calculations are similar to the previous case and, similarly, the SU(2) C gluons have vanishing Meissner 
mass and exhibit partial Debye screening. However the dispersion law of the gluons is different. 

As a matter of fact we write the one loop lagrangian for the SU(2) C gluons as 

C = ±(E?E?-B?B?)+6£ , (5.137) 



with 



5C = -\KlA1Al . (5.138) 



In the approximation \p\ <C A, SH^ is again given by Eq. (|5.12U|) . but now A e ff is given by (|5.78|l . One gets 

,2„2 r J^OJA „,„,. „2,,2 



«£. J**} ^^L + ASAUJ-^J <!«.«. (5.139) 

u " Jpairing °" eff J blocking 

Evaluating the integrals one finds 



C = l -{e l3 E?Ej - BfBf) + \M D 2 A a A a , (5.140) 



with the tensor e lJ given by 



and 



c 



y = Sij [1 + kt(R)} . (5.141) 



M D = wJl + 8 C0Sda - C ° Sdb . (5.142) 

7T V 2 

where cos^q^ are solutions of Eq. H5.119|) . The tensor e u is isotropic. This result can be easily explained noticing 
that the lagrangian should be a quadratic function of the field strengths and should also satisfy th e cubic symmetry . 
There fore it must be constructed by the invariants ^(^i) and hiBf) that are isotropic. As shown in ijCasalbuoni et all 
l2002d) t(R) is given by 

t(R) = ^[2f(R)+g(R)}. (5.143) 

It should also be noted that the values of the parameter R for the cube and the plane wave can be different. A plot 
of the function t(R) is in Fig. EH Even if the crystalline structure is not isotropic, the dielectric properties of the 
medium will be isotropic and the velocity of propagation of the gluons will be the same in all the directions. 



VI. INHOMOGENEOUS SUPERCONDUCTIVITY IN CONDENSED MATTER, NUCLEAR PHYSICS AND 
ASTROPHYSICS 

As observed in the introduction, the main focus of this review is on the theoretical methods rather than phcnomcno- 
logical consequences. However, for completeness, in this Section we give a review of the different approaches developed 
so far to detect the inhomogeneous phase in superconductors. The LOFF is expected to be ubiquitous, therefore one 
might expect to find it in completely different physical systems. For obvious reasons research in solid state physics is 



G4 



t(R) 



1.0 




0.8 



0.6 



0.4 



0.2 



10 



15 



20 



25 



30 
R 



FIG. 23 Plot of the function t(R). 

much more advanced and indeed signals of formation of the LOFF phase have been reported in the literature. In the 
first part of this Section we review them. In order to produce the effective exchange interaction of Eq. (|2.f|) a suffi- 
ciently high magnetic field is needed to produce an appreciable difference between the chemical potentials. This could 
be done in type I superconductors, but the required fields are likely to destroy superconductivity altogether. This 
issue is discussed in Subsection I VI. Al To overcome this difficulty, type II superconductors must be used; moreover 
they should be free of impurities and have large electron mean free path; the needed requirements and the associated 
phenomenology are discussed in Subsection IVI.BI and IVI.CI Another way to overcome the effects of high magnetic 
fields that are detrimental to electron superconductivity is to use layered superconductors and magnetic fields parallel 
to the layers. Superconductors of this type are rather different from the ones considered in the original LOFF papers; 
in particular organic superconductors are compounds with these features and are therefore good candidates. They 
will be discussed in Subsection IVI.DI while in Subsection I VI . El we briefly discuss future possible developments in the 
area of atomic physics. The final part of this Section is devoted to phenomenological implications of the LOFF phase 
in nuclear physics (Subsection IVI.F|) and QCD ( Subsections I VI . CI and I VI ,H|l . In particular, in this last Subsection 
we discuss a possible role of the QCD LOFF phase in the explanation of a peculiar phenomenon of pulsars, i.e. the 
periodic glitches in their angular velocity. 

A. Type I superconductors 

In the original LOFF papers l|Larkin and Ovchinnikovlll964h . l|Fulde and Ferreil Il964'l the difference in chemical 
potentials between spin up and spin down electrons arises from an interaction of a magnetic field with the electron 
magnetic dipole moments. The magnetic field can hardly be the external field H ext ., which exerts a stronger influence 
on the orbital motion than on the electron spin. Therefore the inhomogeneous phase was thought to arise in nonmag- 
netic metals in presence of paramagnetic impurities. Under an external field the host impurities display ferromagnetic 
alignment; decreasing the temperature the material becomes a superconductor while the ferromagnetic alignment 
persists, leading to a constant self-consistent exchange field, proportional to the average spin of the impurities. This 
field is at the origin of the modulated order parameter. The value S/ii « a bove which LOFF phase can exist 

corre sponds to a critical value of the magnetic field that can be derived as follows l)Ck>gstonlll962|) . l)Chandrase khar. 
1962). The susceptibility of an electron gas in the normal phase at T = is 



where p = gp 2 F / (2it 2 vf) is the density of states at the Fermi surface, \xb is the Bohr magneton and g = 2 is the 
electron degeneracy factor. On the other hand the susceptibility in the superconducting phase at T — vanishes: 
Xs — 0, because to polarize the superconductor one has to break the Cooper pair, which costs energy. The free energy 
per unit volume f s of the superconductor, in absence of paramagnetic effects, is: 



Xn 




(6.1) 



fs fl 



1 1 




(T) 



(6.2) 
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where f n is the free energy of the normal phase and H C (T) is the critical field. Including the Pauli paramagnetism 
implies adding the hamiltonian a term 



while 1)6. 2)1 becomes 



(iB^o-'Hf; (6.3) 



A = / „ + fc^_SCT). (6 . 4) 
Z on 



Therefore the BCS superconductivity will survive at T = for magnetic fields satisfying 



H<^&1^H P (0). (6.5) 
Now H^(0)/8n = pAp/4 and therefore the Pauli limiting field at T = is 

iJ P (0) = 4=— ■ (6.6) 
V2 Ms 

The identification <5/ii = iipHpiQ) arises by the comparison between 1)6.3(1 . 1)6. 6JI and 1)2.1(1 . 

For a type I superconductor it is difficult to reach the Pauli limit 1)6.6)) because, while -ffp(O) is typically of the 
order of 300 K Oe, H c (0) is of the order 1 K Oe. Therefore superconductivity will be broken by the magnetic field 
well before the Clogston-Chandrasekhar limit is achieved. This i mplies the LOFF phase is unlikely to be produced by 
these materials and one has to turn to type II superconductors l|Saint- James et q/1 Il969fl because for some of these 
superconductors the upper critical field H C 2 can be very high. 



B. "Clean" and strongly type II superconductors 

To evaluate the possibility of the LOFF sta t e one has to take into account not only Pauli paramagnetism of the 
electrons but also the orb ital effects I GinzbureL Il957|) . Before doing that, let us first distinguish between "clean" and 



1959j). One calls "clean" the superconductors in which the electron mean free 



'dirty" superconductors ijAndersonl 
path I is much larger than the superconducting coherence length £o : 

I » £o , (6.7) 

while "dirty" superconductors are characterized by the opposite condition I <C £o- In "clean" superconductors electrons 
at the Fermi surface move with velocity vj?, while in "dirty" superconductors the electron motion is described by a 
diffusion equation. "Dirty" su perconducto r s are characte r ized b y the presence of impurities, which can narrow and 
even destroy the LOFF state ijAslamazovt Il969t iTa kadal Il 97dl). Therefo re materials with small /, e.g. PbMogSs, 
should not display the LOFF phase, see e.g. ijDecroux and Fischer! Il982j) . On the other hand the so called heavy- 
fermion superc onductors are favored: these materials are indeed characterized by the small Fermi velocity of their 
quasi particles ijRauchschwalbel Il982|) : since 

& = 5? , (6-8) 

for small enough vf the condition 1)6.7)1 is satisfied. For example in the heavy fermion compound UPd2Al3 the 
superconducting coherence length is £o ~ 85 A, much smaller than the electronic free mean path I = 700 A; therefore 
it can be considered as a "clean" superconductor. 

To evaluate the possibility of the LOFF state in type II superconductors one has to take into account not 
only Pauli paramagnetism of the electrons, but also the orbital effects. This ana lysis was first performe d by 
i|Gruenberg and Guntherl |l966). These authors followed the variational method of ijWerthamer et all Il96rj) . by 
making an ansatz for the condensate. In general there is a competition between the orbital and the paramagnetic 
effect, the former trying to organize a structure of Abrikosov vortices and the latter a periodic LOFF structure; 
therefore the orbital effect reduces the possibility of the LOFF state that can exist only for sufficiently high H C 2. The 
quantitative criterion at T = for clean superconductors with isotropic dispersion law is as follows. The LOFF state 
can persist in a type II superconductor provided 



(6.9) 
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Here a is the parameter first introduced in l)Makil . ll964|) . H C 2(0) is the Gor'kov upper critical field at T — in absence 
of paramagnetic effects llGorWll960h and H P (0) is the Pauli limiting field defined in 1)6.60 . 

In conclusion good experimental conditions to observe the LOFF state should be provided by a clean supercon- 
ductor with a large a value. These features are not easily found in the most common superconductors and therefore 
experimental investigations consider unconventional superconductors, e.g. heavy-fermion, organic or high T c super- 
conductors. As a matter of fact many of these materials are strongly type II superconductors, which means that 
the condition l|6.9|l can be satisfied. Moreover they have often a layered structure, which implies that, applying the 
magnetic field parallel to the layers, the orbital effect is minimum and the Zeeman effect, on which the LOFF phase 
is based, is dominant. 

The condition of being very clean and simultaneously strong type II superconductors should be more easily realized 
in d~ wave superconductors like high T c cuprate superconductors and organic superconductors like k-(ET) 2 or A- 
(ET)2 salts. They will be discussed in more detail below; suffice it here to say that in d— wave sup erconductors the 
regio n of the LOFF phase is much more extended than in s— w ave superconductors. The analysis of l)Maki and Wonl 
1996), where this conclusion was drawn, has been extended in dWou and Makl feo^ to the calculation of the LOFF 
free energy, specific heat and magnetic susceptibility; in particular for these layered d— wave superconductors the 
energetically favored structure at T = is found to be 

A(x,y) oc cosqx + cos qy . (6.10) 

Oth er materials where t he possible existence of LOFF phase has been investigated are ferromagnetic metals or 
alloys l|Pickett et all Il999l) , llDvugaev et all 1200 ll). The study of the possible coexistence of ferromagnetism and 
superconductivity was initiated bv ijGinzbur el 119571) who noted that, though the two orderings can in principle 
coexist, their simultaneous presence is practically impossible under ordinary conditions. As a matter of fact the 
presence in ferromagnets of a spontaneous magnetization Mo produces at T = an internal magnetic induction 
Bq = AttMq even in absence of external magnetic field. For superconductivity to exist, Bq should be smaller than the 
lower critical field at T — 0, in absence of ferromagnetism: 

B o <H° cl (0). (6.11) 

However, the induction B at T = is of the order of 10 K Oe (e.g. 22, 18.5, 6.4, 24.8 KOe respectively for Fe, 
Co, Ni, Gd), while the critical field is in general much smaller, of the order of a few KOe or less. Superconductivity 
of ferromagnets is therefore difficult unless special conditions render the condition (|6.11l) possible. They might be, 
for example, a reduced size of the sample, with dimensions of the order of the penetration depth. The formation of 
the vortex phase in type II supercond uctors , howe ver, screens locally the internal magnetic induction, and allows to 
avoid Ginzburg's negative conclusion llKrevl Il973l) As a matter of fact superconductivity was recen tly reported in 
the ferromagnetic alloy RuSr^GdCu^Os (|Bernhardl Il999t lHadiievl. Il999t iPringlel Il999t iTallonl |l999). This layered 
material becomes first ferromagnet at T = 132 K; superconductivity appe ars at T = 35— 40 K and finally, at T — 2.6 K 
Gd ions acquire an antiferromagnetic order. The theoretical study of ijPickett et all Il999j) confirms these reports, 
but suggests that the superconducting phase is of the LOFF type, because the coupling between ferromagnetism and 
superconducting layers appears to be suffici ently weak to permit superconductivity, but strong enough to require the 
inhomogeneous phase. In a similar context l)Dvugaev et adl200l|) consider the possibility of creating the LOFF phase 
using ferromagnetic materials instead of nonmagnetic bulks with paramagnetic impurities as in the original LOFF 
papers. Since the impurities create not only an exchange interaction, but also an electromagnetic interaction, using 
nuclear ferromagnetism, as they propose, would reduce the latter, since the effective field would be proportional to the 
nuclear magneton and not to the Bohr magneton. They show that in some metals, e.g. Rh, W, the BCS condensate 
imbedded in a matrix of ferromagnetically ordered nuclear spins should manifest the LOFF phase. 

All the proposals we have discussed so far are rather different from the one discussed in the original LOFF papers. 
An extension of the LOFF analysis to these materials and unconventional superconductors is beyond the scope of 
the present review. We will therefore limit our presentation to a brief survey of the experimental results, referring 
the interested reader to the specialized literature l) Agterberg and Yand.l200ll|Gege bwartL]T996llMurthv and Shankait 
119951 ISamohkinL Il997t IShimaharal Il998bl: IShimahara et all Il996t ISvmingtorl Il999t lYand l200l[) . 



C. Heavy fermion superconductors 



Th e first exper imental inves tigations on the LOFF phas e used heavy-fermion compounds such as CePai2 l|Huxlevl 
1993), UPd 2 Ai3 llGloosi Il993h and UBei 3 llThomasl Tl996). For all these materials the conditions for the formation 
of the LOFF state are met. For example CeRu2 is in a metallurgically clean state; moreover it exhibits extreme type 
II behavior, because the Ginzburg-Landau parameter, which discriminates between the two type of superconductors 
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llSaint- James et all Il969ft . has the value k = 16. As another example, the compound UPchjAL; used by ijGloosl 
1993) is characterized by a rather large value of the parameter a in 1)6. 9|) i.e. a = 2.4, while als o being a very clean 
superconductor. To quote another result, in the analysis of a high quality single crystal of UBei3 l|Thomas| . ll996ft . the 
very high value H C 2(0) ~ 140 KOe was reac he d. All these exp erimental results are however inconclusive. In a critical 
analysis of the experiment of ijGloosl fl 993ft . l)Normari 1 1993ft shows that the computed Gor'kov upper critical field 
does not correspond to the experimental resu lts reported there; for fu rther analysi s of the compound UPd2Al3 see 
l|Yin and MakiL 1199^1 and (jSchimanskil. Il994ft . In the case of Ce Ru?. llTenval 1 1999ft shows that the observed effects 
can be explained by flux pinning mechanisms involving disorder. l|Modlerl . ll996ft makes a comparative study of high 
quality single crystals of UPd2Al3 and CeRu2 in the mixed state. The order parameter exhibits a periodic array of 
nodal planes perpendicular to the Abrikosov vortex lines. In the mixed state the pinning force is very weak; however 
the authors find, for H > 10 KOe and T < 0.9 T c , a first order transition to a state characterized by strong pinning, 
which might be interpreted as the formation of a LOFF state. The mechanism by which Abrikosov vortex lines in 
type II superconductors are pinned to the vortex cores is similar to the one that pins vortex lines to non-supcrfluid 
neutrons in a rotating superfluid within neutron stars. It will be explained in more detail in Subsection I VI. HI in 
connection to a possible role of the QCD LOFF state in the physics of pulsars. 



D. Two-dimensional, quasi-two-dimensional and organic superconductors 



As we already mentioned, the paramagnetic effect can dominate if the superconducting bulk has a layered structure 
and the magnetic field acts parallel to it, because in this case the orbital upper critical field can be extremely high and 
the breaking due to the spin interaction is m ost sig nifica nt. The imp ortance for two-dimensionality (2D) to favor the 
LOFF state was first observed in l|Bulaevskiil Il973ft and l)Bulaevskiil Il974ft whe re both the orbital and t he spin effect 
were taken into account and the upper critical field H C 2(T) was calculated; in ijBuzdin and Kulicl Il984ft the analysis 
was carrie d out near the tricrit i cal point. For the same reason also quasi-one-dimensional (Q1D) com pounds were 
discussed llBuzdin and PolonskiL Il987t iBuzdin and Tueushevl Il983t IPumul Il995l fPupuis et all ll993L even though 
the results of (Shi maharai Il998aft indicate that the 2D structures are favored in comparison the ID ones. 

Th ese results were generalized to arbitrary temperature and d— wave superconductivity in ijShimahara and Raineit 
11997ft . The main result of this paper is that the critical field curve H C 2 (T) is non monotonic and consists of different 
pieces corresponding to different Landau levels, characterized by n > 0. On the contrary, the Ginzburg-Landau 
theory would predi ct the pair wave function to be in the lowest energy Landau level, with n = at H C 2- The paper 
( Shi maharai Il998a|) studies the most favored structure for a 2D LOFF crystal with a cylindrical Fermi surface. First, 
the author finds that in general the 2D structures are favored over the ID ones; second, it finds that the favored 
crystalline structure changes with T. For s— wave the results are as follows: at high temperature the antipodal pair 
condensate 



A(r) cx 2cosq- r (6.12) 

is favored. This was the result found by ijLarkin and Ovchinnikovl ll964) in 3D at T = 0. Decreasing the temperature 
other structures become favored: first the triangle, then the square and finally, at low temperatures, the hexagon. 
For d— wave pairing at high temperature again (|6.12|) is favored, while at small temperature the square dominates; on 
the other han d at intermediate temp eratures the phase transition s hould be first order. The result at T — has been 
confirmed by IjWon and Makil 12002ft . see Eq. (|6.10|l . A s shown in |Lebed'L Il986j) . the quasi-2D superconductors can 
be treated as essentially 2D and therefore the results of ijShimaharal Il998af) should hold also for quasi-2D compounds 
pr ovided the externa l field is sufficiently strong and is kept parallel to the superconducting layer. 

l|Klein et all 12000ft consider a layered superconductor in a magnetic field of arbitrary orien tation with respec t 
to the conducting plan e. The calculation is based on the quasi-classical Eilenberger equations ijEilenbergeit Il968). 
((Alexander et aZ.fll985ft and allows to elucidate the structure of the stable states below i?2c minimizing the free en- 
ergy. The stable states are neither pure LOFF states nor pure Abrikosov vortex states, but are two-dimensional 
periodic structures or quasi- one-dimensional structures where LOFF domains are separated by vortex chains. 
( Barz vkin and Gor'kovt 12002ft address 2D surface superconductivity in presence of intense magnetic fields parallel 
to the surface. The spin-orbit interaction at the surface changes the properties of the LOFF state; the authors find 
that strong spin-orbit interactions significantly broadens the range of parameters where the LOFF phase can exist 
and produces periodic superconducting stripes running along the field direction on the surface. 

Organic superconductors are good candidates for the formation of the LOFF state for the reasons mentioned above: 
i) They have narrow electron bands and therefore they are in principle clean type II superconductors; ii) due to 
their low dimensionality the orbital pair-breaking effect is suppressed fo r magnetic fields parallel to the layers they 
form. For these reasons they have been d i scussed by several authors, e.g. llBurkhardt and Raineit 1 19941 iDupuisl . Il995i 
iDupuis et al or'kov and Lebedt Il987t iLebedl Il986t IShimaharalTl994L Il997ft . It is interesting to note that 
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in general to detect the transition from the homogeneous BCS phase to the LOFF phase thermodynamic signatures 
are chosen. This however can gi ve ambiguous results since the signatures can be produced by phase transitions of 
different nature. Therefore Ref. ((Yang and Agterberel l2000l) proposes the use of the Josephson effect. According to 
this analysis, at the Josephson junction between two superconductors, one in the BCS and the other in the LOFF 
phase, the Josephs on current is suppressed. 

As discussed in (|Shimahara and Rained . Tl99^ . the upturn of the upper critical field d 2 H c /dT 2 > is a common 
feature in the LOFF state in quasi-2D systems and is due a Fermi surface effect. Investigation s on the sensitivity of the 
LOFF state to the shape of the Fermi surface are in IjAoi et a/.lll974ilShimah ara. 19941 Il997|) . This upturn and a first 
order transition below the critical field have been observed in the organic compound K-fBE DT-TTF^CufNCS)?- This 
quasi-two-dimensional (Q2D) organic superconductor is examined by a number o f autho rs ijHouzet and BuzdirJ .l2000a; 
llshigurol l2000t [Nam et all Il999i ISingleton etaR 120011 120001: ISvmington et all l200l|) and, for a similar compound, 
llGoddard et all 12002ft . These studies indicate that some evidence of the formation of the LOFF state has been 
reached. For example ijSingleton et all l2000|) studied resistance and magnetic behavior of single crystals of this 
superconductor in magnetic fields up to 33 T and at temperatures between 0.5 K and 11 K. When the magnetic field 
lies precisely in the Q2D planes of the material, they find evidence for a phase transition from t he superconducting 
mixe d state into a LOFF state, manifested as a change in the rigidity of the vortex system. (Ma nalo and Kleirl 

2000) compare the theoretical anisotropic upper critical field H c of a quasi-two-dimensional d-wave superconductor 
with recent H C 2 data for r-(BEDT-TTF)2Cu(NCS)2 and find agreement both with regard to the angular and the 
temperature dependence of H c . According to these authors this supports the suggestion that the LOFF phase exists 
in thi s ma t erial fo r exactly plane-parallel orientation of the magnetic field. 

In Jul l20f)l field induced superconductivity was reported in an organic superconductor A— (BETS^FeC^ 
(BETS=bis(ethylenedithio)te traselenafulvalene). A pos sible mechanism to create field induced superconductivity 
is the Peter- Jaccarino effect llJaccarino and Peteii 11962ft . However the upwards convex nature of th e lower critical 
field as a function of the temperature casts doubts on this interpretation. Therefore some authors, e.g. ijBalicas et all 

2001) , have proposed tha t these results can b e int erpreted as eviden ce of the formation of the LOFF state. These 
results were reviewed in l|Houzet et all 12002^) and l)ShimaharaL 12002]) . In particular in the latter paper, an experi- 
mental phase diagram of the field-induced superconductivity in this organic c ompound was t heore tically reproduced 
by a combination of the LOFF state and the Jaccarino-Peter mechanism. ijTanatar et all |2002|) discusses wether 
LOFF state has been observed via thermal conductivity k(H) in quasi-two-dimensional organic superconductor A- 
(BETS)2GaCl4. For clean samples the behavior of k(H) is similar to the one expected by a second order phase 
transition and is consistent with the formation of a LOFF phase. 



E. Future developments 

The superconducting LOFF phase might be realized even if the difference in chemical potentials of two species 
were not generated by a magnetic field acting on electron spins. Apart from nuclear physics and pulsars, to be 
discussed below, another context might be offered by ultracold quantum degenerate Fermi gas of atoms comprising 
two hyperfine states. The experimental investigations of ultracold gases were first ded icated to the study of the Bose- 
Einstein condensation ijAndersorl Il995t iBradlev et al. L ll997t Davi sl ll995l:lFriedl. l 1998|) : subsequently these techniques 
have been extended also to magnetically trapped u l tracold alkali Fermi gases or to gases wit h coexisting Bose-Einstein 



condensate and Fermi gas ( Modug no et all l2002t iRoati et all l2002t ISchreck etal. , 2001 ) . In particular two s tate 
mixt ures of ultracold gases have be en employed, with 40 K vapors llDeMarco and Jir. L 199E : iDeMarco e^"aD . l200l|) . or 
6 Li ()Granadel 12003 lO'Haral l200l|) . or a mixture of 6 Li and 7 Li l)Mewes et all ll999j )~~Future developments could 
lea d to the observatio n of superconductivity and Cooper fermion pairs condensation in these systems. As discussed 
in ijCombescotl l2000j) it is quite likely that the two hyperfine states would have different atomic populations, since 
at the moment there are no known fast relaxation mechanisms to equalize the two atomic population s. Therefore 
super conductivity for two-state ultracold Fermi gases is likely to be of the LOFF type. The author of Ref. l)Combescotl 
2000) has performed a theoretical study of 6 Li under the above mentioned conditions; he considers not only s— wave 
interactions, but also an anisotropic term induced by density fluctuation exchange and shows that the range where 
the LOFF phase is realized increases with the increasing role of the anisotropic term. This is an interesting theoretical 
development, which adds new interest to the experimental investigations of ultracold atomic Fermi gases. It remains 
to be seen, however, if such possibility is actually realized in Nature. 
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F. LOFF phase in nuclear physics 

Neutron-proton pair correlations and the possibility of n—p Cooper pair condensation are presently studied in several 
different contexts, from heavy ion collisions to astrophysics. They have been investigated, using the BCS theory, in 
infinite nuclear matter l|Alm et all 1 19931 1 199(1 ll99rA fealdo et all Il992t ISedrakian et all 119971 IVonderfecht et all 
119911) . and by mean-field effective interactions in finite nuclei. In several cases nuclear matter is highly asymmetric, 
with proton concentration at most 30-40% in supernova matter and 10% in neutron stars. These asymmetries are 
detrimental to nucleon superconductivity; on the other hand, for weakly asymmetric states, fermion condensation is 
indeed possible. For example, weakly isospin asymmetric nuclear matter favors the formation of Cooper pairs in the 
3 Si- 3 -Di channel, due to the presence of a tensor force; gaps are of the order of 10 MeV. Condensation in this channel 
might be relevant for low density bulk matter such as dilute nuclear matter in supernovas. On the other hand there 
is no evidence of larg e gap isospin singlet pairing in ordinary nucle i, which might be explained by the presence o f 
spin orbit interaction l)Goodmanlll999MMartinez-Pinedo et a/.L IT999]) . The authors ijSedrakian and LombardoL l2000() 
study the dependence of the gap as a function of both the isospin asymmetry a np = (p n — Pp)/ P and the temperature, 
using realistic nuclear interactions. For small asymmetries the gap develops a maximum at a certain intermediate 
temperature; for large asymmetries the superconducting phase exists only at finite temperature, because the smearing 
effect of the temperature on the Fermi surfaces favors condensation. At higher values of a np (— 0.11 in their model) 
pairing is no longer possible. 

Also in the context of isospin asymmet ric nuclear matter it is possible to have a transit i on fro m the BC S state to a 
LOFF phase instead of the normal state l|Isavevl l2002HSedrakiarl I200H ISedrakian et a0.ll997[) . In Ref. ijSedrakiarl 
l200l|) the possibility of spatially inhomogeneous condensate in asymmetric nuclear matter is studied. Condensation 
is possible in different channels. The isospin triplet channels are favored for large enough asymmetries; more exactly 
the channel 1 Sq dominates at low densities and the channel 3 P2~ 3 F 2 (or 1 P 2 ) a t high densities. For weak asymmetries 
the dominant channels are the isospin singlets 3 Si- 3 Di (low densities) and 3 D 2 (high density). The author considers 
the case of low density; as the isospin singlet 3 Si- 3 Di has a strength much larger than the isospin triplet 1 So, he 
neglects the latter. The interaction is modelled by the Paris nucleon-nucleon potential. The gap equations are solved 
numerically and have non trivial solutions for non vanishing total momentum of the pair P. The LOFF phase is 
favored for a np > 0.25 and P — Q.ipp, independent of a np . For a np > 0.37 pairing exists only in presence of non- 
vanishing P. The maximal values of a np and P compatible with the LOFF state are 0.41 and Q.3pf respectively (the 
actual values are indicative, as a refinement of the treatment of the nuclear interaction may change them by a factor 
as large as 3). The results are obtained at T — 3 MeV. From the BCS to the LOFF phase the phase transition is first 
order, while one passes from the LOFF to the normal state by a second order phase transition. No attempt is made 
to determine the most favorite crystalline structure. 

Under hypotheses similar to those of the previous paper ijlsavevl [2002) studies the effect of coupling between 
the isospin sing l et an d isospin triplet, since at low densities pairing between these two channels may be important 
(Akh iezer et adll999j) . Besides, the autho r goes beyond the app roximation of "bare" nucleon interaction, by using the 
Fermi-liquid phenomenological approach ijAkhiezer et ^Zl. ll994|) . By these chan ges one finds interesting peculiarities 
at T — 0. First, the triplet-singlet channel turns out to be energetically favored; second, the phase transition from the 
LOFF to the normal state can be of first order, depending on the nature of the nucleon interaction. While still model 
dependent, these investigations of the LOFF phase in nuclear interactions are interesting as they offer, in principle, 
a different way to the LOFF phase. To be closer to phenomenology one should consider however more complicated 
structures such as, for example, hyperon rich matter. Alternatively the modulation of the order parameter might be 
caused by Pauli paramagnetism due to strong magnetic fields in highly magnetized neutron stars ( magnetars) . In thi s 
case one could have a splitting in the Fermi surfaces of a nucleon pair in the / = 1, L = channel IjSedrakiari 12001 1 . 



G. Why color LOFF superconductivity could exist in pulsars 

In this Subsection and in the next one we will be interested in some numerical estimates of the values of the 
parameters needed for the LOFF phase in pulsars to occur. In general, color superconductivity in quark matter might 
be realized in compact stars. This expectation follows from the following two facts. 

First of all the BCS critical temperature is given by 

T c = 0.57A BC s (6.13) 
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high chemical potentials fi 3> 10 8 MeV (iRaiaeop al and Shusterl bOQOl). and from mode l s with parameters ad- 
justed to reproduce the physics at zero densit i es iAlfordl 120001: lAlford et all Il998l fl"999bt iBer ees and Raiaeopa l 
19991: ICarter and Diakonovi Il99flt lEvans et all Il999albt lR.ai agonal Il999t iRaop et all 119981 1200(1 ISchaferl 12001: 
Schafer and Wilczekl . 11999a IWilczekl l2000|) . None of these calculations is valid at chemical potentials around 400 



MeV, which correspond roughly to the density of the inner core of a neutron star, as we shall see below. However in 
all these cases one gets values of the gaps of the order quoted before. 

The second fact has to do with the thermal history of a pulsar. The general belief is that compact stars such 
as pulsars are formed in the core of a supernova explosion. The temperature at the interior of the supernova is 
about 10 11 K, corresponding to 10 MeV (1 MeV = 1.1065xl0 10 K). The star cools very rapidly by neutrino emission 
with the temperature going down to 10 9 -10 10 K in about one day. The neutrino emission is then dominating the 
cooling for one thousand years. In this period the temperature reaches about 10 6 K. After this period the star cools 
down due to X-ray and photon emission and in a few million years reaches a surface temperature around 10 5 K. 
Therefore, for the greatest part of its existence a neutron star has a temperature lower than the critical temperature, 
with the possibility of forming color superconducting condensates. It follows that also in this context a compact star 
can be considered at zero temperature because its temperature is much smaller than the typical BCS energy gap, 
T n . s ./A BCS « 10- 6 - HT 7 . 

We have seen previously that QCD favors the formation of BCS condensates in idealized cases, e.g. two or three 
massless flavors of quarks. However in realistic cases the three quarks have different Fermi momenta due to the mass 
difference. It is interesting to have an idea of the order of magnitude of the scales involved in the description of a 
neutron star with a quark core. We begin with a very crude example of a free gas of three flavor quarks, taking up 
and down massless and the strange one with mass M s ijAlford et all |2000() . Requiring that the weak interactions are 
in equilibrium it is easy to determine the chemical potentials and the Fermi momenta for the quarks. We find 

_ 2 

1 d 
fi d = /i + -/i e , p F = (i d , 

M + ^Me, P F = Vrf - Ml , (6.14) 



Ms 

where fi is average chemical potential 

A* = ^ (fhi + fJ'd + AO (6.15) 
and fi e the chemical potential of the electrons. Notice that 

HiNi + (j, e N e = fiN q - n e Q, (6.16) 

i— u,d,s 

where 

N q = V N u Q = lN u -\(N d + N s )-N e . (6.17) 
,i 6 6 

■i—u,d,s 

The chemical potential for the electrons is fixed by requiring electrical neutrality, corresponding to the following 
condition for the grand potential f2 at zero temperature 

Q = ^- = 0. (6.18) 

il is obtained from Eq. Q2.6()|l (omitting the volume factor) 



1 



PF 



In our case we get 



n = — / p 2 (E(p)~fi)dp. (6.19) 
n Jo 



V = ^ 1^ P 2 {E l {p)-fi l )dp+\ [' p 2 (p~fi e )dp, (6.20) 

K j Jo ft Jo 
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where 

K,d(p)=P, E s (p) = V^TMl . (6.21) 

Although the integral is feasible its expression is algebraically involved and it is easier to do all calculations numerically. 
In particular the result for the chemical potential of the electrons for different values of /i as a function of M s is given 
in Fig. |21 




FIG. 24 The chemical potential of the electrons vs. M a for three values of the average chemical potential. 



We can get an analytical expression by performing an expansion up to the order M^//i 4 . One gets 
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The baryon density is obtained as 
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(6.22) 
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(6.24) 



The plot of the ratio of the baryon density to the nuclear baryon density is given in Fig. |23 The nuclear baryon density 
has been assumed as the inverse of the volume of a sphere of radius about 1.2 fm. Within the same approximation as 
before one finds 



PB ~ 2 



1 (Ms 



(6.25) 



We note that densities in the core are of the order of 10 15 gj cm 3 , corresponding to a chemical potential of the order 
of 400 MeV, as shown in Fig. ESI 



In particular let us discuss the range of values around 400 MeV of the average chemical potential, with a strange mass 
of the order 200-300 MeV (the strange mass here is not the current mass but an effective density dependent mass). 
With M s = 300 MeV one finds /i e = 53 MeV (56 MeV from the approximate equation) with Fermi momenta 



p u F = 365 MeV, p d F = 418 MeV, p s F = 290 MeV , 



(6.26) 
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1*1= 300 MeV 
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FIG. 25 The ratio of the baryon density of the free quark gas to the nuclear baryon density vs. Ms, for three values of the 
average chemical potential. 



and a baryon density about 4.4 times the nuclear matter density. With M s = 200 MeV the result is fi e = 24 MeV (25 
MeV from the approximate equation) and 

p u F = 384 MeV, p% = 408 MeV, p s F = 357 MeV , (6.27) 

and a baryon density about 5.1 times the nuclear matter density. To go to baryon densities relevant to the central 
core of the star, i.e. densities from 6 to 8 times the nuclear matter density, one needs to go to higher values of \i and 
lower values of M s where the difference among of the Fermi momenta is lower. This can be seen from Fig. 1261 or 
using our approximate expression for /i e : 

l4«M+^, 1*«M-^, (6-28) 



with 



p F -p F ttp F -p F n^l. (6.29) 



The previous results are rather general, but in order to discuss the possible astrophysical applications we need to fix 
a value for A . Notice that we can trade the coupling constant G ( see Eq. (14.551) ) for A n s i nce G is fixed once we 
give the cutoff 5. On the other hand, the equation for the chiral gap ijRaiagopal and Wilczekll200l|) gives a relatione 
between the NJL cutoff A = 5 + n, the coupling G and the constituent quark mass. By taking the constituent mass 
around 300 — 400 MeV and fixing A, one has still a parameter to play around and it is possible to get values of Ao 
from about 20 MeV up to about 100 MeV. In the present case, since the typical value of Sfj, inside the LOFF window 
is 0.7Ao and 

fy= -^(Vd-Vu) = ^Me, (6-30) 

we can reproduce approximately the situation illustrated at the beginning of this Subsection with M s = 300 MeV by 
choosing A = 40 MeV. With this choice the LOFF grand potential at Sfi = 5fii is of the order 10~ 7 GeV 4 which 
as we shall see in Section lYl.HI is of the right order of magnitude to give rise to the glitch phenomena ijAlford et all 
2000). Notice also that the LOFF condensate evaluated at S/ii 

AiOFF^Mi) « 0.25A = lOMeV (6.31) 

is much larger that a typical temperature of neutron stars, of the order of keVs. 
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FIG. 26 The Fermi momenta of the three quarks vs. Ms. 

H. Astrophysical implications of the QCD LOFF phase 

While a great experimental effort is devoted to the search of the LOFF phase in condensed matter, so far nothing 
similar happens for the crystalline phase of QCD. The reason is that it is difficult to create the experimental conditions 
of high density and low temperature for hadronic matter. The crystalline superconducting phase of quarks may 
however result relevant for astrophysical dense systems, in particular in the explanation of pulsar glitches. Pulsars are 
rapidly rotating stellar objects, characterized by the presence of strong magnetic fields and by an almost continuous 
conversion of rotational energy into electromagnetic radiation. The rotation periods can vary in the range 10 -3 sec up 
to a few seconds; these periods increase slowly and never decrease except for occasional glitches, i.e. sudden increases 
of the rotational frequencies, when the pulsar spins up with a variation in frequency of the order of 5CI/Q, « 1CP 6 or 
smaller. Glitches are a typical phenomenon of the pulsars, since probably all the pulsars experience them. 

Pulsars are commonly identified with neutron stars; these compact stars are characterized by a complex structure 
comprising a core, an intermediate region with superfluid neutrons and a metallic crust. With a chemical potential of 
the order of 400 MeV, as we have seen, the conditions for color superconductivity in the CFL or the LOFF versions 
might be reached in the core. Before examini ng this possibility, let us however describe the standard explanation of 
glitches, in the form originated by the papers ijAlpar et all 1 1984at I Anderson and ItohUl973l) . This model is based on 
the idea that the sudden jumps of f2 are due to the movements outwards of rotational vortices in the neutron superfluid 
and their interaction with the crust. Crucial ingredients of the model are therefore the existence of a superfluid and 
a crystal (the metallic crust). This is one of the main reasons that allows the identification of pulsars with neutron 
stars, as only neutron stars are supposed to have a metallic crust. The LOFF state can be relevant in this context 
because, if there is a LOFF phase inside the pulsar, the superfluid might interact with the LOFF crystal instead of 
the crust, thus providing an alternative or complementary mechanism for the glitches. Thus far, there is no developed 
model for the pinning of the superfluid vortices to the QCD LO FF crystals within compact stars. Therefore we lim it 
our survey to an introduction to the subject, along the lines of l|Alforol EoOOt lAlford et a/1 l2001rilNarduHill2002c|) . 

Let us consider a compact star whose metallic crust rotates with angular velocity ft. The superfluid inside the 
star should not rotate because, in absence of friction, the crust cannot communicate its rotation to the superfluid 
component. The velocity of the superfluid is v s = h V^/m where $ is the phase of the superfluid condensate wave 
function. The consequence of this formula would be § v s ■ d£ = 0. This would imply the absence of rotation in the 
super fluid, which however does not correspond to the state of minimal energy (for a discussion see ijLandau et all 
1980)). The correct condition is 



where k is the quantum of vorticity: k — h/m. For Eq. I|6.32|) to hold the curve 7 must wind a singular point; the 
integer n is the winding number which counts the number of times the curve goes around the singular point; the most 
energetically favorable condition is realized by n — 1. If 7 is in a plane the condition H(j.32|l holds for any plane and 
the locus of the singular points is a vortex line (v.l.). In absence of rotation there are no v.l.'s; the minimal angular 




7 



(6.32) 
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velocity f2 m j n for the formation of the first vortex line is 

Qcrit = — %^ In — . (6.33) 

Here we are assuming a cylindrical configuration with radius R; a is a cutoff of the order of the interatomic distances. 
By increasing f2 also the number N of vortex lines per unit area in the superfluid increases according to the formula: 

N = ^ (6.34) 

and one gets, instead of l|6.32[l . 

<f d£-v s = NA2ttk , (6.35) 

where A is the area encircled by 7. Eventually the v.l.'s tend to fill in all the space. As a numerical example 
one can estimate iV for the pulsar in the Crab nebula. Here m = 2mjv (the condensate is formed by neutral 
bosons: pairs of neutrons) and fl = fl pu isar gives N ~ 1.9 x 10 5 cm~ 2 with an average distance between vortex lines 
d ~ N^ 1 / 2 ~ 10~ 2 cto. If the vortex line is a straight line, v s is perpendicular both to the vortex line and to the 
radius joining the singular point and the point at which we compute v s . At a distance r from the singular point one 
has 

v g = — , 6.36 
r 

as can be immediately seen from (|6.32|) . More generally: 

K f d£AR 



2 J,, * ■ ( " 7) 

where R is the distance vector from the vortex line to the point at which we compute the superfluid velocity. 

During the rotation the vortex lines follow the rotational motion of the vessel, which is clear because they are pinned 
at the boundary of the superfluid; in particular, for rotations around an axis, the vortex lines are, by symmetry, straight 
lines parallel to the rotation axis. Their motion imitate the motion of the liquid as a whole and, as a consequence, 
also for the superfluid one can use the hydrodynamical formula 

r? = ivAv s , (6.38) 

which in principle would be valid only for the fluid normal component. 

Let now v{r) be the number of vortices per unit area at a distance r from the rotation axis; if v = v s is the 
superfluid velocity, one gets 

<fd£ v=[ dS ■ V A v = 2ttk f 2vr' v(r')dr' . (6.39) 

Jj Jo Jq 

We put k = 2itk = h/2m n and write l|6.39|) as follows: 

2irr 2 n{r) = k [ 2-kt' v{r')dr' , (6.40) 
Jo 



which implies 



kv(r) = 2Q(r) + r — . (6.41) 
or 



Since the total number of v.l.'s is conserved, one has 

dv/dt + V • (i/vy) = , (6.42) 
where v r is the radial component of the superfluid velocity. We write (|6.4U|) as 



27iVQ(r) = k / vdS (6.43) 
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and take the time derivative, using (|6~4l|) to get 

2nr 2 — = -k dSV-{vw r ). (6.44) 
ot J 

2 dn 

Using the Gauss theorem one gets 2nr — — = —k2irrvv r , i.e. 

dQ, kvv r x dCl\ v r .„ , „. 

*=-— = - { m{r) + r ^)v ^ 

Eq. (|6.45|) shows that the only possibility for the superfluid to change its angular velocity (Cl ^ 0) is by means of a 
radial motion, i.e. v r =/= 0. 

Let us now conside r a rota ting superfluid in contact with rotating normal matter on which an external torque 
is acting l)Alpar et all Il984a(l . We denote by I c and fl c the moment of inertia and angular velocity of the normal 
components that, in a neutron star, includes the crust and possibly other normal components. The equation of motion 
of the normal component is 

I c (l c (t) = M^t + M in t . (6.46) 

Besides the external torque M ext , basically related to the spin down of the pulsar (or the steady accretion in binary 
pulsars), we have included the internal torque Mint due to the interaction with the superfluid: 

M int = - J dl p n(r,t) (6.47) 

where dl p is the infinitesimal moment of inertia of the superfluid component. Eqs. H6.45I6.47|) are the equations of 
motion for the angular velocities £1 and tt c (superfluid and crust). The two velocities are coupled not only through 
Mi n t, but also by v r , because we will show below that v r depends on the difference £1 — Q c . We note again that 
fundamental for this model is the existence of radial motion, for, if v r = 0, then = const, and only f2 c changes, due 
to the external torque alone. 

In the neutron star, superfluid neutrons (in Cooper pairs) coexist with nuclei of the crust. Also in the crust there 
are superfluid neutrons, but they are characterized by a different (and smaller) A. Computing the difference in the 
free energies between the two phases one obtains the difference of pressures and, consequently, the force per unit 
length of vortex line. Let b be the average distance between the nuclei; b is also the average distance between two 
consecutive pinning centers. Let us assume 

2£o < b , (6.48) 

where £o is the superconducting coherence length, which also gives the dimension of the vortex core, since £o is of the 
order of the spatial extension of the Cooper pair. The maximum pinning force is obtained when the vortex passes 
through one layer of the lattice; therefore the maximum force per unit length of vortex line is 

f P - ^ , (6-49) 

where 

8E P = F s -F c (x pA* , (6.50) 

where F s and F c are the free energies of the superfluid neutrons and the nucleons in the crust; A s is the gap for 
superfluid neutrons and one can neglect A c , the gap of superfluid neutrons in the crust since A c <C A s . Eq. (|6.5U|) 
implies that neutrons tend to remain in the volume V of the vortex core because they experience a force repelling 
them from the superconducting phase (if neutron rich nuclei are present, the repulsion will be less important). Typical 
values for the pinning energy per nucleus SE P at densities 3 x 10 13 — 1.2 x 10 14 g/cm 3 are 

SE p = 1 - 3MeV , (6.51) 

while b — 25 — 50 fm and £o = 4 — 20 fm give 

/ p = 40 - 1200 MeV 3 . (6.52) 
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On the basis of these considerations let us now ske tch a possible mechanism for the formation of glitches ijAlparl 
119771 lAfoar' et al .1 Il984albl lAnderso n and ItohL Il973|) (for further references see below) . We consider the rotating 
neutron star with superfluid neutrons in its interior and a metallic crust, which is a simplified model, but adequate 
to our purposes. As stressed already, we distinguish between the superfluid velocity Q and the crust velocity f2 c . Let 
us suppose that they are initially equal, which is a consequence of the pinning. Due to the spinning down of the 
star, f2 c decreases; as long the vortex cores are pinned to the crust lattice, the neutron superfluid cannot spin down, 
because the radial motion is forbidden. There is therefore a relative velocity of the superfluid with respect to the 
pinned vortex core because £1 > fl c : 

5v = (O - fi c ) A r . (6.53) 

The interaction between the normal matter in the core of the v.l. and the rest of normal matter (nuclei in the lattice, 
electrons, etc.) produces a Magnus force per unit length given by 

f = pkA<5v, (6.54) 

where the direction of k coincides with the rotation axis and its modulus is equal to the quantum of vorticity. / is 
the force exerted on the vortex line; as it cannot be larger than f p there is a maximum difference of angular velocity 
that the system can maintain: 

^ = ^-^* = ^ = ^- (6-55) 

If bj < u> cr the vortices remain pinned at the lattice sites instead of flowing with the superfluid as they generally do 
inside it (sec discussion above). On the contrary, if lu > co cri the hydrodynamical forces arising from the mismatch 
between the two angular velocities ultimately break the crust and produce the conditions for the glitch. A possible way 
to get it is by the observation following Eq. I|6.45|l . If a bunch of vortex lines are unpinned and move outwards then 
cq. H6.45fl implies that the angular velocity (and the angular momentum) of the superfluid decreases, and, therefore, 
the angular momentum of the crust increases, which is revealed from outside as a spin up of the star, i.e. a glitch. A 
numerical analysis would imply solving the set of Eqs. I|6.45I6.47|I . but this is outside the scope of the present review 
5 . Let us instead discuss the possible role of the LOFF phase in this context. The QCD LOFF phase provides a 
lattice structure independently of the crust. Therefore it meets one of the two requirements of the model for glitches 
in pulsars we have outlined above, the other being the presence of a superfluid. On the other hand the only existing 
calculations for the inhomogeneous phase in color superconductivity have been performed for the case of two flavors, 
which however, in the homogeneous case, does not present superfluidity, since there are no broken global symmetries. 
Superfluidity is on the other hand manifested by the CFL phase of QCD. Therefore a realistic application to QCD 
superfluid has to wait until a calculation of the LO FF phase with three flavors will be completed. For the time being 
one can give some order of magnitude estimates llAlford et oil . l200ibD . Let us assume the following choice of the 
parameters: A 2 sc = 40 MeV, Alqff ~ 8 MeV, corresponding to the Fulde-Ferrel state; since q w 1.2£/z « 0.7A 2 sc, 
one would get for the average distance between nodal planes b = 7r/(2|q|) f=s 9 fm and for the superconducting 
coherence length £o = 6 fm. From l|3.28[l . with 6fj, — 6fj,% and an extra factor of 4 to take into account the two flavors 
and the two colors, we get the free energy per volume unit as follows 6 : \Floff\ = 8 x (10 MeV) 4 and therefore, from 
(|6.50() . the pinning energy of the vortex line is 

SE P = \F LO ff\ x b 3 = 6 MeV . (6.56) 

To get the pinning force we cannot use l|6.49|l since l|6.48|) does not hold in this case. For an order of magnitude 
estimate one can use 

U - ^ , (6-57) 
giving a pinning force per unit length of the vortex of the order of 

f p « 3 x 10 3 MeV 3 . (6.58) 



Models differ in the mechanism by which angular mome ntum is released; instead of pe r formi ng outward movement s for example v.l.'s 
might break the cr ust or rearrange it. For reviews see lAlparl Il995t iPines and Alpa.it Il985ft and, more recently, iAIpar et all Il993t 
lEpstein and Bavmt ITiffll iLink and Epsteirlll99ftlRudermarlll99lllRuderman et aZ.LIl998ft . 

Using the exact expression instead of 13.281 . that is valid only in the weak coupling limit, one would get \F^Qpp\ = 5 X (10 MeV) 4 . 
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Comparing these numerical values with Eqs. I|6.51|l and (|6.52|l one can see that these order of magnitude estimates 
give figures similar and therefore some of the glitches in neutron stars may be generated well inside the star by vortices 
related to the LOFF phase of QCD. 

As we already stated these conclusions are tentative because the analysis of the QCD LOFF phase needs extension 
to the three flavor case; moreover the true ground state is likely to be different from the Fulde-Ferrel one plane wave 
structure. Nevertheless they are encouraging and leave open the possibility that neutron stars might give another 
laboratory where to study the inhomogeneous superconducting phase. It can be useful to stress that even in quark 
stars, in the QCD superconducting LOFF phase, one would get a crystal structure given by a lattice characterized 
by a geometric array where the gap parameter varies periodically. This would overcome the objection that pulsars 
cannot be strange stars. This objection is based on the following observation: If strange matter there exists, quark 
stars should be rather common; however, in absence of metallic crusts, strange stars can hardly develop vortices, 
at least by the model we have described here. On the contrary, if the color superconductivity is able to produce a 
crystalline structure it could also give rise to glitches and the argument in favor of the existence of strange stars would 
be reinforced. 



VII. CONCLUSIONS 

Inhomogeneous crystalline superconductivity was predicted forty years ago by Larkin, Ovchinnikov, Fulde and 
Ferrell, but realistic conditions for its experimental investigations became available only a few years ago. In condensed 
matter the existence of the LOFF phase, with its characteristic space modulation of the energy gap, still awaits 
complete confirmation. This is due to the fact that it is indeed a subtle effect. It arises when the Fermi surfaces of the 
two species participating in the Cooper pairing are different. However for large separation pairing is not possible at 
all and superconductivity disappears altogether. In condensed matter the separation of the Fermi surfaces is obtained 
by a Zeeman splitting due to an exchange interaction due to a magnetic field. However the needed field strength 
is such to destroy superconductivity due to diamagnetic effects. As we discussed in the paper, the way to avoid 
the problem was to use unconventional superconductors such as organic compounds. These materials have in fact a 
layered structure and therefore, if the magnetic field is parallel to the layers, the orbital effects can be controlled. 

New opportunities have recently arisen to detect the LOFF phase in atomic physics (ultracold atomic gases) , nuclear 
physics and especially quark matter. This last development is a consequence of the recent excitement generated by the 
study of QCD at high density and small temperature. Inhomogeneous crystalline superconductivity in this context 
could be generated by the difference in quark chemical potentials induced by weak interactions in the inner core of 
pulsars. Their main phenomenological effect might therefore be to provide a mechanism for the explanation of glitches 
in pulsars. If pulsars are neutron stars with a core made up by color superconducting matter, this mechanism would 
be complementary to the standard models of glitches. If pulsars are strange stars, then the crystalline structure of 
the condensate would provide the possibility for pinning the superfluid vortices and eventually creating the glitches. 

This paper was mainly aimed to the presentation of a unified formalism to describe the LOFF phase both in 
condensed and hadronic matter. The simplest way in our opinion to describe superconductivity effects, including the 
LOFF state, is by the effective lagrangian approach. Since they are based on the general mathematical ground of the 
Renormalization Group, effective lagrangians allow the conditions for this unification. The existence of a common 
mathematical basis should allow experts of one side to fully appreciate and take advantage of the progresses made in 
the other. We would be gratified if this paper turned out to be useful to this end. 
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APPENDIX A: Calculation of J and K 

We give here an outline of the calculation of the integrals J and K appearing in the GL expansion at T = 0. Using 
the definition of J, Eq. (|3.5|) . and K, eq. (|3.6[1 . we have 

J = J(qx,q 2) q3,q4)=+^f | ^ | + °° g f[ ^M) 

igp /" dw / >+ ' 5 /•+ 00 d^A r 1 



' 2 J 4vr y_ 5 ^7_ 00 2tt W \~E + iesignE + £ - 8^ + 2w -ki 
E + iesignE - £ - tf/x - 2 w£i J ' ^ 



x = -ftT(qi,q2,q3,q4,q5,q6) 



i<7P /" dw 



47T 



d£ 



-5 



igp /" dw 

47T 



-5 
1 



+ iesignE — £ — 5a — 2 w-£; 



i=l 

;}■ 



1 



E + iesignE + £ — <5/i + 2 w • ki 



where we have defined 



ki = 0, k 2 = qi-q 2 , k 3 = qi q 2 + q 3 q4 , 

t-i = qi, = qi -q 2 + q3, 4 = qi-q2 + q3-q4 + qs, 



(A2) 



(A3) 



with the conditions qi — q 2 + q3 — q4 = and qi — q 2 + q 3 q 4 + q 5 — q 6 = for J and K respectively. We 
introduce the Feynman variables Xj, yj (j = 1,2 for J and j = 1,2,3 for K) to form the vectors k = J^Xiki and 
£ = yi£i', after rotation of the energy integration contour E —f ip^ we get: 



dE y-r 1 

2^ -1.1 J B + £-<5 A t + 2w-k i 
1 



«(i-£a*) 



2tt 11 £ - i - 5/x - 2 w£i 



2tt [zp 4 +£-(5 ) u + 2w-k] 2 

idp 4 S(l-J2Vk) 
2tt [ip 4 -£-^-2w-€] 2 



n 



next we perform the £ integration by the residues method and the angular integration; for J the result is 



J = — 



dp 4 e(p 4 ) (Sfi - ipi) 



Y[ dx n dy n 



,71=1 



g(i-E^)^(i-Si/fc) 

[(^-ip^a-^lk-^aja 



After the energy integration we remain with 

" 2 

dx n dy n 



71=1 



*(i-£^Wi-£ 



, Vk 



v 2 F \k-£\ 2 -5a 2 



This expression is general; we can specialize it to the various crystal structures, as explained in the text. 
For K we get, instead of I|A5J| . the result 



K 



Zigp 



dpie{pi) 



Y\ dx n dyn 



(l-5»*(l-$» 



x [5fjL — ipi) 



(5a-ip 4 ) 2 +v F \k-£\ 2 
[{5^-ip A ) 2 -v 2 F \k-i\ 2 } 



which, after energy integration, becomes 

3 

U> 



K=- 9 -l 



1 { dx n dy n 



(i-E-^^-E 



Vk 



5a 2 + 3v 2 F \k-£\ 2 
[v 2 F \k-£\ 2 -Sfx 2 f 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 
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APPENDIX B: Expansion of II around the tricritical point 

Let us consider the expansion of n(g) in Q = qvp, at finite T and /i, which can be obtained from Eq. I|3.4|) after 
introducing the Matsubara frequencies: 

1 f d\v ~t~oc ^ 

m = - 2 9pTj -^j_ s d£ £ — -J^^^^— -g-^y (Bl) 



Expanding the first denominator in the momentum q we find 



where, as in Eq. (|2.70() . 



(2w • quf)^ 



(B2) 



uj n = uj n + iS/i . (B3) 

Notice that we have inverted the sum over the Matsubara frequencies with the integration over £. In this way, as 
we did for the homogeneous case, we are converting the divergence in £ into a divergence in the series, which can be 
treated as before by introducing a cutoff in the sum. Performing the angular integration and the integration over £ 
with the help of the following integral 



1 1 7T 

^u;l+e (icon - 2m = { ~ ir 2^Q 2 n m+1 ' (B4) 



we get eventually 

+oo oo (_-\\m n2m 
n— — oo 

By using the definition of the first term in the grand potential as 1 — n(q) multiplied by 2/g we recover easily the 
expression i|3.42|l for a. In analogous way, to get f3 and 7 one proceeds expanding J (see Eq. (|3.5|l and K (see Eq. 

(El))- 
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